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On Liouville’s Function 


By R. Sherman Lehman 


1. Introduction. Liouville’s function A(n) is defined by the equation A(n) = 
(—1)’ where r is the number of prime factors of n, multiple factors being counted 
according to their multiplicity. Pélya [6] conjectured that the function 

L(z) = © Xn) 

nsz 
is negative or zero for all z = 2, and in fact this is true within the range where 
this function has been previously calculated. Calculations connected with the 
present study show that L(x) < 0 for 2 < x < 10°. Nevertheless Haselgrove [3] 
has shown that the Pélya conjecture is false and that there are infinitely many 
integers x for which L(z) > 0. However, his method does not furnish explicitly an 
x for which the conjecture fails; and in fact it does not give an upper bound for the 
first counterexample. In the present paper we shall describe calculations which lead 
to the result that L(906,180,359) = +1. We have not found a smaller value of 
x for which the conjecture fails, but also we have not proved that this is the smallest 
x greater than 2 for which L(z) is positive. 


2. Background and Heuristic Considerations. Liouville’s function is connected 
with the Riemann zeta function by the equation 


t(s) i a-1 n° 
Let the zeros of ¢(s) on the line Re s = 3 be p, = 4 + ty, (n = +1, 42,--- ) 
with y_, = —vy, and let yo = 0. If it is assumed that these zeros are all simple, 
then the function ¢(2s)/(s¢(s)) has simple poles at 4 + ty, forn = 0, +1, +2, --- 
with residues 


$(2s) _ 5 A(n) 


1 ¢(2p,) 
a= —; an = - ya. (n = +1, +2, +3,---). 
(3) pnt (pn) 

Fawaz [1] has obtained an explicit formula for L(x) which is valid if the Riemann 
hypothesis holds and the zeros of {(s) are simple. Under these assumptions he 
| showed that there is a sequence of numbers T, , with k < T, < k + 1, for which 
(1) L(z) =lm > a, x **™ + O(1) 

koo |7,/STe 
for x > 0. 

Let 

(?) A,r(u) - Ane 

lynlst 
Fawaz’s result suggests that one study numerically the behavior of Ar(u) for 
| different values of T. Since Ar(u) should be an approximation to e “L(e"), one 
might expect L(e") to be positive for a u for which Ar(u) > 0. 
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In §5 we shall show that if the Riemann hypothesis holds, if the zeros of [(s) 
are simple, and if a conjectured estimate for 1/¢’(p,) holds, then 





tw . 
(3) lim ‘| Kr(u — we “L(e") du — > Qn ene = 0, 
ere AC 1Y¥nlST 
where 
, _ sin Tt 
K,(t) = oa 


This suggests that the sum A;(u) will represent a smoothing of the func- 
tion ¢ “L(e") which effectively operates over an interval of width about 2x/T. 
However, since the kernel K,7(¢) is not always positive, the sum does not represent 
a true smoothing of e™L(e"). Thus, if we find a maximum of the sum A7(u) we 
cannot always expect that there will be a value of e ’“L(e") in the vicinity which 
is as high. Instead it is possible for airapid falling off of e “L(e") at a distance 
away of from x/T to 2x/T to be translated into an extra high peak of Ar(u). On 
the other hand, high values of A7(u) for several different choices of T will make 
such behavior appear less likely. 

Haselgrove’s disproof of the Pélya conjecture was based on a numerical study of 
the similar sum 


A;*(u) = - me, (1 i is ') eit 


lynls? 
Ingham [5] had shown that for any wo 


Ar*(uo) < lim sup e"L(e"). 
Hence to disprove the Pélya conjecture it was sufficient to find a 7 and wz for which 
Ar*(uo) > 0. Haselgrove found that Ato(831.847) = 0.00495. 

Using an IBM 701 at the University of California, Berkeley, we have inde- 
pendently computed approximations to the zeros of ¢(s) and the residues a, and 
have confirmed Haselgrove’s result. We have also obtained a smaller value of w for 
which A foo(t%) is positive. Our values to 4 decimal places are 


A tow(831.847) = 0.0050, A oo(831.847) = 0.0526 


Afow(814.492) 


0.0782, A jo00( 814.492) 


0.1102. 
As a result of a search for smaller values of u for which Ao0(u) > 0 we found 
Afooo(79.28) = —0.0418, — Aiow(79.28) = 0.0075. 


The number e””* is still a very large number, and there does not seem to be 
any more hope of calculating L(x) for « = e”™ than for z = e'™’. On the other 
hand, it is possible to find a method for calculating L(a) at isolated values which is 
quite feasible with present-day computers for x as large as 10°. Therefore, we com- 
puted Ar(u) for wu = 12.5(0.01)20.69 with T = 100 and T = 200, covering ap- 
proximately the range 2.7-10° < x < 9.7-10°. The vicinities of several high points 
were then selected for further study. The most promising of these was near 
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Fic. 1.—The functions A7(u) for 7 = 100, 200, 500, 1000. The maximum of Ayww(u) is 


near u = 20.623 = log (9.05-105). 


u = 20.62. Figure 1 shows a graph of A,r‘u) for u in this vicinity with four different 
values of 7’. Although none of the functions Ar(u) graphed is positive there, the 
strong upward trend suggested the ccnjecture that L(x) is positive for some x 


near 9.05- 10°. 


3. A Formula for Calculating L(x). A direct calculation of L(x) forall x < 10° 


by factoring each number would require far too much machine time to be feasible 
on the IBM 704—at least 1000 hours. Fortunately it is possible to devise an ap- 


propriate method for calculating L(x) at isolated values. 
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We begin with the known formula 
(4) > (2) =tva 


where by [y] we mean the greatest integer less than or equal to y. This formula is 
not suitable for calculating L(a2) for two reasons. First, there are [x] terms in it; 
and second, to calculate L(x) requires knowledge of values of the function L at 
large arguments such as x/2, 2/3, etc. The formula can, however, be modified to 
get around both of these difficulties. 

THEOREM 1. Let k, 1, and m be variables ranging over the positive integers. Let 
u(n) be the Mobius function. Let v, w, and x be positive real numbers with v < w < x. 
Then 


ua) = Bam /5]- Er (5) -[5))} 


pa A (i) 3 a, u(m). 


z z 
-<ls- 
. = 


Ly 


w 


Proof: Replacing x by x/m in (4) and breaking the sum into three parts, we 


obtain 
Vz] 5 (e)+.2,4(e)+2,4(2) 


<ns— 
re . “mm 


We multiply each of these sums by u(m) and sum for m S x/w, obtaining 


> wm) + L (=) - i (;) > u(m) = L(x), 
; ns 2 


p ok ra a ‘ L (=) = 2 L () > u(m), 


ms= a <ns= — <<i<- 4 
Dd u(m) - u(2)- u u(m) dD Dd Ak) 
"s, a<" Ss a ™s, mi" Sn "San 
= > u(m) D> Xk) Dd 1 
x x 
= i p(m) XL Mk) (Fa — |=}). 


Rearrangement of terms then yields the theorem. 
We observe that if v & 2'” and w & 2" 
portional to «”” if one has available a table of \(n) for n < w, a table of 


e(l) = 2 u(m) 


z 
ms— 
=w 


(5) 


for | < x/v, and a table of u(m) for m S x/w, 
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In practice it turns out to be more efficient to use a somewhat more complicated 
formula. Let k range over positive integers, and let I’, m’, and n’ range over positive 
odd integers. From (4) we obtain 


%4(5)- wva-[y/5), 


If x is replaced by x/m’ and the sum is treated as in the proof of Theorem 1, we 
obtain 


L(x) = & ulm’) 


o~8 
m s— 
= 


‘ ry oa xz deg j x i) 
© -{/2]-[/s]+[s 3|E - Daw [52,- 3) 


- © 1) & wm, 


a , 
i 
w ? 





m\i 
ms 


Bin 


We remark that similar formulas can be obtained for 
M(z) = > u(n) 


by modifying the formula 


> u(t) - | for x = 1. 


nsz 


4. Numerical Computations. The computations described in this section were 
performed on an IBM 704 at the University of California, Berkeley. 

Formula (6) contains two parameters v and w which can be chosen to minimize 
computation time. We chose to fix z/w = 1000 and takev = (10-*r)"” in order to 
make the program suitable for values of x near 10°. Two preliminary programs 
were then used to compute tables of A(n) for n S 10° and £(1), defined by (5), 
for odd 1 < 10°. 

In 1955 using the ORDVAC at Aberdeen Proving Ground, W. G. Spohn and 
the author computed Liouville’s function for n < 802,000 and verified that the 
Pélya conjecture holds up to this limit. In the present computation, the same 
method was used to obtain a table of \(n) for n < 10°. If one is given a table of 
\(n) for n S N/2, then the following sieving process will allow the extension of 
the table to N. One begins by entering —1 as the value for each integer from N/2 
to N. One then considers in turn each of the primes p S +/N, and one runs through 
the multiples of p among the integers from N/2 to N. If n is such a multiple, one 
sets A(n) = —A(n/p) after erasing the value already recorded for n. When this 
is done for all multiples of primes < +/N, the table is complete. For a machine 
with enough storage space to hold the table of A(n) for n < N/2, this method is 
much more efficient than factoring each integer in succession. 

The table of \(n) was placed on magnetic tape with each value taking up one 
bit. This table was then summed to compute L(x) for x = 100(100)1,000,000. 
The values for « = 1000(1000)802,000 were compared with the ORDVAC com- 











316 R. SHERMAN LEHMAN 


putation and all were found to agree. Also a comparison for x = 100,000( 100,000)- 
600,000 was made with a computation of D. H. Lehmer. Finally as a further check, 
the formula (6) was later used to compute L(x) for « = 200,000( 200,000 ) 1,000,000. 
The values found agreed with those obtained by summing the table of \(n). (The 
circularity here is only apparent; in these cases the formula (6) makes no use of 
the table for \(n) beyond n = 1000). To compute the table for \(m) and sum it to 
obtain L(x) required approximately 30 minutes of machine time. 

A table of the function £(/) was also needed for odd I < 10°. By a combinatorial 
argument, which is easy but requires consideration of a number of cases, one can 
show that if 1 < 10°, then —7 < &(1) S 7. Hence each value requires just 4 bits of 
storage. Again the table was put on magnetic tape. The method for computation 
was quite straightforward. Each odd square-free number m < 1000 was considered, 
and for each of its multiples u(m) was added to a corresponding accumulator. A 
check of the accuracy of the computation was made by using the formula 


| ets P | v | 
Pe 8) = 200m) | ome + 9 
with « = 10° and I’ and m’ running over positive odd integers. The program for 
&(1) required approximately 20 minutes. 

Next a program for computing Z(a) by the formula (6) was constructed. The 
odd square-free numbers m’ < 1000 together with the values u(m’) were stored as 
constants. Newton’s method was used to compute [+/] with especial care taken to 
avoid error due to round-off when ~+/y is an integer. The tables of \(n) and (1) 
on magnetic tapes were used as inputs, and L(2/l’) was obtained by summing the 
table of A(n). To obtain L(x) for a value of x near 9-10° required approximately 
16 minutes. 

Table 1 contains values of L(x) computed in connection with the search for 
a positive value. The order of computation is indicated in the last column. The 
values of « were chosen partly by guess and partly by heuristic considerations based 
on (3). There seem to be two separate peaks which for the functions A 7(u) were 
smoothed into a single maximum. After we found a positive value on the seventh 
trial, it did not seem worthwhile to pursue an investigation of the other peak. 
Hence we do not know whether the maximum for the other peak is also positive. 

















TABLE 1 
x L(x) | Order of computation 
| 
903 000 000 ~952 | 3 
904 000 000 —1144 2 
905 000 000 —1902 | 
906 000 000 — 584 4 
906 170 000 — 230 | 10 
906 200 000 90 8 
906 300 000 648 | 9 
906 400 000 | 708 7 
906 470 000 | 226 | rl 
906 500 000 — 120 | 5 
907 000 000 | — 1920 | 6 
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Some other values of L(x) computed by means of the main program are listed 
in Table 2. The values for z < 10° are all confirmed by agreement with those ob- 
tained by directly summing the table of \(n). 

The values of L(x) for z = 10° and 4- 10° are confirmed by results of a hand com- 
putation of D. H. Lehmer using the fermula 


L(x) = uM (;) a ~ u(1) 1/7 ~ u (4) | 4/5 | 


with g chosen to minimize computation. In this computation the values for M (z/k’) 
were taken from a corrected version of von Sterneck’s tables of M(x) (see [7]) with 
the following values differing from those given by von Sterneck: M (444,444) = 
—37, M(10°) = 212, M(4-10°) = 192. The first two of these values were obtained 
by factoring all numbers < 10°. The value for 4-10° was obtained by making small 
adjustments of von Sterneck’s computation which are required because of errors 
in the earlier tables. 

After finding positive values of L(2) we next took up the question of deter- 
mining zeros of L(x). The results given in Table 1 indicated that such zeros must 
occur between 906,170,000 and 906,200,000 and between 906,470,000 and 
906,500,000. Consequently a program was constructed to factor a!l numbers in 
these ranges. This program required approximately one minute for each 1000 
numbers factored. As a byproduct of this computation we obtained a further check 
of the program for computing L(x) at isolated values. 

One of the results of this computation was a listing of all zeros of L(z) in the 
intervals from 906,170,000 to 906,200,000 and from 906,470,000 to 906,500,000. 
In all, 167 zeros of the function L(2) were found in these intervals. We list the 
first 16 occurring in the first interval: 

906180358, 906180362, 906180364, 906180366, 906180370, 
906180374, 906180376, 906180386, 906180388, 906180390, 
906180518, 906180520, 906180524, 906180534, 906180536, 
906180554. 
There are 34 zeros from 906,192,698 to 906,193,478 inclusive; 22 zeros from 906,- 
194,914 to 906,195,298; 19 zeros from 906,195,986 to 906,196,098; 15 zeros from 
906,477,702 to 906,477,936; 43 zeros from 906,486,640 to 906,487,288; and 18 zeros 
from 906,487,932 to 906,488,080. 








PaBLe 2 
x L(x) x L(x) 

200 000 — 294 10 400 000 — 394 
400 000 — 460 10 410 000 — 330 
600 000 — 802 10 420 000 — 384 
800 000 — 600 10 430 000 — 300 

1 000 000 — 530 10 440 000 — 292 
4 000 000 — 1098 10 450 000 — 522 
10 460 000 — 588 


453 200 000 —27088 
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The first value of x greater than 906,170,000 for which L(x) is positive was 
found to be 906,180,359. We, of course, are not able to say whether this is the 
smallest x greater than 2 for which L(2) is positive. To decide this question without 
the use of essentially new ideas might very well require an enormous amount of 
computation. 


5. Derivation of an Explicit Formula. In this section we give a derivation of 
equation (3), which was used heuristically in finding where L(x) is positive; this 
derivation proceeds from several unproved assumptions. We shall assume that the 
Riemann hypothesis holds and that the zeros p, = } + iy, (n = +1, +2,--- ) of 
¢(s) are all simple. In addition we shall assume that there is a real number v < | 
such that 


(7) = O(p,’) (n = +1, +2,---). 


1 
£"(pn) 
Numerical evidence makes this conjectured estimate appear quite plausible. The 
twelve largest values of | 1/¢’(p,) | for0 < y, < 1000 are listed in Table 3. 

As in §2 let Kr(t) = (sin Tt)/(at). 
Lemma 1. If R > 0 and T > 0 and y is a real number, then 


| ( : ) | f 
RL, Sa | | 
e “ j)1 +0 RT — 17) ifiy| <T, 
K,(t)e” dt = ) 
LR 1 “a , . 
(° (erst =n) veienitian 

Proof: We have 


“ 2 f* sin Tt 
[ K,(t)e' dt = = = - cos yt dt 
R T 0 





* sin (T — y)t 
- -— ri 
0 t 


1 f* sin (T + y)t , 
t 


a+! 
T 40 T 
a * 8i(R(T +)) + * si(R(T - )), 


where 
fs +0 () for x > 0, 
7 oa x 
Si(x) = [ aa} 
6 t ~, l : 
| -2 + O{- forx < 0. 
v4 Hy 


The conclusion of the lemma follows immediately. 
We shall also use the following estimate for R > 0 and T > 0: 


m 9 pRT |<: . a 
| | K(t)| dt = tt | | sin ¢ | dt < ] dt + dt 
LR T Jo t 4 i 


= 1+ log (RT). 





(' 


> 0, 


a} 
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TABLE 3 

s v. | 11/8 + fy) | 
1 14. 134725 | 1.2608 
2 21.022040 0.8796 
72 185 .598784 0.8109 
135 294 .965370 0.8029 
298 540 .213166 0.8357 
299 540 .631390 0.8892 
363 630 . 473887 0.8334 
364 630 .805781 0.9106 
436 728 .405482 0.8371 
437 728 . 758750 0.8491 
606 946 . 765842 0.9744 
607 947 .079183 0.9914 





We are now ready to derive (3). We shall assume throughout that R > 0, 
T > 0, and T # y, for n = 1, 2, 3, --- . The remainder of the notation is explained 
in §2. 

Fawaz [2, p. 284] has shown (assuming the Riemann hypothesis and the sim- 
plicity of the zeros) that if u is restricted to a finite interval, then there is a positive 
constant C independent of k and u such that 

‘by age T < 6. 

1Y¥nlSTe 
Also | Kr(u — w)| S T/x. Consequently we can apply the Lebesgue bounded 
convergence theorem to obtain 


ot 


R * 
lim (Kr(u—w) >> ane) du 


kwao Yw—R lynlSTe 


w+R 


= K-(u — @) (lim > On e'™") du. 
w—-R kow |%nISTe 


Equation (1) can be written in the form 


e™L(e") =lim D> ane™ + O(e). 


koo |YnlSTe 
Hence we obtain 
o+R - ” 
[Kew =e ML (€*) du = tim (ane [ Ko(t)e™ at) 
— koo |¥nlSTe —R 


(9) w+R 
+ / Kr(u — w) O(e™) du. 
w—R 


If T;, = 27, then by Lemma 1 we have 


R ° ° 
an” | K,(t)e"™ dt= >> on™ 
-R 


l\Y¥nlST 


} . | dn | > | an | ) 
' ' O , 
lias om R\T — |. :) ad Lz. R( | yn |/2) 


lY¥nlSTe 


(10) 
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Since ¢(1 + i) = O(log t) for t > 2, we obtain the following estimate from the 
assumption (7): 


| an | | ¢(2p,.) | v—2 
oon = Oy, log 7,) 
| Yn | | Y¥nPat (p,.) | . 
for n = 1, 2, 3, --- . Since the series d, vy," converges for 6 > 1, (see [4, p. 57]), 
n> 


it follows that the series 


! 
| on| 
Yn>0 Yn 


converges; hence the last term on the right-hand side of (10) is O(1/R) uniformly 
in k. For fixed T the second term on the right-hand side of (10) is also O(1/R) 
since it is a finite series. 

Estimating the last term on the right-hand side of (9), we obtain 


otR R 
/ K+r(u — w) O(e™) du = O(e4#"™) i] [Ky (t) | dt 
7) —R 


—R 
= O(e*™ log (RT)) 


for R > 2/T. 
Combining (9) and (10) and letting k — ~, we conclude that for fixed T’, 


w+R 
/ Kr(u—w)e"L(e")du= > e™ +0 (; r) + 0(e**™ log (RT)), 


o—R Weer 


provided R > 2/T. Equation (3) can now be obtained by taking R = w/2 and 
letting w— ©. 
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On the Conjecture of Hardy & Littlewood 
concerning the Number of Primes of the 
Form n’? + a 


By Daniel Shanks 


1. Introduction. In a famous paper, [1], Hardy and Littlewood developed a 
number of conjectures concerning the twin primes, the Goldbach problem, and 
other unsettled questions. One of these, Conjecture F, concerned the number of 
primes of the form Am? + Bm + C. We reword this conjecture, and at the same 
time reduce its generality somewhat, as follows: 


Consecrure. If a is an integer which is not a negative square, a # —k’, and if 
P.(N) is the number of primes of the form n’? + a for 1 <n < N, then 


(1) PAN) ~ 3 he [ ines 


where the constant h, is the infinite product 


0 n= he-@)et) 


taken over all odd primes, w, which do not divide a, and for which (—a/w) is the 
Legendre symbol. 

In the trivial cases, a = —k’, since (k’/w) = +1 for every w, we have h, = 0 
on the one hand, and on the other there can be at most one prime of the form 
n — k = (n — k)(n +k). For any other a, h, > 0, and the conjecture indicates 
that there are infinitely many primes. But for no a has this been proven. 

In particular, for a = 1, since (—1/w) equals +1 or —1 according as w = 
4m + 1 or 4m — 1, we have 





(3) m= (1+ $)(1 — 4)(1 4+ 4)(1 + Ao) (1 — wy) --- = 1.37281346 --- 
and therefore (1) implies that 


(4) P,(N) ~ 0.68640673 I © dn 


log n° 
A. E. Western [2] verified that the number of primes of the form n* + 1 agreed well 
with the right side of (4) up to N = 15,000. 

In a recent paper [3] a sieve method was developed for factoring numbers of 
the form n’ + 1, and more generally of the form n® + a, and it was 
shown that the good agreement in (4) continues to hold out to N = 180,000; 
(N? + 1 = 32,400,000,001). This verification, however, was not applied to (4) 
directly but to the related formula, (7), given below. 

Let #.(N) be the number of odd primes, g, which are =< N, which do not divide 
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a, and for which (—a/q) —1. These are the primes which never divide n’* + a. 


It is well known that 


N 
(5) z.(N) ~! I dn 
2/2 logn 
and therefore (1) can be rewritten as 


P,(N) 
(6) FAN) ha. 


Likewise (4) can be rewritten as 


P,(N) 
#,(N) 


Since, in [3], we had P;(180,000) = 11223, #,(180,000) = 8178, and 11223/8178 = 
1.37234, the agreement with the right side of (7) was even better than could be 
expected. 

It is clear that the #,(N) in (6) could be replaced by the asymptotically equal 
4x(N) or by #.(N), (for the latter number we count the p’s such that 
(--a/p) = +1). But (6) as it stands is to be preferred for two reasons. First, #,(N) 


(7) ~ 1.37281346 ---. 





N 
is generally much closer to 4 / dn/log n than are either of the other two counts. 
2 


See [4, sec. 10 and Table 7] for a discussion of the case a = 1. Second, the ratio 
in (6) has a simple geometric interpretation in the algebraic number field R(+/ —a). 
See [3, p. 82] for a discussion of the case a = 1, the Gauss plane. 

In the present paper [5] we first develop an interesting and rapidly converging 
formula for computing the h, and we tabulate these constants for a = —4(1)4. 
We then present short tables of P.(N) and #,(N) fora = +2, +3, +4, and for 
N = 10,000(10,000)180,000 which show that (6) also gives good agreement in 
these five cases. Finally we present an elementary (sieve) argument which makes 
it plausible that the Hardy-Littlewood conjecture is true for every a. Further, an 
analysis of this computation enables us to isolate the essential difficulty in obtaining 
a proof. 


2. The Right Side of (6). To compute the ha. we will want the following 
Lema. For | x | <3, 


(8) 5 - (42) 


s=1 1 — Fog 





where the exponents b(s) are given by b(1) = b(2) = b(3) = 1, 6(4) = 2, b(5) = 3, 
b(6) = 5, and, in general, if d is an odd divisor of s and u(d) is its Mébius func- 
tion, then 


(9) b(s) = >? u(d)2", 


Examples of (9): A.) If s = p, an odd prime, d = 1 ord = p and [6] 
(9a) b(p) = (2? — 2)/2p = (2? — 1)/p. 
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B.) If s = 2", then d can only equal 1 and 
(9b) b(s) = 2°"/s. 


Therefore b(7) = 9 and b(8) = 16. 
Proor oF THE Lemma. After taking the logarithm of both sides of (8), 


(10) —In (1 — 2z) = > b(s) In [(1 + 2°)/(1 — 2*)], 


we expand both sides in Maclaurin series and identify the corresponding coefficients. 
This yields the condition, for s = 2“m, with m odd, 


1 moys ‘). 
(11) 2 > 7° ( 7 
Now applying the Mébius inversion formula we obtain (9). Since from (11) we 


also have b(s) < 2°/2s it follows that (10) converges if | x | < } and the steps may 
be reversed to yield (8). 


Now for any a ~ —K let p; be the odd primes such that (—a/p) = +1, let 
qi be the odd primes such that (—a/q) = —1, and let r; = 2,172,173, ---, Te be 
the (finite number of) primes which divide 2a. Further, for s = 1, 2, 3, --- , let 


mo =[ Qe 


the product being taken over the p’s and q’s in numerical order. Finally 
for s = 2,3, 4, --- , let 
(13) go(s) = (8) IT (1 — ri") 


where ¢(s) is the Riemann zeta function. 
THEOREM. If 


(14) fa” = £.(2)/La(1) and Ke” (s) = £(28)/La(8)ta(8) 
for s = 2, 3,4, --- , then 


(15) i = f™ F II [K.(s)~, 


where b(s) is given by (9). More generally, for more rapid convergence, we may select 
a positive integer u and define 


1 (w) - 0) “ (1 = 2 )-  @) . ( Le >) (z+ ). 
(16) f f I pi(pi — 1) f I] Pi pi — 1 


and 

(17) K,'(s) ae K.(s) II (1 + Pok.. ) ae K,(s) II (2 tt) j 
i=] pi hai 1 t=] pi a 1 

Then for every u = 0, 1, 2, --- 


, 


18) he = fe TT (K.'(s)). 
s=2 
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Proor. For every s = 2, 3,4, ---, 


to =[B0-3)0-HO-H)] 


and we easily verify that 


(19) 1 = K,(s) [J (? = ‘). 
Se gla 1 
We likewise find that 


> 9G) 


p—il 
so for any positive integer m, we have from (19) and (20) 
m C) _m CJ 3 b(s) 
he = ” Il [Ka (s) 4 II (1 i 2) (2: + ') ; Il II (2: + ‘) ; 
3=2 isu+ Pi] \pi — 1 end ime ti \p,? — 1 
Since m is finite the order of the products may be changed to give 
m c) me 8 b(s) 
ha = a Il [K.(s)? ’ II (1 - 2) ; IT (2: ‘) ; 
s=2 t=u+l Di amit \D* — 1 
Now every p > 2, and we may therefore use (8) with x = 1/p; to obtain 
m ° C) a | b(s) 
he — as II [Ke (s) : II II (2 ') , 
s=2 i=u+l1 s=m+1 ‘" + 1 


But it may be readily seen that the double infinite product on the right converges 
(monotonically increasing) to 1 as m — o, and it thus follows that the right side 
of (18) converges (monotonically decreasing) to h, as m— ~. 

The computation of the h, from (18) requires knowledge of the L,(s). Now 
every L,(s) has a Dirichlet series 


L,(s) = > d.(a)n 














with real periodic coefficients. Specifically we have 
Lis) =1-—3°+5*°-—7°+-4-, 
Lele Ne BHR AT? ht, 
L(s) =1—3°-—5°+7" +--+, 
L,(s) =1—5°+7° — 11° +-4-, 
L_»(s) =1—5° —7* + 11° +——4, 
Is) = 1 — 3° 45° — T° +-4-. 


The L,(1), which enter into f,°° as defined by eq. (14), may be obtained in 
closed form by use of Gauss sums and Fourier series, [7]. Specifically, for a > 0 
we have the simple 


(22) L.(1) = 


(21) 


aoeutl 
2-/a . 


where the q, for i < a < 100 are listed in Table 1. 
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TABLE 1 
a da | a qa | a Ga a Va 
1 gley ogg 6 | 51 6 76 6 
2 f) +4! bee 3 | 52 4 77 s 
3 1 28 2 53 6 78 4 
4 1 29 6 | 54 6 79 5 
5 2 30 4 55 4 80 8 
6 2 31 ef 8 81 6 
7 1 32 4 | 857 4 | 82 4 
8 2 33 s: Se 2 83 9 
9 2 34 Oe ced 9 84 8 
10 2 35 = 4 85 4 
11 3 36 4 61 6 86 10 
12 2 37 2 | 62 s 7 6 
13 2 38 6 63 4 88 4 
14 4 39 4 64 4 89 12 
15 2 40 4 65 s 90 8 
16 2 41 8 66 8 91 6 
17 4 42 4 67 3 92 6 
18 2 43 3 68 s 93 4 
19 3 44 6 69 s 4 s 
20 4 45 4 70 4 95 s 
21 4 46 4 71 7 96 8 
22 2 47 5 72 4 97 4 
23 3 48 4 73 4 98 S 
24 4 49 4 74 10 99 6 
25 2 50 6 75 6 100 4 
TABLE 2 
a he 
andl 0 
=3 1.38342429 
—2 1.85005441 
~~ 0 
0 0 
1 1.37281346 
2 0.71306310 
3 1. 12073275 
4 


1.37281346 





The L,(1) for negative a are a little more complicated and will not be listed 
here. As regards L,(s) for other values of s, Z,(s) is a well known function, but 
except for a few scattered results, [8], values of the other L’s do not seem to have 
been published. J. W. Wrench, Jr. has computed unpublished tables of L,(s) for 
a = +2 and +3. With his permission the author used these tables, together with 
(18), to compute the four corresponding values of h, in Table 2. The remaining 
entries, hg = h.. = ho = 0 and hy = hy, are trivial. 

The variation of the h, in Table 2 is notable. For example, there should be 
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more than two and one-half times as many primes of the form n’ — 2 as of the 
form n’ + 2. As a side remark, we note from (15) that fa” = 2¢.(2)+/a/mq, is the 
leading factor of h, . Thus for a > 0, n’ + a will therefore have few or many primes 
according as q. is large or small (relative to 2+/a/r). From Table 1 we see that 
there will be few primes for a = 2, 5, 11, 14, 26, 41, 89, and 194, (qi, = 20) and 
there will be many primes for a = 7, 37, 58, and 163, (qe; = 3). The famous func- 
tion of Euler, n’ + n + 41, equals 3{(2n + 1)* + 163] and its well-known richness 
in primes is thus closely related to the small value of qs; . This, in turn, is related 


in class number theory to the unique factorization of the integers in the algebraic 
number field R(«/— 163). 


3. The Left Side of (6). Tables of P.(N) and #,(N) fora = +2, +3, +4, and 
N = 100k (k = 1, 2, --- , 1800) were computed with an IBM 704 program based 
on the sieve method and the p-adic square roots of —a, [3, sec. 9]. At the same 
time the prime divisors of n’ + a which do not exceed N were counted, and from 
these counts the values of #,(N) are easily obtained. Summaries of these results 
are given in Tables 3, 4, and 5. In the last of these, the results for a = 4 are com- 
pared with the previous results [3] for a = 1. 


4. Both Sides of (6). In Figure 1 we plot P.(N)/#.(N) versus N together with 
the conjectured limits, h, , fora = +2 and +3. The casesa = 1 anda = 4, (which 
should be asymptotically equal since hi; = hx), are not included in this figure for 
clarity. If included, these two graphs would intertwine that for the case a = —3, 

5. An Elementary Interpretation. The over-all impression of the foregoing results 
is that (6) and its equivalent (1) are almost surely true for a = 1, +2, +3, 4. 











TABLE 3 
N | PN) | #2(N) P2(N)/#2(N) | P_2(N) #-2(N) | P_2(N)/#-2(N) 

10000 446 | 622 | 0.6737 | 1153 625 1.8448 
20000 817 | 1134 | 0.7205 | 2140 1140 | 1.8772 
30000 1180 | 1632 0.7230 3087 1631 | 1.8927 
40000 1494 | 2117 | 0.7057 3977 2112 | 1.8830 
50000 1821 | 2580 | 0.7058 | 4824 2587 1.8647 
60000 2160 | 3051 | 0.7080 5643 3041 1.8556 
70000 2489 | 3478 | 0.7156 6464 3481 1.8569 
80000 2823 | 3942 | 0.7161 7296 3927 | 1.8579 
90000 3139 | 4378 | 0.7170 8083 4374 1.8480 
100000 3422 4798 | 0.7132 | 8888 4808 1.8486 
110000 3721 | 5229 | 0.7116 9681 5242 1.8468 
120000 4027 | 5649 | 0.7129 10500 5682 1.8479 
130000 4347 | 6090 | 0.7138 11304 6117 | 1.8480 
140000 4652 | 6516 | 0.7139 12086 6533 | 1.8500 
150000 4966 6945 0.7150 | 12828 6956 | 1.8442 
160000 5250 | 7347 | 0.7146 | 13628 7362 | 1.8511 
170000 5522 | 7767 | 0.7110 | 14397 7763 | 1.8546 
180000 5847 | 8192 | 0.7138 | 15134 8184 1 


.8492 
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- TABLE 4 . ae ep 
the N pan) | eum) | Paanyeum | Pao | ay | Parveen) 
nes | | 
nat 10000 711 | 616 | 1.1542 | 850 | 620 | 1.3710 
aa 20000 | 1302 | 1136 1.1461 | 1569 | 1139 1.3775 
a, 30000 1851 | 1633 1.1335 2238 1637 1.3671 

4 40000 2378 | 2112 1.1259 2903 2108 1.3771 
—- 50000 2920 | 2575 1.1340 3550 2577 1.3776 
ted 60000 | 3428 3041 1.1273 4168 3030 1.3756 
-aic 70000 3967 3490 1.1367 4796 3466 1.3837 

80000 4463 3937 1.1336 5442 3935 1.3830 
90000 4941 4373 1.1299 6049 4374 1.3829 

d 100060 5426 4806 1.1290 6664 4819 1.3829 

- 110000 5917 5233 1.1307 7253 5247 1.3823 
sed 120000 6410 5665 1.1315 7874 5673 1.3880 
me 130000 6873 | 6105 1.1258 8491 6097 1.3927 
om 140000 7337 | 6532 | 1.1232 9073 | 6524 1.3907 
ats 150000 7823 | 6940 | 1.1272 9663 6950 1.3904 

160000 8302 | 7361 | 1.1278 10236 7363 1.3902 

_ 170000 8781 | 7768 | 1.1304 10799 7765 1.3907 

180000 9240 | 8195 | 1.1275 11354 8200 1.3846 
vith 

nich TABLE 5 
for | N | Pa) FAN) = n0W)| PAN)/s(N) | PAN) | Pr)/eN) | Px)/PA) 
—3, : | 

10000 | 870 619 | 1.4055 841 | 1.3586 | 0.967 

lee 20000 | 1554 1136 | 1.3680 | 1559 | 1.3724 1.003 
3 4 30000 | 2216 | 1633 | 1.3570 | 2268 1.3889 1.023 
wil 40000 | 2838 | 2117 | 1.3406 | 2952 | 1.3944 1.040 

50000 | 3459 2583 | 1.3391 3613 1.3988 1.045 

ae 60000 4083 3038 | 1.3440 4252 1.3996 1.041 
ww) 70000 | 4690 3485 1.3458 4888 1.4026 1.042 
ee 80000 5281 3933 | 1.3427 | «5513 | 1.4017 | 1.044 

90000 5903 4364 | 1.3527 | 6084 | 1.3941 1.031 

. 100000 6517 4808 1.3554 | 6656 | 1.3844 | 1.021 

110000 7099 5247 | 1.3530 | 7239 | 1.3796 | 1.020 

f 120000 | 7700 | 5675 1.3568 | 7795 | 1.3736 | 1.012 
D - 130000 | 8300 | 6103 | 1.3600 | 8369 | 1.3713 | 1.008 
7 140000 8893 6531 | 1.3617 | 8944 1.3695 | 1.006 

4 150000 | 9442 | 6941 | 1.3603 | 9505 | 1.3694 | 1.007 

9 160000 | 10008 | 7361 | 1.3596 | 10072 | 1.3683 1.006 

9 170000 | 10565 | 7770 | 1.3597 | 10658 | 1.3717 | 1.009 

; / 180000 11143 8178 | 1.3626 11223 1.3723 1.007 

8 ——s 

: We now offer a theoretical argument in favour of these asymptotic equations for 

0 all a. We will specifically carry it through for a = 1, but the argument is easily 

2 generalized. The case a = 1 is the only one which Hardy and Littlewood treated 
1 } in detail. Their computation, however, was deep and function-theoretic. In con- 

o trast, the present argument is elementary, [9]. It will be assumed that the reader is 


acquainted with the n* + 1 sieve which is described in detail in [3]. 


| 
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P(N) 


FaiN) Me 


h_,° 1.850054 


h_,* 1.383424 


h, = 1.120733 


h, = 0.713063 





- 50,000 100,000 150,000 200,000 - 


N-— 


Fig. 1.—The Hardy Littlewood Conjecture. 


Consider the infinite product (3) for h; , not in the form in which it was given 
by Hardy and Littlewood, (2), 


m= (14 52a)(1~52a)(+ ea) + aa) ea) 


since this masks its true nature; but in the equivalent form 


poe Gea ae. line) 
(3) (3) (3) G-a) O-a) 


or, even better, as 


1 2 2 
(i-3) 1 (-3) (1-%) 
(23) hy = ae ae | a te ee ee 137... 
pret (1-5 1-1 1-3) +. 1-7) 
2 3 5 7 li 13 
Now for a suitably large N let w* be the greatest prime satisfying w < N and let 


p* be the greatest prime of the form 4m + 1 which satisfies p < N. We write the 
corresponding partial product of (23), which approximates h, , as follows: 


1 2 2 2 
eee! pe y=) 0-5)0- a) 0-3) 











=O-DE-DO=d C8) 


Now this approximation to h; is in turn seen to be approximated (and we will 
inquire later as to the degree of the approximation) by N times the ratio of the 
primes which remain in two sieves, the Eratosthenes sieve (for all primes) from 
n = 1 ton = N’ in the denominator and the n> + 1 sieve from n’ + 1 = 2 to 
n’ +1 = N’ + 1 in the numerator. 
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Without attempting precision at this point—that is, without bounding the 
error—we note that in the Eratosthenes sieve one first strikes out the multiples 
of 2. This leaves N*(1 — 4) numbers (with an error of 0 or }). One then strikes 
out the remaining multiples of 3 leaving N’*(1 — 4)(1 — 4) numbers (again except 
for a possible end-effect correction.) Continuing with the primes 5, 7,--- , w* 
creates the denominator of (24). The latter therefore equals 


x(N*) — x(N) + E(N), 


the number of primes up to N* minus the number of primes up to N, with an end- 
effects error, E(N), which is not yet bounded. We note that 


oft) ali x x(N) ~ Na,(N) 


by the prime number theorem. 

In the n® + 1 sieve we first factor a 2 from all numbers where n = 2m + 1 
leaving N(1 — 4) of the numbers (except for an end-effect error). We then factor 
a 5 where n = 5m + 2 and where n = 5m + 3. This leaves N(1 — $)(1 — #) 
numbers (except for the end-effect error). Continuing with all primes of the form 
4m + 1; 13,17, --- , p* generates the numerator. The latter therefore equals 


P(N) — P(VN — 1) + e(N), 


the number of primes of the form n’ + 1 up to N* + 1 minus the number of such 
primes up to N with an end effect e(N ). 
Therefore, we may write 








_ 1, P(N) — P(N — 1) + e(N) 
(25) hy = lim AN) + E(N)/N 
while what we would like to write is 
_ ,. P(N) 
ae a 


Now by Merten’s Theorem the denominator of (24) is asymptotic to N *e*/log N 
where 7 is Euler’s constant [10]. Therefore the end effect, E(N)/N, is not negligible 


compared with #,(N). Instead we have 
E(N)/N jin arm 
(26) FAV) | 0.1229 = 2e 1. 


If we could show 


ood —— ~ 2" —1 
P(N) — P(N — 1) 
all would be well, but the difficulty of the problem is such that we cannot even 
prove that the left side of (27) is bounded from above. If we could do that, we 


would at least have P(N) — © but even this “weak” result eludes us. 
It is of interest to analyze this difficulty. Let 


(28) D(N) = P(N) — P(VWN — 1) 





(27) 
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and 


(29) sin) = w(1-3)(1-2)--- (1-3). 


Then the conjectured relation (27) is equivalent to the conjecture 


S(N) _ 
D(N) 
Now from the sieve for n’ + 1, [3], we can obtain an exact formula for D(N) 


by using the “integer part of x” function, [2]. Consider the set of numbers obtained 
from 


(30) 2e” = 1.1229. 


da FB Ie... p* 


by assigning (in all possible ways) 0 and 1 to the exponents a, b, c, --- . For each 
such d, let A; be the solutions of 


A’? = —1 (mod d) 
which satisfy 
0S A<d. 


Then if d is a product of a primes, we have 
: e N +A; 
(31) D(N) = > (—1) > ee 


It may be seen that if there are M primes of the form 4m + 1 which are <N, 
then there will be 2-3” terms in this sum. Even for a very modest N, say 15, we 
have p* = 13, M = 2, and there are already 18 terms. Specifically, 


» — an _[N+1]_[N+3]_[N+2 N+7 N+3 
pow) = im — [V4] [N48] _ [N42], (N47) [NS 
N+8]_[N+5 N +21 N+5 N +57 
-| +8] \"4 |+[ 26 |+| +| 65 | 


N + 47 N+18 N+8 N +83 N +73 
+| 65 |+| 65 | + | 65 | x Ee | 7 [| 


_|N+57)_ | N+47 
130 130 ‘ 
In general, it is easily seen, the formula for S(N) may be obtained from that 


for D(N) by deleting the A; and the square brackets. Thus for N = 15 in the ex- 
ample, we have 
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TABLE 6 
N SW) DW) S(W)/DW) 
100 16.261 15 1.016 
200 28 . 252 28 1.009 
300 39.800 42 0.948 
400 50.696 51 0.994 
500 61.344 62 0.989 
600 71.763 68 1.055 
700 81.656 78 1.047 
800 91.345 87 1.050 
900 101.075 92 | 1.099 
1000 110.901 102 1.087 
1100 119.913 112 1.071 
1200 129.451 122 | 1.061 
1300 138.223 128 1.080 
1400 147.754 140 1.055 
1500 156.790 150 1.045 





For N small, S(N) and D(N) are nearly equal; e.g., S(15) = 3.81, D(15) = 
As N increases, S(N) gradually pulls ahead of D(N), as is seen in Table 6. 
The end effect 


e(N) = S(N) — D(N) 
is given by 


(32) e(N) = 5 (-" E48 - [=++}. 
a “\d d ) 

Since the quantity in each brace is smaller in magnitude than unity, it is easy 
enough to bound e(N). What is difficult to obtain is a sufficiently good bound—that 
is, to prove in general, the extensive cancellation of terms of opposite sign which 
occurs in the sum of (32). The essential difficulty stems frozn the very rapid increase 
in the number of terms, 2-3”. 

Techniques of deleting or combining terms, in sieve formulations of related 
problems, have been devised by Brun and others [11] but to date nothing sufficiently 


sharp has been developed. A general assessment of sieve techniques given by Selberg 
[12] is not encouraging. 
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On the Propagation of Errors in the Inversion 
of Certain Tridiagonal Matrices 


By Arnold N. Lowan 


Abstract. When the differential equation of heat conduction is replaced by the 
implicit difference analog, one is led to the solution of Ay = b where A is a tri- 
diagonal matrix whose elements on the principal diagonal are =2 + 2r and whose 
elements off the principal diagonal are = —r. 


The system of equations may be solved by the following algorithm: 
B=u="Br, =m; n=—-B > w=(kht+raa)h, a= bu; 
Yr = Ze — VeYe » Yu = 2u- 


An upper bound of the round-off errors in the computed values of the y,’s is ob- 
tained. An actual test case showed that the theoretical upper bound is about four 
times larger than the true round-off error. Moreover, the theoretical upper bound 
does not seem to vary appreciably with r. 


When the differential equation of heat conduction 


aT eT 
— = ¢— szs 
x os? eSs86, t>@ 


is replaced by the “implicit’’ difference analog 


Fant = lan o 








i = (Az) [Tn-1,0+1 - y > papery + T m+i.n+i + ae _ SF i. ts a Tau 
a 
m= 1,2,3,---M, Ae= 57S 


or 
(2 + 2r)T mins — 1(T mans + Topix) = (2 — 2r)T an 


+ nT eas + : — 


where T'm.n = T(mAz, nAt) andr = cAt/(Az)’ it is a known fact that the difference 
scheme (1) is unconditionally stable [1]. If the desired solution is required to vanish 
on the boundaries x = 0 and z = a, the system of equations (1) may be written 
in the compact formt 


(1*) ATax: = BT, =b (say) 


where A is a tridiagonal matrix whose elements on the principal diagonal are 
= 2 + 2r while the elements off the principal diagonal are = —r. 





Received July 31, 1959; revised form February 9, 1960 
+ When the temperature is prescribed on the boundaries, equation (1*) is essentially 
unchanged except for the fact that the first and last components of b are slightly altered. 
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The system of equations (1) may be easily solved by the following algorithm [2] 
2 


(2) ep Sh se eee k = 1,2,3,---,M; B, = 
Bea 
(3) n= —— k = 1,2,3,---,M; 
Bi. 
1 bi 
(4) Ze = — (be + ree) k=1,2,3,---,M; a=- 
Bx u 
(5) Yu = Ze — Ve Yer k = 1, 2,3,---,M; Yu = 2u 


where we have written u for 2 + 2r and we have denoted the components of T,., 
by yz . The question arises: if the computations involved in the above algorithm 
are carried to p decimals (i.e., if products and ratios are rounded to p decimals) 
what is the upper bound of the round-off errors in the computed values of y;? 

In the derivation of the desired upper bound we shall require a lower bound 
of the 8;’s and upper bounds of +; , z, and y, . If in (2) we put k = 2, 3, --- we get 


2 2 


r 
ee ee 

2 

(6) B=u-e 


From the first of the above equations it is clear that 62 < 8,. From the first two 


equations it follows that 
p- a= (2-2) <0. 
1 Be 


Thus 6; < 62. Similarly it may be shown that 6, < 83, --- , Be < Bx-+. Thus the 
B,’s form a monotonically decreasing sequence. It may be readily shown that the 
lower limit of the sequence, to be denoted by 6 , is the larger of the two roots of 
the quadratic equation 


(7) e —(24+2r)¢+r =0. 
Accordingly 
(8) Be =1l+r+ V1 + 2r. 


From (3) it follows that | y. | < r/8; . If then y* denotes an upper bound of | +: | 
we may put 





(9) == as 
Be l+r+~V7i+2r 
If in (4) we put k = 2, 3, --- and subsequently eliminate z., z;,---, z:. we 
ultimately get 
b; rbyp-4 by» r* by 





4=— 


a Oa © ow wees » weer 


wl 


wl 
bo 


Fi 


(1 


Zk 


er 


wl 


= 
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whence 


? ay b* 
a= 1 o@e r) |=e. = 
|2s| Bs “| +5, zat (F Bs Sande i —¢ 
Bs 


where b* is the largest of the absolute values of b, . If then 2* denotes an upper 
bound of | z, | we may put 








b* 
(10) z* = aon 
Finally from (5) we readily get 
Yu => 2m 
Yuu >= 2Zu-1 — Yu-%u 
Yu—2 = 2u—2 — Yuu + Yu-2¥ uz 
Yr = 21 — Vide + yevets — +°* (1) eve -- + Y-te- 


From the above system of equations it is clear that 

y* = 5°(1 + 7* + 7% --- + y™”) 
(11) ee er h aa S: ee 
Bs —T a" : (Bs — r)? 


Bs 





is an upper bound of the absolute values of the y;’s. 

We now turn to the evaluation of upper bounds of the errors in the §;’s, 7:’s, 
z’s and y,’s. It will be convenient to denote by E(6;,) the absolute value of the 
error in 8, and by E*(8) an upper bound of the errors in the ;’s. A similar notation 
will be used for the +;’s, z,’s and y;’s. From a we have 


— E(f) +6 





: a 
E(B) = gp EA) +és = 
where 6 = 3 X 10°” is the maximum round-off error. Similarly 


2 
=. o(B2) + 8 


-5 “la E(6) +a] +3 


1 5 =) E( 
(+m) +) 


Proceeding in this manner we ultimately get 





2 
E(63) = 5p E(b) +és 








2 2\ M-—2 2\a-1 
E(B.) < E +(4 :) + ss (*) Je+(S) E(6,) 
* ~ * 
hol Bx” 
oa beeen a8 
By? 
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where we have neglected the second term of the above inequality since r < By . 
Thus 
8 2 
(12) E*(g) = —~— 8 
Bs 


2— 72 


is an upper bound of the absolute values of the E(6,)’s. 
Consider now the evaluation of E*(7). From (3) it follows that 


E(y) = re E(B.) +8 


< 33 B*(8) +6 
whence 
— 7 ae rT ’ Bx” 
E*(y) = 3,3 B*(8) +i= Ba? a2 — 7° +4 
(13) 


=(t+53t)* 


Consider next the evaluation of E*(z). From (4) we get: 





E(z) = x { (bi + rzr+) E(B) + BLE (bi) + rE (z4))} + 6 
1 b* *) P* 1 ns Zz. 
$ 7a (0 + rz*)E is, th, were 
- oa (b* + rat) Pe be + hE) + 2 Blas) +3 
b* -2 rZ e 
=(1+5% reel BM(b) + 2 Bla). 


Proceeding as in the evaluation of E*(8) we ultimately get 








(14) E*(z) = fi (1 + b* = s etme _ ee 
By — 7 -_ 
If in the last equation we replace 2* by its expression from (10) we ultimately get 
Bx | b*Bx | E*(b) 
(15 E*(z) = ——_] 1+ 1+ 
(15) (2) By — 7 (8, — r)(By? — 7°) By — °° 


If on the other hand we replace z* in (14) by Z, the largest absolute value of the 
z's we obtain 


' Be b =tit E*(b) 
* E*(2) = —"* E. 
) (2) (ith . A et 


While the expression in (15) is an upper bound of the errors in the 2,’s, it is rea- 
sonable to refer to the expression in (15*) as the least upper bound of the errors 
in the z;’s. 

Finally, consider the evaluation of E*(y). From (5) we get 


E(yx) = E(ee) + veE (yer) + ye (ve) + 6 











get 


the 


rea- 
rors 
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whence 
(16) E(y.) S E*(z) + y*E(yeus) + y*E*(y) + 6. 


Substituting for E*(z), y*, y* and E*(y) their expressions from (15), (9), (11), 
and (13) the last inequality becomes 


Be b*B E*(b) 
En) & 5 lt =n nit tks a 








(17) 
, 
Proceeding again as in the Dries, of E*(8), the Rib aha ultimately yields: 


*(,) — Bx {Be | b*Bs ] 
BX) = le 7 + & one, 








(18) ae 1 " 
om 
+ Gea + gtap) tft gape. 
If, on the other hand, we substitute for E*(z) in (16) its expression from (15*) 
and replace y* by Y the largest of the absolute values of the y,’s, while y* and 
E*(7) are replaced by their expressions from (9) and (13), we obtain as the counter- 
part of (17) 


Be b* + a4 E*(b) 
=< 
E(y) S 5 +h ial Ee 


+¥(1+57" 











(17*) 





r 
6+i+—E . 
—— 3) +842 Blum) 
Proceeding again as in the evaluation of Z*(8), the last inequality ultimately yields: 
Bs { Bx | b* + a4 (1 arn) hs 
E*(y) = 1+2—5/+Y(1+ +1 
= Rare TL te Bs 
= 
— E*(b 
@ —7 (b). 
It will be convenient to rewrite the last equation in the form 
(19) E*(y) = So(r) + b*S,(r) + YS.(r) + ZS;(r) + E*(b)S,(r) 


where 





(18*) 
+ 

















S49) = mo rt act r? 
We? aoe he =F) 

woaad Br) = (1+ 555) 
we? ie es a 
Si(r) = 


(6, — vr)?” 
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In (19) in conjunction with (20) we have an upper bound of the round-off errors 
in the values of the y,’s—the solutions of Ay = b. In case the },’s are exact we 
must of course put E*(b) = 0. 

To test the formula (18*) the exact components of b were computed from 
Ay = b where A is a 20 X 20 tridiagonal matrix of the type above considered 
with r = 1 and r = 2 and the 20 components of y were arbitrarily assigned; the 
values of the components were then calculated by the above algorithm in terms 
of the exact values of b, . The computations were carried to eight decimals. The 
maximum discrepancy between the exact values of the y,’s and the corresponding 
computed values was two units in the last place. The upper bound of the round-off 
errors evaluated from (19) in conjunction with (20) was eight units in the last 
place in the case r = 1 and seven units in the last place in the case r = 2. The 
estimated upper bounds of the round-off errors must be considered as indeed very 
close to the actual round-off errors. 


The writer wishes to express his appreciation to Roy Steeves who carried out 
the above-mentioned test. 


Yeshiva University 
New York, New York; and 
AVCO Corporation 
Wilmington, Massachusetts 


1. The implicit scheme was first suggested by J. CRANK and P. Nicouson in the paper en- 
titled, “‘A practical method for the numerical evaluation of solutions of differential equa- 
tions of the heat conduction type,’’ Camb. Phil. Soc. Proc., v. 43, 1943, p. 50-67. Its stability 
is discussed by G. G. O’Brien, Morton A. Hyman, and Sipney Kap an in “A Study of the 
numerical solution of partial differential equations,’’ Jn. Math. and Phys., v. 29, 1950/51, 
p. 223-251. It is also discussed in the writer’s book The Operator Approach to Problems of Sta- 
bility and Convergence, Scripta Mathematica, 1957. 

2. See for instance M. Lorx1n’s ‘‘The numerical integration of heat conduction equations,” 
Jn. Math. and Phys., v. 37, 1958, p. 178, and R. D. RicutmyEr, Difference Methods for Initial 
Value Problems, Chap. VI., Interscience, N. Y., 1957. 
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Tabulation of Coefficients for Operations on 
Taylor Series 


By Daniel C. Fielder 


1. Introduction. When certain operations are performed on polynomials or power 
series, the results often assume the form of a series of ascendingly indexed terms. 
A typical term, the kth, appears as the product of two distinct types of factors. 
One of these factors is a multiplier peculiar to the particular operation. The other 
factor is a group of literal coefficients which forms a kth order symmetric function 
of weight k and degree h in those literal coefficients. A detailed examination reveals 
that the second factor is composed of sums of combinations of appropriately sub- 
scripted and exponentiated literal coefficients. This second factor, which will be 
called a literal group, has a form common to many operations on polynomials and 
power series. The multiplier of each combination within a group is observed to be 
either a positive integer or a positive common fraction. 

As is well known, the weight of a term in a symmetric function is the sum of 
the products of corresponding subscripts and exponents. The degree is the sum of 
the exponents of a product combination of literal coefficients: As an example, for 
k = 6 and h = 3 the kth term might be {6} (gigs + [2]g.geg; + [4]g2°). The weight 
of a product combination of literal coefficients, say the gg, term, is 1 X 2 + 
4 X 1 = 6, and the degree of a product combination, again using the g,°g. term, 
is 2 + 1 = 3. For uniformity, the identification of enclosing symbols with alge- 
braic terms used in the above example is used throughout the paper. 

Several direct and inverse operations are investigated, and methods for finding 
general terms for groups and multipliers are discussed. For convenience, literal 
groups are tabulated. Because of the common properties of the literal groups, a 
single table of literal groups suffices for all operations described herein. 


2. Power Series of Power Series. The first operation to be considered is 


(1) (diy + bey? + bay? +--+) = (max + mex” + mgt + ---), 
in which 
(2) y = qx t+ gx +gn +---. 


This is essentially 6.362 of Adams [1]. 
A direct calculation of the first few terms leads to 
mi = {bi} (gx), 
mz = {bs} (gr) + {bi} (ge), 
(3) ms = {bs} (gr) + {2b} (gig) + {bi} (gs), 
ms = {ba} (gr') + {3bs} (gr'g2) 
+ {2bs} (gigs + [3]g2') + {01} (gs), 
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Table 1 presents a general form for the literal and numerical (braced) coefficients 
up through index 12. For any {t,,} of Table 1 except {t,,}, the numerical coefficients 


for the operation described by (1) and (2) are hb, 


. For {ti}, the numerical co- 


efficients are b; . 


TABLE 1. Values of literal groups, k = 1(1)12 





m, 


= {tin} (m), 

= {te} (gi?) + {tn} (g2), 

= {tas} (g:*) + {ts2} (gig) + {ta} (9s), 

= {ta} (gi*) + {ts} (g'g2) + {tee} (gigs + [3]gs*) + {tu} (gs), 


= {tes} (gi>) + {tea} (g:°g2) + {ts3} (gigs + gigs") 
+ {tse} (gigs + gogs) + {ter} (Gs), 


(lio) + {tes} (gi‘ge) + Uh Ga) gs + [3]g:°g2") 
+ {tes} (gigs + [2loigegs + [4]g2*) 
+ {tee} (99s + g2g4 + [3lg:?) + {ter} (ge), 


{tr} (gi) + {tre} (i592) + {ts} (gitgs + [2]g:*92") 
+ {ta} (gi'gs + [3lorge9s + gigs*) 
+ {ts} (orgs + [2logegs + 92’9s + 919s") 
+ {te} (gige + gogs + gags) + {tn} (gz), 


{tes} (gi®) + {ter} (g:° go) + {tes} (g:'gs + [§]gx*g2") 
+ {tes} (gitgs + [4]o:°929s + [2]9:°92*) 
+ {tes} (gi'gs + ([3lgr°g2gs + [3]ox*gs? + [Blogs gs + [4]gs*) 
+ {tes} (9296 + [2]oigegs + [Zlowgags + go’gs + gags") 
+ {tse} (9ig7 + g2ge + 939s + [3]gs”) + {ter} (gs), 





tos} (g:5ga + 
+ {tos} (gigs + 
+ {toa} (gi°ge + 


(tal tor) + {ts} (gr'g2) + {tz} (gi.gs + [Bloi°g:") 


5lgi‘gags + [s"]g:°92*) 
4lg,*gogs + [2]gx°gs” + [6]g:°g2"9s + gigs") 
3]9:°929s + [3]g:°g2g4 + [Blgigegs? + [3]oig2’gs + 92*gs) 





+ {tos} (91°97 + [2]orgege + [2loigsgs + gigs? + 9279s + [2]orgsgs + [4]gs°) 
+ {ta} (gigs + gog7 + gage + gags) + {tor} (go), 


= {tor} (g! 1 ©) + {too} (gi8ge) + {tos} (gigs + [$]g1°g2”) 
+ {ho.7} (gi°9s + [6lgs®°g2gs + [5]gi*g2*) + {to,6} (gi°gs + [Slor‘gags 
+ [$lgi*gs? + [1O]gi*go*9s + [$]g:°ge*) + {tos} (gitge + [4]9:°9296 
+ [Algs*gag + [6]g:*g2?gs + [6]g:*gs’g2 + [4lqo*gagi + [4]92°) 
+ {toa} (gi'g: + [Blosg.ge + [8]orgs9s + [$]gx°92 + [Sloi92"9s 
+ [6]gige92gs + go°gs + [8]ge’gs? + gigs®) + {tos} (gr gs 
+ [2]gigeg7 + [2loigsg6 + [2loigags + [2loegsgs + goge + 92°Qe 
+ gs ?g4) + { tio, 2} (9199 + gogs + 9397 + gage + [3 ]9s”) + {tho, 1} (910), 


= {tuu}(g't) + {t.10} (gi%g2) + {t.9} (gi8gs + [4]gr7g2) 
+ {tus} (gr'gs + [T]gi°gegs + [Tlgs®g2®) + {tu,7}(gi'9s + [6]loi®gegs 
+ [S]gi'gs? + [15]g:*g’gs + [5]g:°g*) + {tu,6} (gi'gs + [5]gi‘gegs 
+ [5]gx*gsg1 + [1O]g:?go*gs + [LO]g:°go7g4 + [10]g:*g29s" + gigs?) 
+ {tus} (Qi‘gr + [Algi’gogs + [4] gi°¢s9s + [2]g:°g.? + [6]9:°92"95 
+ [12]9:?gegsgs + [4]gigo*gs + [6]gigsgs? + gotgs + [2]9:°9s°) 
+ {tus} (gi'9s + [B]o:*g297 + [3]g:°osg6 + Blox’ 949s + [8]g:927G6 
+ [6]oigegsgs + [3loigegs? + [3lgrgs’gs + go°gs + [Slgs*gsgs 
+ gs*g2) + (fu, 3} (gi'go + [2]lgigegs + [Zloigsgr + [2loigage 
+ gigs? + 9297 + [Zlgegag6 + [2lgogugs + gags? + 9s°Qs) 
+ {ti.2}(gigio + gogo + gags + gagz + gsgs) + {tua}(gu), 











3°) 





eee 
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TABLE 1—Continued 


My = {ter2}(g'?) + {teu} (g,°g2) + {te.0} (gi'gs + [$]os*92") 

+ {t0}(g:'gs + [Bloggs + [Plgi®gs*) + {ts} (orgs + (Flos'oe? 

7]g:°92gs + [*P lox*ge* + [21]gs°92*gs) + {h2,7} (gigs 
3]g:’92 + [6loi°g29s + [6lgigags + [15lgs‘gx*gs + (15]gi*geg:" 
20]9:°92"gs) + {hee} (gi'g7 + [elos‘ge + [YP los*gs® + [5lor*oags 
5)gi‘gags + ([5lgigs'gs + [10]g:*g2*gs + [10]9:792*94 
15]g:°92"gs* + [20]9:*g2gags + [t]g2") + {ts} (gitgs + go'gs 
2]92*93? + [4]gi*g2g7 + [Aloi*gags + [4]gs*gags + [4]gg2"os 
4lgigegs® + [6]g:°g2"96 + [6lor*gege + [6lg:*92"Gs 
12]9:*92929s + (12|gigs*gags) + {he} (gi'gs + g2"Ge 
$)g:7gs? + [8]:*92 + [Slorgags + [3lg:*ga97 + [3lgigs"g: 
3]g:°gage + [3lgi9s"9s + [Slgigags? + [3l92*oa0s 
3]gegs'gs + [6logegags + [6lowgegags + | h1gs*) 
{ties} (gi°gi0 + g2°gs + gogs® + gags + [2logag + [2logags 
[2lgigegr + [2]oigsge + [2lg2gag: + ([2losgags + [2losgegs 
i + gogo + gage + gags + gsr + [4]g0*) 
12,15 (G12). 








+4+++4+44+4+4+++4+4+4+4+ 





3. Reversion of Power Series. (See 6.26 of Adams [1] and Van Orstrand [2].) If 

(4) z2=2— 92 — gr — gr —--- 

is reversed to yield 

(5) w= 2t me t+ me + me t+ ---, 

the result is similar to that noted in the preceding section. The first few terms are 
m, = {1}(g:), 
me = {2}(gr) + {1} (g2), 

(6) ms = {5}(g:') + {5} (gig) + {1} (gs), 
ms = {14}(gr') + {21} (gi ge) 


+ {6} (gigs + [B]go") + {1} (ga), 


Table 1 applies as before. For any {t,,} of Table 1 except {t.}, the numerical co- 
efficients for the operation described by (4) and (5) are the familiar binomial 


: k+h ae ene 2k 2 2k 
coefficients “ 4 a) For {tix}, the numerical coefficients are | ( r) (, r .) 


4. Development of a Power Series Ratio. The four succeeding operations are 
based on a tabulation derived from the ratio of two power series. A short discussion 
of this tabulation follows. The two power series and the quotient P(z) are 


(1) Pe 7 Plt) = Cot Cre t Ca + +. 
0 1 2 see 











342 DANIEL C. FIELDER 


If the indicated division is performed, a triangular tabulation array can be obtained 
as shown below. 

CoBo — ae = 0, 

CBo + CB — am = 0, 
(8) CBo + Cipi + CB2 — a = 0, 

C'sxBo + Cp + C182 + CoB; — a, = 0, 


Thus, if the a’s and @’s of the array (8) are known, the C’s can be found by the 
division indicated in (7). 


5. Direct and Inverse Summation Formulas of Power of Roots. In a previous 
paper, the author [3] discusses means of finding sums of powers of roots of the 
integral rational function 


(9) f(x) = age” + az” + --- +a +4, 


vy elaborating on a method credited to Newton [1]. Newton’s method utilizes the 
tabulation, 
la; _ Sido = 0, 


2a2 ++ Sid; + Sedo = 0, 
3a3 + Sia, + Sea; + Spa = 0, 


(10) 


where the sum of the bth powers of the roots of f(x) are 
(11) Ss = ret re tre + eee 


A comparison of (8) and (10) suggests that it should be possible to find a 
P(z) whose coefficients are functions of the a’s of (10) and are simply related to 
S,’s. Such is indeed the case for the determination of the direct summation formulas. 
For convenience, (10) is rearranged and restated as 


m + _— = 0, 
2m: + mgi + 292 = 0, 
(12) 3m; + 2mogi + mMige + 3g; = 0, 


where g, = (a;/do) and m = (S,/k). 


For the direct summation, it is seen that Bp = i, C, = (kK + 1l)meu, Be = oe, 


and a = —kg,. Thus, P(z) becomes 


(13) P(z) = —“2— 2gi2 — 3gaz — 4giz’ --- 


l+nztmetge+t+--- - 
The first few terms of the indicated division are 


— gt (gr — 2g2)z + (— gr + 39:92 — 3gs)2 
+ (gr — 4gr'ge + 4gigs + 2ge’ — 4qs)e + ---, 





(14) 
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where —g, = Co, 9° — 292 = C,, etc. Thus, 
m = {—1/1}(g:), 
mz = {1/2}(gr) + {—1}(g), 
(15) ma = {—1/3}(gr) + {1} (gg2) + {—1} (gs), 
ms = {1/4} (g:*) + {-1} (gigs) 
+ {1} (gigs + [1/2lg2") + {—1} (gs), 


It is apparent that (13) not only yields the desired summations but serves as a 
device for generating the literal groups of Table 1. For any {t,,} of Table 1 except 
{tix}, the numerical coefficients for the direct summation are (—1)". For {t,,}, the 
numerical coefficients are (—1)*/k. 

It is not possible to use (12) directly to find a P(z) whose coefficients lead to 
the inverse summation forms. However, a two-step solution leading to the inverse 
summation forms can be found by rearranging (12) as 


g + B; = 0, 
(16) ge + my91 + 2B, = (), 
gs + mige + 2mg, + 3B; = 0, 


where 
1B, = lm, + 09, , 
(17) 2B, = 2m. + Igo, 


The relationships g, = (a;,/do) and m = (S;,/k) apply as before. 
For the inverse summation, it is seen that B® = 1, Ck = gi, Be = km, and 
oa. = (k + 1)Bys4,. Thus, the P(z) becomes 


—B, —_ 2B.z om 3B;2 = 4Byz = 20 
P(z) = ——— 


(18) Lt mz + 2m2 + 3m;22 4+ --- 
q+ gz t+ gst t+ met. 





The first few terms are 


( 19) —B, oe (Bym, —/_ 2Be)z + (—Bym; + 2Bom, + 2Bym2 “> 3B;)2 + oO? . 


Through use of (17), the B terms are eliminated and the result is 


f—1/1!}(m), 
{1/2} (mm) + {—1/03(m), 
(20) gs = {—1/3}(m)*) + {1/1} (mm) + {—1/08 (ms), 
gs = {1/49 (my") + {—1/28 (mm) 
+ {1/11} (mms + [1/2}ms) + {—1/03 (m,), 


2 
= 
II 
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Except for different numerical coefficients in braces and the interchange of the roles 
of the m’s and g’s, (20) is identical with (15). The discussion foliowing (15) thus 
applies for finding the literal groups of (20). Table 1 is immediately adaptable with 
appropriate changes of notation. For any {t,.} of Table 1 except {i}, the numerical 
coefficients for the inverse summation are (—1)"/[(k + h) — (k + 1)]!. For 
{tx}, the numerical coefficients are (—1)*/k!. 


6. Logarithms and Exponentials of Power Series. Expression 6.364 of Adams 
[1] is essentially 


(21) n(l + giz + gee” + gee? + ---) = mez + me + me + --- 


Through a comparison of Adams’ table accompanying 6.364 with (8) and (12) and 

an application of the reasoning of the last section, the arrangement of the m’s and 

g’s of (21) becomes that of Table 1. For any {t,} of Table 1 except {t,:}, the nu- 

merical coefficients are (—1)"**. For {t,.}, the numerical coefficients are (—1)***/k. 
The exponential form is given as 


(22) exp (giz + gee’ + gsz’ + +++) = 1 + mez + me + me + paw 


Again through use of the methods of the last section, the arrangement of the m’s 
and g’s becomes that of Table 1. For any {t,.} of Table 1 except {t,}, the numerical 
coefficients for the exponential form are 1/[(k + h) — (k + 1)]!. For {it}, the 
numerical coefficients are 1/k!. 


7. Calculation of Numerical Coefficients within a Literal Group. It becomes 
apparent that in any literal group except those which are coefficients of the {t,,}’s 
the term gi ‘g:-n+: always appears. This term has weight k, degree h, and multi- 
plier unity. The letter combinations in any literal group consist of all the combi- 
nations of g’s which simultaneously have weight k and degree h. The quotient 
formed by (h — 1)! divided by the products of the factorials of the exponents of 
a combination of g’s determines the bracketed numerical coefficient for that com- 
bination of g’s. The single term of the g;* group (i.e., the coefficient of {t,}) has 
a multiplier of unity. This method for finding the terms of a literal group is essen- 
tially that outlined in Fielder [3]. 


8. Conclusions. It can be concluded that there are algebraic operations on poly- 
nomials or power series which have common forms for the literal groups. The 
braced coefficients are peculiar to particular algebraic operations. For a given 
operation, all numerical coefficients except {t,,} follow one combinatorial rule, 
while the {t,,} coefficients follow a different combinatorial rule. Since available 
tables of many algebraic operations which follow patterns similar to those dis- 
cussed herein are terminated at index four or five, the information contained 
herein should be useful in extending the available tables. It is anticipated that 
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many other algebraic operations on polynomials can be found which aré members 
of the operation class of this paper. 
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Explicit Solutions of the One-Dimensional 
Heat Equation for a Composite Wall 


By Marcia Ascher 


1. Introduction. Explicit numerical solutions of the equation of heat conduction 
in a wall of one material have been widely discussed in the literature. Consideration 
of the forward difference equation studied in references [2], [3], [4], and [6] suggests 
a variety of ways to handle the solution for a composite wall. This paper is a study 
of the convergence, stability, comparative accuracy and comparative computing 
time of three explicit numerical solutions of the heat equation for a wall composed 
of two materials. 


2. System of Equations. The equation for the one-dimensional flow of heat is: 


(1) Cs Ps > = < («. =) where a,_; S x SQ; s=i1,2 


0O<tStp 
with the condition at the interface 


_fau\ _, fdau 
(2) ky (2) = ke (), 


where p,, c,, and k, are constant with respect to time and temperature but may 
be different for each material. 
We will assume the boundary conditions: 


u(do,t) = constant, ¢ 


IV 


(3) 
u(a,,t) = constant. t 2 0 
and initial conditions: 


(4) u(x, 0) = constant; a < 2 < a 


Let each material’s thickness, a, — a,-1, be divided into N, equal parts of Az, , 
and ty into equal parts of At, . Let i denote the subscript associated with the space 
variable and j the subscript associated with the time variable. Let the solution of 
(1)-(A4) be called T'(2, ¢). 
Taylor series expansions of 7; j4: , Ti4:,;, and T;.,; , about 7;; are used to 
obtain 
r, At 


(5) Tint = [Ting — 2764+ Tad t+ Ti + hh 





where r, = k,/p.c, and 


A’aT _r, Ata,” oT a 
E, = whe : sed < + terms of higher order. 








Omitting E, , equation (5) gives a difference equation for finding the approximate 
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solution of (1)—(4), 7,541, when x; , 2:4 , and 2,4; are in the same material. This 
is the same as the forward difference equation for a one-material wall. 
The following equation for 7’; ;,: at the interface is derived in a manner similar 


to that used by M. Lotkin [5] in his discussion of an implicit method for a wall 
of two materials. 








[ (ress — Ti) e + (Tiny - Tu) ) ge | aa 
+2 - 1 = * 4 
adiatinoand wees + Antipr— ienalnarens 
where 
Ds 2At mm [~* C; p, At a Py _ Anh, (= :) 
re Axe Co po a Ax Ci pr 4 Fo “e ox? 

24 Ox* 4 OF 





an 2k ke Axe’ : 
a ler + mT (ery | + terms of higher order. 


Omitting EZ, , equation (6) gives a difference equation for finding 7';,;,, when 2; 
is at the interface, x;_; is in the first material, and z;,; is in the second material. 


3. Definition of Methods. Stability is maintained in the explicit numerical solu- 
tion for a wall of one material by choosing 


a — pchx” 
(7) At 3 => 





Using equations (5) and (6) as our basic computing equations, three different 
means of choosing At will be defined and thereby different computational schemes. 
We will arbitrarily assume that r; > r2 and, for simplicity, will confine the dis- 
cussion to cases in which r,/r2 is an integer R, and +/r,/r2 is an integer. 

If Az is specified as the thickness of each lamina within the wall, equation (7) 


gives two different maximum usable time increments depending on the properties 
of each material; 





2 2 


Ax Ax 
, 2r; and 7 2re 





Meruop |. In the first method an attempt is made to circumvent the difficulty 
of having two At’s by letting Ax, = Az and redefining Az, such that Ar. = Az,/+/R. 
This increases the number of laminae in the second material but yields only one 
time increment, At; = At, = At. The computations would take place as follows: 

Given Ar = Ax, r; and rp 
1. Compute Az. = Ax,/+/R 
. Find At = At, = Ake 
. Set time equal to At 
Use equation (5) to find the temperatures in material 1. 
. Use equation (6) to find the interface temperature. 
6. Use equation (5) to find the temperatures in material 2 


oe Ww bo 
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7. Advance the time by At. 
8. Repeat steps 4 to 7 until the temperatures at ¢y have been computed. 

Mertuop 2. In the second method At is chosen by evaluating Af, and Af, and 
using whichever is less. The computations would then be: 

Given Ar = Az, = Am, mr and rz 
1. Choose At = min (At, , Ate) 
2. Proceed as in steps 3 to 8 of Method 1. 

Meruop 3. In the third method both time increments are utilized by using the 
smaller increment only for those points at which it is necessary and the larger 
increment for the rest. The smaller time increment must be used for all points in 
the first material, at the interface, and for enough points in the second material 
to enable a smooth transition. For example, when R = 4, those points on the grid 
in Figure 1 denoted by dots are computed using A/,; and then those denoted by 
crosses are computed using Af, . The computations would proceed as follows: 

Given Ar = Az, = Azz, 1; and rz 

1. Compute At, 

2. Compute At 

3. Set time = 0 

4.SettQ=R-1 

&. Advance time by At, 

6. Compute the temperatures in material 1 using equation (5) and At, . 

7. Compute the interface temperature using equation (6) and Af, . 

8. If Q = 0 proceed to step 12. 

9. Compute Q points in material 2 using equation (5) and At, . 

10.Q-—1-@. 

11. Repeat steps 5 to 10 until indicated by step 8. 

12. Compute temperatures in material 2 using equation (5) and Af. 

13. Repeat steps 4 to 12 until the temperatures at tp have been computed. 


4. Convergence of Solutions. 
THeEoreM. If there exists a solution of the system of equations (1) to (4) which has 
bounded derivatives &°T /at’, a°T /dx*, and &T/ax* in 0 S t S ty, a S x < a; and 





Fie. 1.—Grid points for R = 4. 
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a < x S a, then the solutions obtained with methods 1, 2, and 3 converge to the true 
solution. The rate of convergence is O ( Az’). 

Proof. Let B, = | upper bound on #T/at |, Bs = | upper bound on 37'/az" | , 
and B; = | upper bound on 4°7'/az* | . A barred 2 Babe denotes that it is evalu- 
ated somewhere within the intervalO StS tp,amSr<am,au<2S Mm. 
Define the error at the point 2; , t; to be e,; = Ti; — ui. Here u;; is the true 
solution of (1)—(4). The error arising from the use of equation (5) satisfies the 
following equation: 


2 At mT . OT 
(8) Gin = = leiss.3 — 2:3 + ec-asl + 65 + af _ 7 Athz, T 





























12 dat 
r, At 2r, At At Af r, Atha,” 
(8’) ein S Cia + i -_ i + ia 615 +> x Bs ae va, B,. 
Let a; = max; | e¢;;|, then 
r, At| 2r, At .AtAz,” 
(8”) coos 5 || eA +|1 — H+ (eA a +4 5 ne 12 B.. 





In methods 1 and 2, At < par (s = 1, 2) and therefore r,At/Az,’ < 4. This 
causes each of the terms within the absolute value signs in equation (8”) to be 
positive and so they may be eliminated giving 


At 4 
(9) €iin, S aj + — By + -— Bo. 


In method 3, the larger time increment, At, , is only used at points for which 
At, < Az,’ /2re is satisfied and whenever At, is used, in the first material or for the 
transition values, At, < Az,’/2r,(s = 1, 2). Therefore, whenever equation (5) 
is used At < Az’/2r and so equation (9) also applies to method 3. 


The error in the evaluation of the interface temperature from equation (6) 
satisfies 


k k 
| (esas — €;) — + (e415 — ei) P| 2a 
(10) jan = 63,5 + Avr, a — EB, 
in 7 A%2C2 pz + AX C1 pr 


eo. By BY] 2At | 
eh =" LBae © Ax] LArece pe + Aries pr 


| ke b | 
(11) + 4 Len rg Fe 7. Any, | ™ + At’ B, + AtB;(Ax,*k, + Az,*ke) 


| Ave C2 p2 + Ax, C1 pi | " = 3(Are¢: pe + Axe: mr) 
AtB.(ky Ax;* + ke Ae’) 
12(Ave C2 p2 + Ax Ci pi)” 























Since in all three methods the Az used in the interface equation is min (At, , At), 


Atecepe , Aticipr ~ | A | 
At + At . Vin 
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and, therefore, (11) can be rewritten as 


i At’ B; ky Aaro( Aa’) + peste) | 
ein Say +> Bi + =I Ax kz + Ark, 





(12) 





24 Am ke + Areky 


Let Ax = max (Az, Azz), then a becomes 


Ax‘ 
(13) Ci i+] < aj = eB = By + Oy Re. 


Comparing ve (9) and (13), it can be seen that 
(14) ajay S a; 44 a  B, ref B+ B, where At = max (At, At.) 


and Ar = max (Az, Az). At any point ¢ = ng 


Ax’ 
(15) w Su +j[% B, + 4 aa. + x Bi}. 


The rate of convergence is, therefore, of the order 0(At & Az’). Since At is of the 
order 0( Az’), the rate of convergence is 0( Az’). 


5. Analytical Example. In order to examine the performance of the three meth- 
ods, a test case will be used for which some analytical solutions are known. The 
equations for the composite wall will first be reduced, by transformations of the 


variables, to the equations for a wall of one material. To do this we will impose the 
conditions: 


(i) kicipi = kecepe 
(16) 
(ii) a = 0. 
Define the transformations 
[o+a— 0%] 0Osrt384q 
(i) y= at 
(17) (ba + (1 — b)ae a4irsa@ 
(ii) rt = mbt OstSt, 
where 6 = kid, 


ky(a2 — a) + koa” 
This reduces equations (1) to (4) to 


Ou au 2 
(18) —=-, OS ySm,0 S37 S trrd 


u(O, r) = constant; r = 
(19) 
u(ad2, 7) = constant. r 2 0 


(20) u(y, 0) = constant; 0 < y < ae. 


+=: B, [ Bante ) + ke Ary (Axe }. 
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The solution to this set of equations is found in reference [6] for r = 0(.005).1, 
y = A, constant; = constant. = 0, constant; = 1, a = 1. 


6. Test Case Results. Each method was programmed for the IBM 704 EDPM. 
The test data used was k, = 10, q. = 5, p: = 2, ke = 5, = 4, pe = 5, ao = O, 
a; = .5, and a, = 1. Each method was run for three cases: Case A, Ar = .1; case 
B, Ax = .05, and case C, Ar = .025. 

To examine the rate of convergence, the maximum differences for a given time 
point were found between the results for case A and case C, and the results for 
case B and case C. The ratios of these maximum differences ranged between 3.8 
and 5.0 for each of the three methods. Since the ratios of the Az’s were 2.0, this 
would seem to corroborate that the rate of convergence is 0(A:’). 

To compare the accuracy of the three methods, the data presented in reference 
[6] was used. Their values correspond, according to the transformation presented 
in equation (17), to x = .55 and x = .7 for t = O (.01125).225. These values, 
as well as those obtained for case C for the three methods, are presented in Table 1. 
It can be seen from the table that, although they all showed close agreement, 
method 2 gave the most accuracy with a maximum of .08% error, method 1 the 
next with a maximum of .14% error, and method 3 the least with a maximum of 
.28% error. The symmetry of the transformed equation indicates that for this 
case the temperatures should be the same for « = .15 and x = .925, x = .3 and 
x = 85,2 = 45 and x = .775, and for = .55 and x = .7. When comparing 
the results at these points at ¢ = .1, method 1 has them all the same, method 2 
has a greatest difference of .04% , and method 3 has a greatest difference of .53% . 
A comparison at t = .225 shows method 1 has them all the same, method 2 still 
with a greatest difference of .04%, and method 3 with a greatest difference of 
.20% . These differences are reasonable in terms of the methods of choosing Ar for 
computation and illustrate that in method 1 the results for both materials is equally 
accurate, for method 2 the results for the second material is a bit more accurate 
than for the first material, while for method 3 the result for the second material 
is less accurate than for the first but their differences decrease as more time steps 
are taken. 


TABLE 1 











Time | Bxact | Methodi, | Methodi1, | Method2, | Method2, | Method3, | Method 3, 
Solution | #=.55 z= .7 z= 55 2=.7 z= 55 z= .7 
0225 | .9953 .99536 | .99536 . 99551 .99499 | .99579 . 99557 
.0450 | .9518 | .95145 | -95145 .95188 .95127 . 95254 . 95055 
-0675 | .8832 | .88261 | .8826! 88306 | .88275 | .88354 | .88136 
-0900 | .8088 | .80807 | .80807 | .80851  .80840 | .80868 | .80686 


.1125 | .7363 | .73556 73556 . 73598 . 73597 . 73590 . 73448 
.1350 | .6686 .66782 .66782 .66821 . 66826 .66797 .66689 
.1575 | .6063 .60561 .60561 .60597 . 60605 . 60562 .60481 
.1800 | .5496 .54891 .54891 . 54924 .54933 . 54882 .54821 
.2025 | .4981 .49739 .49739 .49770 .49780 .49724 .49678 
. 2250 .4513 45067 .45067 . 45095 .45105 .45047 .45011 
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Although the time increments were chosen on the basis of equation (7), it is 
pointed out in reference [6] that a slightly larger increment is possible, namely 


Almax & cp Aa" 


, 2k sin* [a a | | 





(21) 


Using case A, the programs were run until instability appeared in an attempt 
to see what the maximum increment actually was for two materials. For method 3, 
if the maximum Ai?’s are computed separately for each material of 5 laminae with 
equation (21), the result is At; = .005528 and Af, < .022112. The experimental 
results corroborated this since it was stable up to At; = .0055 and Ah = .0220 
but unstable for At; = .0056 and Af = .0230. For method 2, computing the maxi- 
mum time increment for each material of 5 laminae and then choosing the smaller, 
one gets At S .005528. The test showed the same result as it was stable up to 
At = .0055 but unstable for At = .0056. In method 1, when using equation (7), 
the adjustment of the thicknesses of the laminae in the second material led to 
equal Ai’s. However, when using equation (21), the A?’s obtained are Af, S .005528 
and Af. = .005125 since the number of laminae in each material is different. The 
test runs showed that it remained stable until At = .0053 but was unstable with 
At = .0054. This might indicate that a maximum time increment was used which 
is the average of Ai, and Af, but no conclusion is possible since the stability con- 
dition states that it should be stable below the computed Afmax but it is not neces- 
sarily unstable for a At above it. However, all these stability test runs seem to 
indicate that when there are two materials in a wall, the maximum usable time 
increment is quite closely related to the maximum increments computed for each 
material separately. 


7. Comparison of Computing Time Required. The amount of computing time 
required for each method can be compared by comparing the number of tempera- 
tures that must be evaluated. 

Let us assume that the thicknesses of the first and second materials are equal 
and that ty = PRAt, (where P is any integer). It should be noted that given Az, 
r;, and re, each method will compute the same value for Aft; and Az, . 

In method 1 there are N,(1 + ~/R) laminae and PR time points. Therefore, 
the number of temperatures computed equals P[RN,(1 + ~/R) — R]. In method 
2 there are 2N, laminae and PR time points and so P[2N:R — R] temperatures 














are computed. In method 3 each of the N; — 1 temperatures in the first material 
TABLE 2 
Case A | Case B Case C 

Method 1 P (56) | P (116) | P (236) 
Method 2 P (36) P (76) P (156) 
Method 3 P (30) | P (55) P (105) 

Ratio | —. 

M1:M2:M3 1.87:1.20:1.00 | 2.11:1.38:1.00 | 2.25:1.49:1.00 
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is computed PR times, the interface is computed PR times, and for each step of 
At, , which occurs P times, N, — 1 values are computed plus the additional number 
R-1 


of transition values = >> i. Method 3, therefore, computes 
i=1 


P| NR +™ — 145-2] 
values. From this it can be seen that for all R > 1 method 2 is faster than method 1. 
For 2N, — 2 > R > 1 method 3 is faster than method 2. This comparison is il- 
lustrated in Table 2 for cases A, B and C of the test data used. 


8. Conclusions. The analysis and test cases used considered constant boundary 
and initial conditions. Since the stability and convergence depend also on the 
boundary and initial conditions, as has been pointed out in references [1], [3], and 
[4], it is quite possible that the introduction of varying conditions would lead to 
different results as to the usefulness of each method. 

From this study it seems that method 1 is the least acceptable since it takes the 
most computing time, gives less accuracy than method 2, and presents the most 
difficulty when R is not an integer. Depending on the amount of accuracy desired, 
methods 2 and 3 seem of equal usefulness since method 2 gives the most accuracy 
but method 3 takes less computing time. 


9. Acknowledgements. The author wishes to thank Dr. Peter Henrici for his 
suggestions concerning this study. 
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A Method of “Alternating Corrections” for the 


Numerical Solution of Two-Point Boundary 
Value Problems 


By David A. Pope 


Abstract. In this paper a method of “alternating corrections’ is defined and 
analyzed for the numerical solution of the two-point boundary value problem 


y” = f(x, y) 

(0.1) y(0) =a 

y(1) = b. 
The case where the first derivative does not enter explicitly into the differential 
equation is chosen for simplicity of treatment. The alternating corrections method 
can easily be modified to treat the more general case. The function f(x, y) is as- 
sumed to have continuous second derivatives, but the differential equation may, 
of course, be non-linear. 

The method to be described is essentially a relaxation technique suitable for 
an automatic digital computer. The main feature of the method is that most of 
the “correcting” is done in the early stages of the computation, using a small 
number of points; thus a rough approximation to the solution is obtained quickly. 
This approximation can then be made more accurate in the later stages of the 
computation, as the number of points is increased. 

In Section 1 the method is described. Section 2 gives a rigorous truncation and 
stability analysis. Section 3 contains the proof of the convergence of the method 
giving an estimate of the rate of convergence, and in Section 4 some experimental 
results obtained on a digital computer are examined. 


1. Definition of the Method. In the following, we will denote by R a closed 
and bounded region of the x-y plane, in which we will assume both the solution to 
(0.1) and the approximations to that solution are known to lie, a priori. (See Collatz 
[1] p. 188 for a sufficient condition for the existence of the solution to this problem.) 
The function f(x, y) is assumed to be continuous and to have continuous first and 
second derivatives in R. The method to be discussed consists of two stages, as 
follows: 

A. Interpolation by Halves. Suppose the interval [0, 1] is partitioned into n 
equal parts by n + 1 equally spaced points, and an approximation y; to the solu- 
tion of (0.1) is defined at these points. We then refine the partition by subdividing 
[0, 1] into 2n equal parts, by 2n + 1 equally spaced points x; ,7 = 0, 1, 2, --- 2n. 
Then the points of the original partition are given by the x; with even index. We 
then interpolate y; for the odd indices by using the explicit formula 


yi = » iyi + ya) — 4 [f( visa, Yi) + (aja, yd), 
(1.1) - 4 
3 = 1,3,5, --- 2n — 1. 
Received May 3, 1960. 
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We shall abbreviate the right-hand side of (1.1) by the operator K(y;). Here, and 
in the following, h = Ax; = 1/2n. 

B. Alternating Corrections. After step A, the values of y; with odd and even 
index are corrected alternately by using the same formula (1.1). Thus we have 


(1.2) ys’ = K(y;’) for j odd, 
(1.3) ys; = K(y"s’) for j even, j] * 0,7 # 2n, 
(1.4) y’.' =a, 
(1.5) ys. = 5, 


where y;° is the value of y; at the s-th iteration of step B. Then, as we shall show 
in Section 3, as s — «, the values y;° approach the solution of the system of differ- 
ence equations 


yi = K(y;) j = 1,2,---,2n — 1, 
(1.6) yo=a 
Yon = b. 


To start the computation, we usually will set n = 1, and yo = a, y, = b. Then 
step A interpolates a value at z = 3, and we renumber the values yo , 7: , y2 . Here 
step B is not needed, so we perform step A again, getting now five values. At this 
point we perform step B a number of times, until sufficient convergence to (1.6) 
for our purpose is obtained. We continue in this way, doubling the number of points 
with step A, then following this with a number of iterations of step B, until we have 
the desired accuracy. In Section 3 we will consider some estimate of the number of 
iterations of step B necessary for a given accuracy. 


2. Stability and Truncation Error. In this section we shall give a rigorous esti- 
mate of the truncation error in and stability of the difference equations (1.6). We 
note that the global truncation error is of order h~” times the local truncation error, 
rather than h™ times, as might be expected from a naive analysis. 

In the following, we let Y(a) be the exact solution to the differential equation 
(0.1). Let Y; = Y(a;), and let y; be the exact solution to the system of difference 
equations (1.6). Let the error e; = Y; — y; . Then, using the law of the mean, we 
obtain from (1.6) the system of difference equations 


1 h’ 
3 < lejss + ej-al = 7 lei Sy (tin, ni+1) 


(2.1) 
+ ea fy(aja, a)) + &, j=1,2,---2n—1, 
and 
(2.2) G@ = & = 0, 
where 7; is between Y ; and y;. Here the local truncation error 
~24 
j= —- a y""(&), Ta < § <u, 
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which is obtained from Taylor’s formula. Rewriting (2.1) in matrix form, we have 


2 
(2.3) Ae = —t re+t 
where ¢ is the column error vector whose transpose is (¢ , €2 , ~*~ €zn-1), ¢ the trun- 
cation error vector with transpose (t; , 2 , --~- fens), A is the second difference ma- 
trix 
- 1 - 
= abe 
5 0 0 
1 1 
Hap Qidwuomg ¥ 
1 1 
0 -~- 1 —= 
2 2 
ne" 4 








and F is the matrix 


i g2 0 O 
gu 0 923 0 


F= 
i 932 0 933 


where gi; = f,(zi, 7;), evaluated at the intermediate points given in (2.1). It is 


well known that the matrix A has eigenvectors v; with components sin . m=1, 


2, --- 2n — 1, and corresponding eigenvalues \; = 1 — cos x j=1,2,--- 2n-—1. 


_1 
Hence A exists, and its largest eigenvalue is (1 — cos z) . Therefore, multiply- 


2n 
ing (2.3) by A”, we have 
2 
(2.4) (1 ~ 2 A“P) e=A't. 


Now we can prove two lemmas giving estimates of the error. 

Lemma 1. Suppose f, > 0 in the region R. Then the components of the error vector 
e satisfy the inequality 
5h? max | y" | 


jl< 
(2.5) lels 12 min f, 


t= 1,2,---2n—1, 


where the extreme values are taken over the region R. 
Proof. Let the norm ||e || = max;|e;|, and the subordinate matrix norm 
\| A || = max; >>; | a; |. (Cf. Faddeeva [2] p. 58.) Rewriting (2.3) in the form 


(2.6) (4+hr)e=U-Bye=t, 


2 


we have defined the matrix B = J] — A — . F. But since f, > 0, we can conclude 
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that || Bi] s 1 — : Min f, < 1, where the minimum is taken over the closed 

region R. Hence the series J + B + B’ + --- converges to (J — B)™, 

and || (J — B)* || < (1 — || B||)~. Also we note that || t|| < Sh pmax | y* |. 
Using these estimates and (2.6) gives us 


pitstieas 
* h? min f, 


1 | || 5h‘ iv 
Nell $ py = ag ml" | 


from which inequality (2.5) follows. 
The second lemma takes care of the case when f, is negative or zero in R. Here 
we estimate the root-mean- square of the error. 


Lemma 2. Suppose | f,| < 2° in the region R. Then, for h sufficiently small, the 
error vector satisfies 


lle 5h’ max | y"” | ‘ 
2.7 — = — + O(h*). 
(2.7) V2n —1  12e°(1 — x ‘max |f, |) ' 


Proof. In this proof, we use the euclidean norm, and the subordinate matrix 
norm ||A|| equal to the square root of the largest eigenvalue of AA” (Faddeeva 
[2] p. 59). We see by inspection that FF” has a maximum row sum not exceeding 
max 4f,’. But we know that its largest eigenvalue does not exceed this maximum 
row sum. This gives the estimate 





(2.8) \|F\| < 2 max |f,|. 

Also, from the eigenvalues of the symmetric matrix A we obtain 
(2.9) \A|| = (1 — cos xh)”. 
Using this together with (2.8) we get 


h? max | fy | 


IE spl 


| 


But for small h, we have the estimate 


lA 


2 


—2 
(2.11) (1 — cos zh) ah 


+ O(n’), 





lr) 


hence, for h sufficiently small, 
| Rh? | a 
(2.12) = AF < (x * + O(h’)) max | f,| <1, 
since max | f, | < 2. Also, with this norm, ||¢|| < mn Max | y” | </2n + 1. 


: er ® ely 
Therefore, by the same reasoning as in Lemma I, (7 + r A 'r) exists, and its 
norm does not exceed (1 - . A“F 


get (2.7), proving the lemma. 





-1 
) . Putting these estimates into (2.4), we 
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It should be noted that the restriction —f, < 7° is a natural one, as can be 
seen from an examination of the boundary value problem 


y” = —Ky 
(2.13) y(0) = 0 
y(1) = 1. 
Here, of course, the solution does not exist when K = —f, = 7, as we have an 


eigenvalue problem. 


3. Convergence of the Alternating Corrections Method. In this section we shall 
give a proof of the convergence of the alternating corrections method (method B 
of Section 1). Following the notation of Section 1, we shall denote the value of 
the approximation at the point x; for the s-th iteration by y;’, and the exact solu- 
tion to the difference equations (1.6) by y; . Then we define the error €;° = y; — y;’. 
Then we have, for 7 odd, 


« \ 8 1 & h? & 8 . 
(3.1) ” 3 lej41 ma éj—1] —~ 4 lejaa Syl tin, nin) + G-ify (aja, m-)), 
where 7; is between y; and y;*, for each 7. For j even, 7 # 0,7 ¥ 2n, 


\ 8 1 8 8 h’ e 8+ } 
(32) g= x ei + ea) — ry (e551 fy (ajar, Ein) + Gti Sy (aja, &-4)) , 


and for the endpoints, 


(3.3) € = €2n = 0. 

4 ; i ’ . : 
Now we define » = max | 1 — 5 Sy , the maximum being taken over the region 
R. Then for 7 odd, we have 
(« s+1 8 s 
(3.4) €j S ull ej41| + | Ga 


and for j even, 7 # 0,7 ¥ 2n 
9. s+l | s+1 s+l1 
(3.5) ej | S dull esta | + | ej |). 


To estimate the error ¢;', we majorize it with a quantity /;*, defined recursively 
as follows: 


(3.6) E? = | 4; | j= 0,1,---«. 
(3.7) ES" = WwiEja + 2Ej° + Ein) = jf = 1,2,---n -— 1. 
(3.8) Eo" = Ex** = 0. 


2n 


In this section we will use the euclidean norms |! « || * = e; and 
} J 
0 


m 
7 12 12 
|E\ = > E}. 
0 


Lemma 3 gives an estimate for || ¢ || in terms of || 2 





ee 
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Lemma 3. For every s = 1 


(3.9) Ne? s (1+ 0°) Bll’. 


Proof. First the inequality | @;| < £;’ is established by induction on s. For 
s = 0, the inequality holds by definition. Assuming the inequality is true for s, the 
proof for s + 1 follows from the inequalities 


et? | < 4ufl| ef. | + | ath |] 
S 4p" || 5-2] + 2) ds) + | Sissel 
< RE 
<n. 


Now for 27 + 1 we have the inequality 


€2j+1 > 5ul| €2j j + | €2 542 \] 
2ulE ;* + E54). 


Combining these, and using the triangle inequality, we obtain 


lA 


s 2 — +a {2 1 = 78 +8 {2 
€ SLES P+ pe Dh Ei + Ein! 
(3.10) 


? 


(1+) |E |’, 


IA 


which proves Lemma 3. 
We now expand £;;’ in a finite Fourier sine series, with Fourier coefficients given 
by 


9 n—l . 
ey 13. Mr - 
(3.11) F,’ =- DE; sin 2 m=1,2,---n—1. 
nN j=l n 
Then we have the expansion 
n—l mj 
(3.12) E;' = a FF, an —, j=0,1,---,n. 
m=1 n 


Substitution of (3.12) into (3.7) now yields the recursion relation 
‘ a m@\ 1, « 

(3.13) ry, = 5M 1 + cos y Po 

From (3.13) we then obtain the estimate 


(3.14) E'\i\< E iw (1 + cos *)| ID 


and applying Lemma 3, we have the final estimate 


(3.15) lis (1 4+,°)" E w (1 + cos *)| Pa 


Therefore we have proved 


Lemma 4. If —1 < p= wi (1 + cos *) < 1, the alternating corrections method 
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will converge geometrically to the solution of the difference equations (1.6), with con- 
vergence factor p. 
For small h, we note that 


(3.16) p = max {i — h'(f, + #) + O(h')] 


and again we have the natural restriction for convergence mentioned at the end 
of Section 2. 

The number of iterations needed at each stage in the method can now be esti- 
mated as follows. Clearly, the convergence factor increases as h — 0, hence the 
convergence is much faster for large values of h. On the other hand, it would be 
futile to carry the iterations so far that || « || is much smaller than || ¢ || , as we are 
interested not in the solution of (1.6), but of (0.1). Hence a useful compromise 
might be to iterate the alternating corrections until || « || and || e || are approxi- 
mately equal, then to interpolate, and start again with interval h/2. If this scheme 
is followed, we would want to cut the error « by a factor of about } by iteration 
after each interpolation, since the error || ¢ || is of order h’. Then we have p’ = 3, 
which gives us the approximation 


ape log4 ss dog 4 
log p W(fy + 7) ; 


This shows that it would take about four times as many iterations for the next 
stage, after interpolation by halves. Since there are about twice as many points, 
the total amount of computational work is multiplied by eight at each succeeding 
stage. Clearly this process cannot be used for very many stages. 

In practice h will probably not be made less than 2~* or 2~°, and if more ac- 
curacy is needed, a more sophisticated set of difference equations than (1.6) would 
be used. The alternating corrections method, however, is excellent for obtaining 




















TABLE 1 

h } error atx = 4 number of iterations 
21 | 052083 0 
27 .013021 27 
2* .003255 100 
2-4 | .000814 | 329 
2-5 | .000203 | 1026 
2* | .000051 | 2948 

TABLE 2 

h | error atx = 4 number of iterations 
ye | .065074 0 
2? | .013935 23 
a .003366 82 
24 | .000834 | 269 
2-5 .000208 | 831 
2 | .000052 | 2339 
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a good approximation quickly, which could be used as a first guess in a more com- 
plicated relaxation scheme. 


4. Some Experimental Results. In this section we shall discuss the results of 
two problems which were computed using the alternating corrections method. The 
compuiation was done using the Univac Scientific 1103 computer at the University 
of Minnesota Scientific Computing Laboratory. 

The first problem was the linear equation 


y” = 2x 
(4.1) y(0) = 0 
y(1) = 1. 


In Table 1 the results of this computation are summarized. 


Formula (2.7) with f, = 0, y” = 4 gives the r.m.s. error < .16887h’, in good 
agreement with the error at x = 3. 
For each value of h, stage B was iterated until there was no change larger than 
2°** in any y;’. The number of iterations necessary to accomplish this is also given 
in Table 1. If we use the estimate || ¢ || = .16887h, we get the relation 

16887h'p* = 2”, 
from which we get the approximation 
s & h*(1.856 + .2026 log h), 


which agrees well with the number of iterations actually performed. 
The second problem tried was the nonlinear equation 


y” = 2y 
(4.2) y(0) = 0 
y(1) = 1 


Table 2 gives a summary of this computation. 

The fact that the global truncation error is of order h’ is again displayed in 
Table 2. The number of iterations necessary at each stage was governed by the 
same scheme as in problem 1, but with the criterion 2™ instead of 2~**. 

Finally, it may be noted that the approximations generated by the alternating 
corrections method could be improved greatly by a “deferred approach to the 
limit”, using the approximations obtained from the last two values of h computed. 

The author wishes to thank Mr. James Rude for his helt in preparing the 1103 
coding. 
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TECHNICAL NOTES AND SHORT PAPERS 


Note on the Asymptotic Expansion of the Mod- 
ified Bessel Function of the Second Kind 


By E. Dempsey and G. C. Benson 


Tables of modified Bessel functions of the second kind, K,(z), for integral order 
n and for special forms of real positive argument z have recently been computed 
|1]. To obtain the desired accuracy throughout the tables it was found necessary 
to look into the methods available for approximating the remainder term R, in 
the asymptotic expansion of the function: 


p-l 
(1) Kn(z) = >. Um + (Rp/Up) Up 


m=0 


where 


—~ 


m wan (ZBatmeD & 
sal 22) T(n — m + 4) m\(2z)™ 


Ups is usually the least term (in absolute value) of the series. Since the converging 
factor R,/u, may itself be developed as an asymptotic series, it is convenient to 
discuss this quantity rather than the remainder term. 

In a treatment due to Dingle a general theory of asymptotic expansions and 
their converging factors [2] is applied to the hypergeometric function [3], special 
cases of which are the modified Bessel functions of the first and second kinds. The 
formula relevant to the latter case is as follows 
(3) Rp _ ST(n+})r(n+p+}3- 8) (92)! 4, 4(22) 

Up = Tin + 3 — dT (n+ p + 3) 
in which u, is now not necessarily the first term after the least term of the series in 
Eq (1) but may be a term some way beyond this. Going to the fifth term in Eq (3) 
the functions AS”, (2z) are defined by 





As +0) = 5-2 (1-%) + (1 — 29 — 46°) 
(4) 
+ jogs (I + 60 + 86 +80) — soa (13 + 226 + 206° + 160° + 166*) 
a) dg: dis, 2 
As (s+0)= 7 iga (1 + 40) + gaa (3 + 80 + 126) 
(5)* 
l 2 ; 
== ate 9 9A 299° 
agg (11 + 206 + 246° + 326") 
. (2) eee | =. a us - > Ag 
(6) AY?(s+ 0) = aa + yqp (8 + 30) - aaa (8 + 126 + 246°) 





Received November 13, 1959; in revised form, February 3, 1960. Issued as N.R.C. No. 5711. 
* The coefficient of (64s*)-! in A given in ref. [3] is incorreet; the correct form is given 
above. 
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1 
7) a® nm 
(7) (s + 0) vie — 34 “1 + #) 


(4) ye = 
(8) A%(s+6) = — 5, 


These equations are derived on the basis that @ is small. @ and p are obtained from 
the relationships 


(9) p-n-4$=8 
(10) 2z=s+0 


By restricting the value of 6 to | @| < 0.5, s is forced to take on the half integral 
value nearest to 2z, since both p and n must be integral. Using this method the best 


results were obtained by calculating K2(z) and K;(z) and using the recursion for- 
mula 


(11) Khe) «© Kesha) + = K,.(z) 


to calculate the K,,(z) required. In the program for calculating K,,(z) a convergence 
test was applied to each term of the series in Eq (3); the summation was terminated 
at the least term and, since the terms alternate in sign, one half of the following 
term was added to the sum. 

An alternative treatment due to Burnett [4] starts from a representation of 
R,/u, as a double integral which is subsequently reduced to a series of products of 
gamma and hypergeometric functions. Expanding this series leads to the following 
expression for the converging factor, which in this case is applied to the term fol- 
lowing the least term of the series for K,,(z) 


1 1 1 —e ; 
Uy 5+ See eee 


oo — n)(30? — 11en + 9? — 36 - 3) +e -35-3] 

(12) 
1 3 92,2 2 a7»? s9.8 » D2 
z~s 4a — 280'n + 48on — 27n — ilo + Ilon—- sn 


27 3 ° 7 4 
— llo + 1l6n + e) +o —o — ot 3° + aI 


ra or 
bol Or 





where 
(13) o=2z2—p+n-— }. 


Previously published versions of this series [4], [5] go only as far as the term in 2~ 
It is easy to show that the successive powers of z' given in Eq (12) may be derived 
completely from the first five terms of Eq (3); in other words, to the extent that the 
equations appropriate to them have been given, the two methods should be roughly 
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. 2.—Deviation plot for calculations of K,(z) 


equivalent. In applying the Burnett treatment the converging factors for Ky and 
K, were calculated from Eq (12) and functions of higher order obtained by use of 
the recursion formula. 

The accompanying figures show graphs of logi| «| against the argument z 
for Ko(z) and K,(z) calculated using the Dingle (D2) and Burnett (B) methods, 
where 


pa K(exact) — K(calc.) 


(14) K (exact) 





The values for K(exact) were taken from a paper by Aldis [6] which gives values 
of Ko(z) and K,(z) in general accurate to more places than are required here. 
Curves S were obtained using Eq (1) with R,/u, set equal to 0.5. The parts of 
the curves below about | «| = 10°” should be ignored since they are influenced 
by the fact that the computer (a Bendix G15D) carried only twelve figures in the 
double precision floating point mode used. 

It is apparent from Figures 1 and 2 that the Dingle method is somewhat better 
than that of Burnett. This may be due in part to the fact that the Dingle equations, 
as given here, contain an infinite number of terms which may be looked on as being 
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incomplete Burnett type coefficients of powers of z* higher than the fourth. Above 
about z = 10 it is also evident that for a computer carrying only twelve figures 


there is nothing to be gained in using a more elaborate converging factor than 
R,/u, = 0.5. 
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On the Factors of Certain Mersenne Numbers 


By John Brillhart and G. D. Johnson 


1. Introduction. For the past 10 months the authors have been conducting a 
search for factors of certain Mersenne numbers on the IBM 701 at the Computer 
Center, University of California, Berkeley. The following is a report on the nature 
and results of that search. 


2. Extent. Prime factors g were sought for the numbers M, = 2” — 1 for primes 
p < 1200 in the intervals indicated: 


p = 101 VPegct 
103 <p 5157, p# 115i << 
157 < p < 257 l<q<2" 
257 <p < 1021, p¥ 397 l<q<2 
p = 397 l<q<2” 
1021 < p < 1200 l<q<2 


No factors <2™ were examined for 101 < p < 157, since these had already 
been investigated [1]. No M, were examined for p < 101 or p = 151, since these 
numbers have presumably been completely factored. Possible factors <2” were 


Received May 5, 1960. 
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also investigated for Mss; , the Mersenne number whose exponent is the “last” 
Fermat prime. M;; was investigated to 2” in the hope of finding more small factors. 


3. Results. 


A. Fifty-five new prime factors were discovered, 6 of which for M, below the 
traditional “limit”? p = 257. These factors are given in the accompanying table, 
and are indicated by *. Also included are all published prime factors, and 6 new 
ones (indicated by t+) of E. Karst, Brigham Young University. Thus, the table is 
believed to be a complete listing of all prime factors of M, for p < 1200 known 
at this time. No factor was found for M53; , whose character is still unknown. 
Since no factor was found to M1 below its cube root, it is the product of two primes. 


B. All known prime factors of M, ,n < 10 000, were tested and found correct, 
with the exception of the two misprints in H. Riesel [2], as noted earlier by J. 
Selfridge [3]. In addition, all factors were tested for multiplicity, but no new mul- 
tiple factors appeared. Hence, to date, only a few multiple factors are known for 
composite exponents n, while none have been found for prime exponents, further 
supporting the conjecture that none exist. 


4. The Program. 


A. Srructure. If d| M,, then d = 1 (mod 2p). Also, since 2 is a quadratic 
residue of M, , n odd, then d = +1 (mod 8). Thus, the divisors, d, lie among the 
common terms ¢, of these arithmetic sequences. 

In production these terms were generated consecutively by the repeated use of 
an increment table, which had also been constructed to produce no terms divisible 
by 3, 5, 7, or 11. (See [1].) 

Divisibility of M, by each t, was tested by examining the remainder of V, 
(mod t,) for 0. 

For 101 S p S 223, M, was reduced mod ?¢, by multiple precision division. 

Example 1. The remainder of Mj mod ¢, was computed for each ¢, by 3 divi- 
sions, until ¢, was >2”, at which time an initial dividend of 67 binary places could 
be used. This change, which produced the remainder in only 2 divisions, was 
actually introduced when ?, was >2* by using a modulus of 2*t,,0 < a S 3, 
instead of #, , the error in the final remainder being removed after the last division 
by an appropriate number of subtractions of ¢, , or multiples of ¢, . This device 
was used consistently in all routines whenever possible. 

When the program was first run for p 2 223, the final remainder was computed 
by residue methods consisting of successive squarings and doublings of the residue 
of some initial power of 2, followed by a subtraction of 1. Later it was realized, 
that in a double register machine like the 701, a residue between the initial and 
final residue could usually be multiplied by a power of 2 greater than the first 
without producing an illegal divide condition in the registers. The magnitude of 
the power that could be used was found to depend on the length of the registers 
(35 binary places) and the length of ?, . 

This discovery decreased the testing time for each ‘,, by about 30% , but greatly 
complicated the programming, since from the many possible programs, one had 
to be chosen that required a minimum number of machine cycles. 














19 


29 
31 
37 
41 
43 
47 
53 


61 
67 
71 
73 


79 


89 

97 
101 
103 
107 
109 
113 


127 
131 
137 
139 
149 


151 
157 
163 
167 
173 
179 
181 
191 
193 
197 
199 
211 
223 
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TABLE or Factors 





3 


7 
31 
127 


| 23-89 


8191 

131071 

524287 

47 - 178481 

233 - 1103-2089 

2147483647 

223 -616318177 

13367 - 164511353 
431-9719 - 2099863 

2351 - 4513 - 13264529 

6361 - 69431 - 20394401 

179951 - 3203431780337 
2305843009213693951 
193707721 - 761838257287 
228479 - 48544121 - 212885833 
439 - 2298041 - 9361973132609 
2687 - 202029703 - 1113491139767 
167 -57912614113275649087721 
618970019642690137449562111 
11447- prime 


prime 

745988807 - 

3391 - 23279 - 65993 - 1868569 
-1066818132868207 

prime 

263 - 


1812] - 55871 - 165799 - 2332951 - prime 
852133201 - 

150287 - 704161 - 110211473*- 
2349023 - 

730753 - 1505447 - 

359 - 1433 - 

43441 - 1164193 - 7648337 *- 

383 - 

13821503* 

7487 - 


15193 - 
18287 - 196687 - 1466449 - 2916841 - 





? 





271 
277 


293 
307 
311 
313 


373 
| 379 
383 
389 
397 








| 

| 

| 401 
| 409 
| 419 
| 421 
| 431 


| 433 
| 439 
443 


_- > 
ao on 
— = 


va 
435 


— 
* 
= 


TELL 


Pactors 


1504073 -20492753* - 

1399 - 135607 - 622577 - 

479 - 1913 -5737 - 176383 - 134000609* - 
22000409* - 

503 -54217- 


| 23671- 


13822297 * - 


1121297 - 


9623 - 


14608903 * - 85798519* - 
5344847 - 

10960009* - 

9511- 


18199 - 2806537 - 


931921 - 

719 - 855857 -778165529* - 
12479-51791041*. 
25569151 *- 


1440847 - 

56478911*- 

2383 - 6353 - 50023 - 53993 
- 202471 -5877983t - 


839- 


863 -3449 - 36238481 * -76859369* 
- 558062249 * - 


104110607* - 
887 - 

1256303 - 
150327409* - 
2767 - 
11113-3407681f - 
121606801 *- 
33385343 * - 
4871- 

983 - 7707719} - 
20959 - 
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TaBLeE or Facrors—Continued 














? Factors p Factors 

503 839 | 26849- 

509 | 12619129f- 853 

521 | prime 857 6857 - 

523 859 | 7215601- 

541 863 | 825891i-169382737*- 

547 | 5471- 877 | 35081-1436527*- 

557 | 3343-21993703*- 881 | 26431- 

563 883 | 8831-63577*- 

569 | 15854617*-55470673*- 887 | 16173559*- 

571 | 5711-27409*- 907 | 1170031- 

577 | 3463- 911 | 1823-26129303*- 

587 | 554129-2926783*- 919 

593 | 104369- 929 | 13007- 

599 937 | 28111- 

601 | 3607-64863527*- 941 | 7529- 

607 | prime 947 | 295130657*. 

613 953 | 343081- 

617 | 59233- 967 | 23209-549257*- 

619 | 110183- 971 

631 977 | 867577-1813313*- 

641 | 35897-49999*- 983 

643 | 3189281- 991 

647 997 

653 | 78557207*-289837969*- 1009 — 3454817- 

659 | 1319- 1013 | 6079- 

661 1019 | 2039-75407* 

673 | 581163767*- 1021 | 40841-795808241*- 

677 1031 | 2063-435502649*- 

683 | 1367- 1033 | 196271-36913223*- 

691 1039 | 5080711- 

701 | 796337 -2983457*-28812503*- 1049 | 33569-459463*- 

709 | 216868921*- 1051 | 3575503- 

719 | 1439-772207*- 1061 | 

727 1063 | 

733 1069 | 

739 1087 | 10722169*- 

743 | 1487- 1091 | 87281- 

751 1093 | 43721-111487*- 

757 | 9815263 -561595591*- 1097 | 980719-4666639*- 

761 | 4567-6089*- 1103 | 2207- 

769 1109 | 

773 | 6864241 -9461521f- 1117 | 53617- 

787 1123 | 

797 1129 | 33871- 

809 1151 | 

811 | 326023- 1153 | 267497- 

821 | 419273207*- 1163 | 

823 1171 | 

827 | 66161- 1181 | 4742897- 

829 | 72953- 1187 | 256393-113603023*- 
1193 








| 121687- 
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In some cases, the initial residue was produced from a comparatively small 
power of 2 by a single division, while in others, it was obtained from a fairly large 
power of 2 by multiple-precision division. 

Example 2. For Mj; , 4 different programs were used, each improving on and 
replacing the preceding, when the length of ¢, permitted. The first divisor used 
was t; = 3-794 + 1 = 2383, which also happens to be the first factor. This is 
shown below by the calculation schemes of the 4 programs, although only the 
first was actually used to test such a small possible divisor. With each scheme is 
also given the interval of ¢, , for which it was used. The letters ir after a residue 
indicate the initial residue used by the squaring part of the routine. 











I:1 < & < 2%. | II: 233 < tn < 2". Ill: 27 << 2. IV: 2” < & < 2@. 
| | 
2% = 1792 (mod 2383) 20 = 1657 (mod 2383)|2% = 1862 (mod 2383)/2 = 299 (mod 2383) 
26° = 1657 25 = 342 ir \2°7 == 1368 ir \2 == 706 ir 
2% = 342ir 2 = 197 [214 = 769 (2198 == 389 
21 = 197 [2196 = 693 [2197 = 1386 29% = 1192 
21% = 693 i208 = 1266 [24 = 298 7 = 86 
2392 = 1266 27= 1 2 = = 
9397 = 1 | 





B. Propuction. The program was run for 60 hours, each p < 223 requiring 
approximately 23 minutes, and each p = 223 requiring from 8 to 18 minutes, the 
larger exponents taking progressively less time. The special number Mj, was run 
for 10 hours. 

The routines used are believed to have been accurate, a fact which will be 
ascertained at a future time, when more rapid computers will accomplish in a few 
minutes, what has now taken many hours. 

The authors would like to thank Professors R. M. Robinson and D. H. Lehmer 
for their suggestions and ideas, and also the staff of the Computer Center for their 
many courtesies. 
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Further Evaluation of Khintchine’s Constant 


By John W. Wrench, Jr. 


In his fundamental investigation of the metric theory of continued fractions 
Khintchine [1] proved that the limit, as n tends to infinity, of the geometric mean 
of the first n partial quotients in the simple continued fraction expansion of almost 
all real numbers is the absolute constant 


. = 1 In r/in 2 
. -HO+445) | 


A different proof, by C. Ryll-Nardzewski, has been recently reproduced by 
M. Kae [2]. 

The numerical evaluation of Khintchine’s constant was considered by D. H. 
Lehmer [3]. In addition to finding an approximation to K to 6 decimal places, 
whose accuracy was subsequently discussed by D. Shanks [4], Lehmer investigated 
the geometric mean of the first one hundred partial quotients of z. 

Recently R. S. Lehman [5] computed the first 1986 partial quotients of z on 
ORDVAC in order to test the applicability of a similar theorem of Lévy [6], which 
asserts that, as n tends to infinity, the nth root of the denominator of the nth 
convergent tends to exp (2°/12 In 2). 

Shanks and the writer [7] have studied the representation of K by infinite series 
and by definite integrals. The computational effectiveness of these series was illus- 
trated by the evaluation of K to 65 decimal places. This calculation has now been 
extended by me to 155 places, using the same series as previously, namely: 


oS 3 f1 AS, 1 A Sx Su, | 
In2inK = In5+In2Ing meets et whe : 


yf? 
where S.; represents 


yn by -1- 7%. 


A preliminary step in this calculation consisted of the formation of a table of 
¢(2k) to at least 155D for k = 1(1) 257. The first 60 entries of this table were 
computed by the formula 


k—1 Bax. ( 2m )*" 
22k)!’ 

where the notation for the Bernoulli numbers is that used by K. Knopp [8]. The 
numerical values of these numbers were taken from the tables of H. T. Davis [9]. 
The requisite decimal approximations to x” /(2k)! were obtained from my manu- 
script table [10] of such data. The remaining entries were computed directly from 
the series defining ¢(2k), a maximum of eighteen terms being required initially. 

From these values of ¢(2k) the approximations to So, and Sox/k were then 
computed to 155D. All these data were subjected to the following check relations: 


§(2k) = (-1) 
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> (r(2k) —  =3, 
k=l 4 


Su/k = 
dd Su/k = 35 

Substitution of the computed values in these formulas resulted in discrepancies 
ail less than 3 units in the 155th decimal place. 


The final results of this calculation when rounded to 155D are as follows: 


In2inK = 

0.68472 47885 63157 12329 91461 48755 77762 04606 75416 33744 
88366 06289 86781 59568 82176 26936 10437 07681 43495 85810 
09970 15696 93974 12470 41578 92227 14396 39612 78766 18097 
72947 ---, 

In K = 

0.98784 90568 33810 78966 92547 27147 07295 43261 99254 96088 
67354 27755 30068 72109 27094 18512 90938 20768 83372 75259 
67479 51231 68801 78544 35925 75519 06227 59695 60965 06769 
43483 ---, 

K = 

2.68545 20010 65306 44530 97148 35481 79569 38203 82293 99446 
29530 51152 34555 72188 59537 15200 28011 41174 93184 76979 
95153 46590 52880 90082 89767 77164 10963 05179 25334 832259 
66838 ---. 
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Note on [ ” ee (") Ji (2) a” dx 


By Henry E. Fettis 


In a recent paper this integral has been calculated numerically by Weeg [1] for 
the cases n = 0, 1. The following analytic expression for this integral when n = 0 
in terms of elliptic integrals is given by Byrd and Friedman (([2], formula 563.01): 


I(p, r, 8) = [ eJo(rt) Ji(st) dt = au — Ao(8, k)] 
(1) 
1 
8 


2 / / / 
[2 “gf {E(k) F(6,k’) + K(k) (E(6,k’) — F(8,k my], 
with 
_W-p-rt+eN-pt+r-s 
Ptr ae +e Fr te 
PS PPA pt+rt+s 
sin 8 = rh + recs at! 
N = V(p? + 1 + 8°)? + 4p's?. 

Ao(8, &) is known as Heuman’s Lambda Function. 

Calculations based on the above expressions do not agree with the ones made 
numerically in [1], and since the latter values have been verified by an independent 
relationship given in Weeg’s paper,* the author was led to attempt to verify the 
above analytic expression independently. A closer examination of Eq. (1) leads 
one to question its validity by consideration of various special cases for which 


simpler expressions can be found. For example, it is known ((3], Chapter III, 
Art. 7) that 


ke 








(2) U(p,0,s) = VE +# =? 





whereas (1) fails to reduce to this simple form when r = 0. It is also possible to 
demonstrate directly the relations: 


a . ee ll 
(3) sl(p,r,8) +ri(p,s,r) = 1 wttrte 1 ee Q.4(x) 
where 

kK’ = 4rs/[p + (r + s)’], r=(p +r + s’)/2rs 


and 
(0, r,s) = 1/s, r<s, 
(4) == 1/20, fr = 8, 
= (), r> 8, 





Received December 8, 1959; revised April 28, 1960. 
1 
* Eq. 2 of [1] should read I(m, £,0) + I(n, 1) = ove (Q(z) — 1Q(z)). 
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neither of which can be obtained as a special case of (1). In fact, for p = 0, the 
right side of (1) reduces identically to zero, irrespective of the relative magnitudes 
of r and s. A simpler expression for this integral has been derived by the present 
author in terms of complete elliptic integrals of the first and third kinds, by making 
use of the representation ([3], Chap. III, Art. 6): 





(5) rJo(rt)J,(st) = [  Sdt\/ BEF — Bre cnn 6) 60. 


Differentiation of (5) gives 











(6) xJo(rt)Ji(st) = [ Silt? + # — 2rs cos 4] 73 re oe o 


Thus 





r/ é Jo(rt)Jx(st) dt = [ e™ at f Silt? + 8 — 2rs cos ¢) 
0 0 ° 


(7) (s — r cos ¢) do 
, | VP + 8 — 2rs cos ¢ 

















-[ (s — r cos ¢) do m 
“LZ 


Pet 2 a dt. 
STS Wank é "J, [tr + 8 — 2rs cos 4] 


Carrying out the integration with respect to “t” gives 





_ * — (8 — rcos ¢) 
Pp = — 
(8) a a I, r+ s — 2rs cos ¢ 





SS 2S 2 ———— 
[1 ~ Vp + 2 + & — 2s cos A 
The first integral on the right side of (8) is known ((4], formula 860.2) : 


" [s—reos¢ldo _ 
fo tet = ais, ren 





(9) = 2/2s, r= 6 
= 0, r> 8, 


and the second can be converted to elliptic integrals in standard form by sub- 
stituting ¢ = « — 26 and by rearranging, thus, 














[ (s — r cos ¢) de as Bethel 

o (r° + 8° — 2rs cos ¢)V/p? + ° + #& — 2scosd 8r/p*? + (r + 8) 
{s - ff [2 de a at — - a ‘SR 1 

0) Vas) acvatovicrant! yo 





. 1 _ je~F ) 
~ sV/p + r+ (ot II (a, &) + K(k)}, 
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where 
bP ws 4rs 
(r+ sy’ 
ag i SE 
P+ (r + 8)’ 


and K(k) and [[(a, &) are complete elliptic integrals of the first and third kinds, 
respectively. An alternative form which expresses the result in terms of Heuman’s 
Lambda Function can be obtained from formula 413.01 of [2]: 














[ sss — TF C08 g) do 
(11) ae (7° + 8° — 2rs cos ¢)~/p? + r? + 8 — 2rs cos ¢ 
aot, ree 
s \2(s—4r) Vp + (r +.8)4)’ 
with 
p 


sin 8B = VP + - = 22° 
The final result i: 


Il 


a\— 





[Pe drt) di(st) at E ohh) — —oee |. ‘ia 


- aV/p+(r+s)? 


i 1 pK(k) | 
12) off oe es 
( 8 - > 2p ote 4r? ° 

if pK(k) | 

= — A ’ k — ——— ---———4 9 8, 

r E o(B, k) « VP + 4s) r>s 

where 
pute Sie “a ee 
eter "Gare 


Equation (12) whes * = s, incidentally, gives a particularly simple expression for 
Weeg’s integral in :hose cases where 7 = 1. 
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Integrals of Products of Laguerre Polynomials 
By R. D. Lord 


In a recent note in this journal Gillis and Weiss [1] have evaluated as a finite 
sum, and found a recurrence relation for, the integral 


(1) aa [ e*L,(x)L,(2)L,(x) dz, 
0 


which gives the coefficients in the expansion 

L(x) L(x) = > a L(x), 
where L,,(2) denotes the Laguerre polynomial, 
(2) La(e) =X (-1"(") s/t 


This attracted my interest since about ten years ago I considered the more general 
integral of the product of generalized Laguerre polynomials 


(3) [ ea LS” (a, x) LSS” (a, x) --- LS” (a, x) de, 
0 


with the object of finding cases which simplified. I tried methods which I had 
used with the corresponding problem for Hermite polynomials [2], but met with 
little success. But the special case (1) is easily dealt with as follows. 

The generating function for Laguerre polynomials is 


= a oe ru 
(4) dX (-1)'L,(2) “u = yo (7), 





and hence that for the product of three polynomials is 


7 > (—1)' L(x) D(a) L(x) uy" us" us’ 


0 s=0 t=—0 


Ms 


r 





- " =) uy U2 _ 
= [(1 + wm)(1 + we) (1 + us)) ep2(; + % + 1+ w% * i+ =): 


On multiplying by e* and integrating from 0 to ~, we get 


@ nm wo 
» +t t —1 
(5) } ® > (1) Cg ttty Ue*ug = (1 — gig — Ugty — Wyte — Zuyugus) 


r=0 s=0 t=—0 
as the generating function for C,,,. 'f we write the right hand side of (5) as 


7 (tog + Ug, + Ute + Zuytwou;)”, 


r+se+t 


we immediately see that C,,, is an integer with the sign (— ) . Further expan- 
sion of each term by the multinomial theorem gives Watson’s formula [3] 


Received 18 March 1960. 
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— (_o)rtett —2n n! . 
(6) — Lz (n — r)!(n — 8)!(n - I(r +8+ 4 — Qn)! 
where in the summation all factorials have non-negative arguments. 
If we transpose the right hand side of (5) and equate coefficients of u"u2'us', 
we have 


(7) Cres -_ ee + Cntieae + Cea o-1,2 + Be ct,ox 0685 


which is the recurrence formula with simplest coefficients, especially as the right 
hand side can be treated as the sum of five terms. Even though all suffixes vary, 
it probably provides the quickest way of computing all values of C,,, up to a given 
set of r, s, t. For some machines at least, it may well give the quickest way for 
calculating a given C,,; , and it provides an easy method for desk-machine computa- 
tion when r, s and ¢ are small. For computations by other methods it provides a 
simple check. Other checks may be obtained by giving w special values in (5) and 
equating coefficients of u2"u;‘. Putting wu. = —1, we get 


(8) y Crt = 1, 
given by Gillis and Weiss; putting u. = —4, we get 


(9) Er Cu = (° + 2, 


§ 
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The Evaluation of Certain Probability Integrals 


By Irwin Greenberg 


A problem which often arises in the field of probability and statistics is the 
following: 

Assume that there are n independent stochastic processes and that the kth 
process has an output distribution f,(2,). The probability that the jth process 
yields a higher output than any of the others is 


(1) P; = [then I [. tee dx, dx;. 
kAj 


In certain special cases, equation (1) is easily integrated; for example, if the 
various f;,(2,) are all uniform or exponential distributions. 
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In certain other cases, equation (1) can be put into a form which is readily 
evaluated by hand computation, using well known statistical reference tables. This 
paper will illustrate one of these cases; for n = 2 or 3 and the z, having a normal 
distribution. 


Under the assumption of normality, 


an 2 
(2) filta) = (Wize) texp[-3(#=—")]; k= 1,2,3 
k 
where », and o; are the mean and standard deviation, respectively. Letting 
(3) X,-=—, 
71 


equation (2) becomes 


2 
(4) fi( Xi) = (VJ 2x ox) ep| - 5 (#*") | 
where 
(5) py’ = eal - o;! = da 
71 1 


The probability that the output of process 1 exceeds the output of process 2 is 


n\2 
P, = (2me;') yee ‘exp| - 5 Xi - 3 (%5"*) ax. ax, 


= (Vx) ay exp (- 5) dt, 
where 


(7) Mz = —p'/V(or)? + 1. 


The simplification of equation (6) is obtained by expanding the exponent, 
grouping terms, completing the square in X, and reversing the order of integration 
after making the transformations: 


(8) Xi = XY’; Xe = X)’ + X,’ + be’. 


Equation (6) is the cumulative normal function and is tabulated in most texts 
on statistics. In a similar manner, it can be shown that the probability that the 
output of process 1 exceeds the output of process 2 and of process 3 is 


P, = [(24)*o0' 03/7 A ¥ [. “exp| - 3 ski - 5 (B52) 
nm (25) ] dX, dX, dX, 
2 o3' 
(9) M4 
= |2 1 vi-a} rs [. exp| - x0 — SU ay (te — Bele te 


+ “| dt, diz , 


(6) 
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where 
Mz = —m'/V(or’)? +1 
(10) Ms = —ws/V(o;')? + 1 
=V(6/? +1 V(o!/) +1. 


Equation (10) gives the volume under the bivariate normal probability surface 
with correlation coefficient r. These volumes are tabulated in [1]. 
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The Congruence 2”* = 1 (mod 7’) for p < 100,000 


By Sidney Kravitz 


Fréberg has previously announced [1] the computation of the Fermat remainders 
corresponding to all odd primes less than 50,000. His results show that p = 1093 
and p = 3511 are the only solutions of the congruence 2?" = 1 (mod 7’) in that 
range. 

The residues of 2”"' (mod p’) have been computed for 50,000 < p < 100,000 
on an IBM 650 system at Picatinny Arsenal. No residue congruent to 1 was found 
corresponding to a prime in this range. 

A copy of the table of residues has been deposited in the Unpublished Mathe- 
matical Tables file. 
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2. N.G. W. H. BEEcER, “On a new case of the congruence 2? ' = 1 (mod p*),’’ Messenger 
Math., v. 51, 1922, p. 149-150. 
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Values of = / (=) dt 
5 t 
By Kasaburé Harumi,{ Shigetoshi Katsura,t and John W. Wrench, Jr.t{ 


In discussing the equation of state for the molecules of one-dimensional square 
well potential [1], the first two authors required the numerical values of the integral 


if sint\”, a ne (") i 
I, = - : (==) dt = 2—“T(n) 2, (—1) (n 2p) . 


Inasmuch as these values seem to be of use in other applications, and apparently 
have not been previously tabulated, the first two authors calculated on SENAC-1 
(Sendai Automatic Computer 1) a six-place table of J, for n = 1(1) 30. 

The last author reviewed this table and recalculated the data, finding the 
corresponding exact rational values. 

The authors have decided to publish the ten-place table of J,, which is appended. 


TABLE oF 7, = - [ (=) dt 
rs t 











” | In n Tn 

1 1.00000 00000 16 | 0.34224 02614 
2 1.00000 00000 17 0.33220 82691 
3 0.75000 00000 | 18 0.32300 93942 
4 0.66666 66667 | 19 0.31453 44009 
5 0.59895 83333 | 20 0.30669 31017 
6 0.55000 00000 | 21 0.29941 02903 
7 0.51102 43056 22 0.29262 26872 
8 0.47936 50794 | 23 0.28627 66141 
9 0.45292 09682 24 0.28032 61985 
10 0.43041 77690 | 25 0.27473 19735 
ll 0.41096 26428 | 26 0.26945 97712 
12 0.39392 55652 | 27 0.26447 98425 
13 0.37884 40845 | 28 0.25976 61480 
14 0.36537 08695 29 0.25529 57845 
15 0.35323 91567 | 0.25104 85320 


{| Department of Applied Science 
Tohoku University 
Sendai, Japan; and 

tApplied Mathematics Laboratory 
David Taylor Model Basin 
Washington 7, District of Columbia 


1. S. Katsura & K. Harumi, “‘A note on the Born-Green linearized integral equation,” 
Phys. Soc. of London, Proc., v. 75, 1960, p. 826-832. 
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380 J. BOERSMA 


Computation of Fresnel Integrals 


By J. Boersma 


Two approximations, one valid for x less than 4 and the other valid for z larger 
than 4, have been established by means of the r-method of Lanczos [1] for the 
Fresnel integrals defined in the form 


f(z) = [ = dt = C(x) — éS(z). 


These approximations are the following: 








11 n 
(1) For Os2zs4 f(z) = ¢™ o> (a, + ib,) () 
n=0 
y 11 n 
(2) Forz =4 f(z) = vo +e” / Dd (en + tdn) (4) ° 
2 Zt n= Zz 
The numerical values of the coefficients a, , b, , c, and d, are given by 
a = +1.595769140 bo = —0.000000033 o& = O dy = +0.199471140 
a, = —0.000001702 bi = +4.255387524 c, = —0.024933975 d,; = +0.000000023 
a2 = —6.808568854 be = —0.000092810 Cc. = +0.000003936 d, = —0.009351341 
a; = —0.000576361 bs = —7.780020400 cs = +0.005770956 d; = +0.000023006 
a, = +6.920691902 b, = —0.009520895 | a= +0.000689892 | d, = +0.004851466 
a; = —0.016898657 bs = +5.075161298 | cs = —0.009497136 | d; = +0.001903218 
ag = —3.050485660 | bs = —0.138341947 | cs = +0.011948809 | d,s = —0.017122914 
a; = —0.075752419 by = —1.363729124 | c, = —0.006748873 | d; = +0.029064067 
as = +0.850663781 bs = —0.403349276 | cs = +0.000246420 | ds = —0.027928955 
a = —0.025639041 bp = +0.702222016 | co = +0.002102967 dy = +0.016497308 
aio = —0.150230960 bio = —0.216195929 Cio = —0.001217930 dio = —0.005598515 
a1, = +0.034404779 bu = +0.019547031 | Cu = +0.000233939 | dy, = +0.000838386 








The derivation of these approximations is given in [2]. 
The maximum error is 1.6 X 10~° for the first approximation and 0.5 x 10° 
for the second approximation. 


Mathematical Institute 


University of Groningen 
The Netherlands 
1. C. Lanczos, Applied Analysis, Prentice Hall, Englewood Cliffs, N. J., 1956. 


2. J. Boersma, ‘‘On a numerical method for the computation of Fresnel integrals’’, Report 
TW 2, Math. Inst., Univ. of Groningen, 1960. 


Received March 2, 1960. 
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REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 





69[A, J, L, M).—I. M. Rysuix & I. 8. Grapsre1n, Summen-Produkt- und Integral- 
Tafeln: Tables of Series, Products, and Integrals, VEB Deutscher Verlag der 
Wissenschaften, Berlin, 1957, xxiii + 438 p., 27 em. Price DM56. 


This volume of tables consists of a translation of the Russian third edition [1] 
into parallel German and English text. That edition has now been improved and 
augmented by the incorporation of corrections listed on a sheet of corrigenda ac- 
companying the third edition, the addition of supplementary remarks in the Ap- 
pendix, and the inclusion of an extensive supplementary bibliography, which 
consists of bocks and monographs pertaining to integral transforms, special func- 
tions, and indexes of mathematical tables. 

The preparation of the third edition was carried out by I. 8. Gradstein, follow- 
ing the death in the Second World War of I. M. Ryshik, who was responsible for 
the first two editions. A description of the contents of the first edition has appeared 
in a detailed review by R. C. Archibald [2]. 

The present book represents an extensive revision of the earlier editions. Major 
changes include the deletion of sections relating to the calculus of finite differences 
(including formulas for numerical quadrature), the addition of a completely new 
chapter on integral transforms, and the enlargement of the chapters on special 
functions. 

Information on each of the special functions—in particular, elliptic, cylindrical, 
and spherical—is presented systematically. Such information generally includes 
definitions; representation by integrals, series, and products; asymptotic formulas; 
functional equations; special values; and theorems relating to characteristic proper- 
ties. 

An introductory section entitled “On the Arrangement of the Formulae” ex- 
plains the arrangement of the contents of the chapters on elementary functions 
and on their definite integrals, and sets forth innovations in the arrangement of 
definite integrals, which in previous editions followed closely the classification es- 
tablished in the classical tables of Bierens de Haan [3]. 

The usefulness of this volume is enhanced by references and cross-references 
for the sources of most of the 5400 formulas presented. Formulas are numbered 
decimally within each chapter, and the chapter numbers are used for the integer 
part, as is customary. Furthermore, a key for the references to the literature cited 
on p. 434 is described in the Preface. It seems appropriate to note here that in both 
the Russian third edition and in this translation the list of numbered references 
consists of 40 items, although reference is made in several places in the book to a 
forty-first item and a forty-second that were apparently omitted inadvertently. 

The first two editions contained a table of 10D approximations to (2n — 1)!!/ 
(2n)!! and (2m — 1)!!/(2n)!\(2n + 1), form = 1(1)15, and to (2n — 1)!!/(2n + 2)!! 
and (2n — 1)!!/(2n + 2)!\(2n + 3), for n = 1(1)14. This numerical information 
is now supplemented by an original table of the Lobatschefsky function L(x) to 
7D for x = 0°(1°)10°, 6D for x = 11°(1°)30°, and 5D for x = 31°(1°)90°, com- 
puted by N. V. Tomantovaia under the supervision of B. L. Laptev. This function 
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is briefly discussed (p. 296-297) in the chapter on special functions. Additional 
numerical data also include exact values of the first 17 Bernoulli numbers and the 
first 10 Euler numbers, 10D approximations to ¢(n) for n = 2(1)11, and Euler’s 
constant and Catalan’s constant to 16D and 9D, respectively. I have examined 
all these data carefully, and the errors detected, together with errors in the for- 
mulas, are enumerated separately in this issue (MTE 293). 

Use of the book is facilitated by an elaborate index of special functions and nota- 
tions on p. 417-422. In addition to supplementary remarks and the bibliographies 
already mentioned, the Appendix contains (on p. 423-429) a discussion of the varia- 
tions in the notation and symbols used for special numbers and functions throughout 
the mathematical literature and a concise list of abbreviations (p. 432-433). 

The lucid expository style employed throughout is exemplified in the Introduc- 
tion. Here, a systematic summary of definitions and theorems relating to infinite 
products and infinite series of various types supplements the list of relevant formu- 
las. Similar explanatory text serves as introduction to several of the subsequent 
chapters and their subdivisions. 

Typographical errors found in the text are minor and do not. detract from the 
intelligibility of the textual material. The typography, especially in a compilation 
of such a large number of formulas, is uniformly excellent, and the appearance of 
the book is attractive. Professor Archibald’s opinion that the first edition was 
“undoubtedly of considerable value for any mathematician to have at hand” 
certainly holds true for this latest version. 

J.W.W. 


1. I. M. Ryzarx & I. S. GrapsHtein, Tabliisy Integralov, Summ, Riadov i Proizvedenit, 
[Tables of Integrals, Sums, Series and Products], The State Publishing House for Technical 
and Theoretical Literature, Moscow, 1951. 

2. R. C. Arcnrpatp, RMT 219, MT'AC, v. 1, 1943/45, p. 442. 

3. BreRENs pE Haan, Nouvelles Tables d’Intégrales Définies, Leyden 1867. Reprinted by 
G. E. Stechert & Co., New York, 1939. 


70(G].—_Evucene Prance, An Algorism for Factoring X” — 1 over a Finite Field, 
AFCRC-TN-59-775, U. S. Air Force, Bedford, Mass., October 1959, iii + 20 p., 
27 cm. 


An algorism is given for factoring X" — 1 over the finite field F, of g elements. 
This can be of use in constructing another finite field over F, , in constructing a 
linear recursion of period » over F,, or in constructing cyclic error-correcting 
group codes. The algorism has two parts: Step 1, the construction of the multi- 
plicative identities of the minimal ideals of F,[X]/[X” — 1]; Step 2, the use of these 
idempotents in the construction of the irreducible factors of X” — 1. 


AvuTHOR’s ABSTRACT 


71[G|.—M. Rotrensere, R. Brvins, N. Merropouis & J. K. Wooten, Jr., The 
3-j and 6-7 Symbols, The Technology Press, Massachusetts Institute of Tech- 
nology, Cambridge, Massachusetts, 1960, viii + 498 p., 29 em. Price $16.00. 


Wigner’s 3-j symbol is closely related to the Clebsch-Gordan coefficients used 
in the coupling of angular momenta. If J; and J. are coupled to give J, with j, j:, 
je as the total-angular-momentum quantum numbers and m, m , mz as the quantum 
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numbers for the z-components, the expansion coefficient giving the coupled states 
in terms of the uncoupled are 


(jxjojm | jrmijem.) = (— 1)*-**-"(2j + 1)” ts nl be 
Here the symbol on the left is the expansion coefficient in the notation of Condon 
and Shortley, Theory of Atomic Spectra; the last symbol on the right is the Wigner 
3-j symbol. The advantage of a tabulation of the 3-7 symbols, 


(4 ji a 
m ™ mM)’ 


rather than of expansion coefficients results from the high degree of symmetry of 
the 3-j symbols. At most a sign change results from an interchange of columns or 
from changing the signs of all the m’s. Thus, from these tables, which are restricted 
to 

paZzp2zj. and m sO, 


all expansion coefficients can be obtained for any j = 0, 3, 1, 3, --- , 8. 

The 6-j symbols occur in the coupling of three angular momenta J, , J: , and 
J; . One can either couple J, and J; first to obtain J’, and then couple J’ to J; to 
obtain J—this coupling scheme results in quantum numbers jj , je , j’, js ,j, m—or 
one can first couple J: and J; to obtain J”, and then couple J” to J, to obtain J— 
this scheme results in quantum numbers j2 , js, 7”, j1, j, m. The overlap integral 
between these two representations is given by 


where the last symbol is a 6-7 symbol. The 6-7 symbol tabulated, 

{jr je js 

\h le Is) . 
has sufficient symmetry that it need be listed only for 

pAZzhpZzijs; AZh, ke2h, R2kh. 
Within these restrictions, it is listed for all half-integral values of the j’s and l’s 
from 0 to 8. 
The tables were computed on the MANIAC II at the Los Alamos Scientific 

Laboratory. An adequate 40-page introduction describes the symbols and their 
uses. Since the symbols are the square-roots of rational fractions, the squares are 


tabulated as powers of primes in a shorthand notation, with an asterisk used to 
denote the negative square root. For example, the entry *1510,2221 is to be in- 


terpreted as 
_f Sw Poem i 
2! x 5' X 11? X 13? X 17? 
GEORGE SHORTLEY 


Booz-Allen Applied Research, Inc. 
Bethesda 14, Maryland 
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72{I|.—Hersert E. Sauzer & GENEvIEVE M. Krimsro, Tables for Bivariate Oscu- 
latory Interpolation Over a Cartesian Grid, Convair-Astronautics, Convair Divi- 
sion of General Dynamics Corporation, San Diego, California, 1958, 40 p. 


Formulas are developed for binary polynomials P(z, y) which agree together 
with the partial derivatives P.(z, y) and P,(z, y), with f = f(z, y), fe = f.(z, y) 
and f, = f,(z, y) at n specified points. They have the advantage over ordinary 
bivariate interpolation of packing 3n conditions into n points. Unlike univariate 
polynomial osculatory interpolation which always possesses a solution for any ir- 
regular configuration of fixed points, a binary polynomial of prescribed form may 
not satisfy those 3n conditions for any choice of interpolation points, or may fail 
for just certain special configurations. Explicit formulas or methods are developed 
for the general 2- to 5-point cases. For interpolation over any square Cartesian 
grid of length h, for suitable 2- to 5-point configurations of (x; , y;), according to 
the formula 


I(x, y) = f(xo + ph, yo + gh) ~ P(xo + ph, yo + gh) 
(1) Space ‘ ay 
= 2 {A (p, afi + MBI (p, fx; + CI, afull, 
tables of exact values of A; (p, q), Bi“ (p, q) and C;(p, q) are given for p and q 
each ranging from 0 to 1 at intervals of 0.1. A closed expression for the remainder 
in (I) has not been found. In its place, formulas are derived for the leading terms 
in the bivariate Taylor expansions for the remainders. These formulas should cut 
down the number of needed strips in the numerical solution of Cauchy’s problem 
for first order partial differential equations by the method of characteristic strips. 


AvuTHoR’s ABSTRACT 


73(K].—D. E. Barton & F. N. Davin, “A test for birth order effect,’”’ Ann. Human 
Genetics, v. 22, 1958, p. 250-257. 


In an ordered sequence of trials it is known that there were r; occurrences and 
2 non-oceurrences of a particular event. The question at issue has to do with the 
randomness of the positions of the occurrences in the sequence vs. a tendency to 
appear either at the ends or in the middle of the sequence. A test criterion is ob- 
tained by dividing the sequence between the Rth and the (R + 1)st event, where 
rT, + r2 = 2R or 2R + 1, and then assigning ranks 1, 2, --- to the events by order 
of position beginning at the point of division and proceeding to the left and then 
again starting with 1 to the right. The sum of the ranks of occurrences as assigned 
is the test criterion S. For r, + re = 4(1)16 with r,; = 72, on the null hypothesis 
of random position of occurrences, the exact distribution of S is tabulated for each 
pair (7; , r2) for re = 2. For r; + re > 16, it is stated that the distribution of S is 
sufficiently closely approximated by a normal distribution. 


C. C. Craic 
University of Michigan 
Ann Arbor, Michigan 
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74(K].—G. D. Bernpt, “Power functions of the gamma distribution,” Ann. 

Math. Stat., v. 29, 1958, p. 302-306. 

If z is a random variable from a gamma distribution with frequency function, 
fo = f(x; B, v) = | B’T(»)| “2” exp (—z/8); 6 > 0, » > O and zx = 0; the fre- 
quency function for éx with 6 > 1 is f; = f(éz; 58, v). To test the null hypothesis 
on the mean, Ho: » = 6v, against the alternate, Hi: « = 68r,5 > 1, one may use 


the statistic a(z) with the critical region defined by a = / je dz. Then the power 


of this test is 7; = i f, dx. The mean of a random sample from a universe whose 
a(z) 


frequency law is fo = f(z; 8, v) obeys a gamma distribution with parameters 6/n 
and nv. For a = .01, .05, .1, charts are given for reading x; for 1 < 6 S 4 and 
v = 3, 1(1)5, 7, 10(5)50. 

C. C. Craic 
University of Michigan 
Ann Arbor, Michigan 


75(K].—H. Cuernorr & L. E. Moses, Elementary Decision Theory, John Wiley & 
Sons, Inc., New York, 1959, xv + 364 p., 24 em. Price $7.50. 

This book is an elementary approach in the theory of statistics through the 
theory of the strategy of games, and as such is a refreshing change from the usual 
run of elementary statistics textbooks. The authors state that only an understand- 
ing of high school (U.S.) mathematics is required, which is possibly optimistic. 
However, it is fair to say that the mathematical content of the book is not ex- 
cessive, the exposition being mostly by example. 

Chapter I gives the principles of decision and an introduction to minimax. 
Chapter II, entitled Data Processing, turns out to be our old friends graphical 
representation and means and standard deviations. No mention is made of grouping 
corrections. There are 38 pages on probability and random variables, both continu- 
ous and discrete, followed by a brisk treatment of utility and descriptive statistics. 
This chapter (IV) will be rather difficult for the beginner. 

The authors have now reached a stage where they can, and do, begin to dis- 
cuss strategies. Chapter V, ‘Uncertainty due to Ignorance of the State of Nature,” 
gives simple Bayes strategies, minimax, and expected regret. (The reviewer liked 
the remark, “‘it is difficult to visualize four-dimensional space.’’) Further chapters 
cover further Bayes strategy and the application to problems which might arise 
in what is termed ‘‘Classical” statistics, in testing hypotheses, and in estimation. 
There is a series of appendices in which some of the statements in the main body 
of the text are proved. 

This is an interesting book and may prove useful to those who see the inter- 
pretation of numerical data as just one more decision to take; it is greatly to be 
doubted whether it is of general utility. In one way, however, it is unique. Fisher 
and Neyman in their several ways might be said to have contributed to statistical 
decision theory, but are not deemed worthy of reference. 

F. N. Davip 


Statistics Departmen: 
University College London 
Gower Street 

London WC 1, England 
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76{K].—R. Doornsos & H. J. Pris, “On slippage tests. I,’’ Indagationes Mathe- 
maticae, v. 20, 1958, p. 38-46 (Proc. Kon. Ned. Ak. van Wetensch., v. 61, Sec. A, 
1958, p. 38-46); “On slippage tests. II,” Ibid., p. 47-55; “On slippage tests. 
III,” Ibid., p. 438-447. 


The tables, which appear in part III, are related to two of the special cases 
included in this series of papers. In the first, from each of k Poisson distributions, 
with means py; , a random drawing Z; is taken (i = 1, 2, --- k). To test the null 
hypothesis that u; = u,7 = 1, --- , k, for which the table is prepared, against the 
alternate that one of the w,’s is greater than the others which have equal values, 
the authors propose the statistic, max Z;. For k = 2(1)10 and the sum of the k 
observations, n = 2(1)25, values of max Z; are given for which the significance 
levels are near 5% and 1%. In each case the actual significance levels are given 
to 3D. 

In the second case, each of k objects is ranked by each of m observers. The null 
hypothesis under test is that each of the m rankings is independently and randomly 
chosen from the set of permutations of the integers 1, 2, --- , k. As a test against 
the alternate that one of the objects has a higher probability of being ranked low 
while the others are ranked in random order, the proposed statistic is min S; where 
S; is the sum of ranks assigned the i-th object (¢ = 1, 2, --- , &). Critical values 
S. of min S; for significance levels near a = .05, .025, .01 are tabled for m = 3(1)9 
and k = 2(1)10. Again in each case true significance levels are shown to 3D. 


C. C. Craie 
University of Michigan 
Ann Arbor, Michigan 


77(K|.—F. G. Fosrmr, ‘Upper percentage points of the generalized beta distribu- 
tion. III,” Biometrika, v. 45, 1958, p. 492-503. 


Let max denote the greatest root of | »B — (A + »B)| = 0 where A and B 
are independent estimates, based on », and v2 degrees of freedom, of a parent dis- 
persion matrix of a four-dimensional multinormal distribution. Define 


T.(4; p,q) = Pr(@max = k) 


with p = 3(» — 3), q = 3(n — 3). Employing methods similar to those used in 
two preceding papers [1], [2] for the two and three-dimensional cases, the author 
tabulates 80 %, 85 %, 90 %, 95 %, and 99 % points of I,(4; p,q) to4D for», = 5(2)195 
and vg = 4(1)11. 

C. C. Craie 
University of Michigan 
Ann Arbor, Michigan 


1. F. G. Foster & D. H. Ress, ‘‘Upper percentage points of the generalized beta dis- 
tribution. I,”’ Biometrika, v. 44, 1957, p. 237-247. [MTAC, Rev. 165, v. 12, 1958, p. 302] 

2. F. G. Foster, ‘‘Upper percentage points of the generalized beta distribution. IT,”’ 
Biometrika, v. 44, 1957, p. 441-453. [MT'AC, Rev. 167, v. 12, 1958, p. 302.] 


78(K|.—W. Herz & H. Kurncer, ‘Untersuchungen zur Frage der Verteilung von 
Objekten auf Platze,” Metrika, v. 1, 1958, p. 3-20. 


For the classical distribution problem in which k indistinguishable objects are 
randomly distributed into n distinguishable cells (as in Maxwell-Boltzmann 
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statistics) the authors take the number, s, of occupied cells as a statistic to test 
the hypothesis of uniform probability over the cells. Let P(s |, k) be the prob- 
ability density for s. The correspondence is noted between this distribution and the 
results of a series of n drawings from a discrete distribution in which the random 
variable assumes only the values 0, 1, 2, --- , and in which the sample sum is k 
and the number of non-zero values is s. In developing a recursion formula for 
P(s|n, k) it is shown that the uniform distribution over cells arises from the Pois- 
son distribution, and the binomial and negative binomial distribution give par- 
ticular non-uniformities. The function tabulated is Z;.., which is defined under 
the hypothesis of uniformity by )°74 P(s|n, k) < a and >-7**** P(s|n, k) 
> a, for a = .05, .01, .001;n = 3(1)20, and ranges of k varying from (3, 15) for 
n = 3to (2, 100) for n = 20. 

C. C. Crate 
University of Michigan 
Ann Arbor, Michigan 


79(K].—A. Hurrson, ‘Further critical values for the sum of two variances,” Bio- 
metrika, v. 45, 1958, p. 279-282. 


Let s;, i = 1, 2, be an estimate of the variance «7 with f; degrees of freedom 
so that f,s;/o;’ is distributed as x’ with f; dif. To assign confidence limits to the 
form 4o;° + soz’, where \; and ), are arbitrary positive constants, the author has 
previously [1] tabulated upper and lower 5% and 1% critical values of 


(As8: + vse") /(ror® + Aoo:”). 


The present tables are an extension, giving upper and lower 23% and 4% 
critical values for the same function to 2D for \y8;°/(A18; + 282’) = 0(.1)1 and 
fi, fe = 16, 36, 144, ~. 

C. C. Craie 
University of Michigan 
Ann Arbor, Michigan 


1. A. Hurrson, ‘‘A method of assigning confidence limits to linear combinations of vari- 
ances,’’ Biometrika, v. 42, 1955, p. 471-479. [MTAC, Rev. 19, v. 12, 1958, p. 71.] 


80[K].—SoLtomon Kuuipack, Information Theory and Statistics, John Wiley & 
Sons, New York, 1959, xvii + 395 p., 24 em. Price $12.50. 


This interesting book, which discusses logarithmic measures of information 
and their applications to the testing of statistical hypotheses, contains three ex- 
tended tables in addition to a number of shorter or more specialized ones. Table I 
gives log, n and n log, n to 10D for n = 1(1)1000. Table II lists values of 


D. log. + (1 — py log, 2 to 7D for pr, ps = .01(.01).05(.05).95 
(.01).99. Table III gives 5% points for noncentral x* to 4D with 2n degrees of 
freedom for n = 1(1)7 and noncentrality parameter 6° for 8 = 0(.2)5. As it is 
stated, this is taken directly from an equivalent table of R. A. Fisher [1]. 


C. C. Craie 
University of Michigan 
Ann Arbor, Michigan 
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1..R. A. Fisuer, ‘‘The general sampling distribution of the multiple correlation,’’ Proc. 
Roy. Soc., A., 1928, p. 654-673. See p. 665. 


81(K].—G. J. LizBERMAN, “Tables for one-sided statistical tolerance limits,” 
Industrial Quality Control, v. 14, No. 10, 1958, p. 7-9. 


Given a sample of n from N(y, o°), it is desired to determine from the sample 
a quantity a (or b) such that with probability 7, the interval (— ©, a) (or the 
interval (b, ~)) will include at least the fraction 1 — a of the population. The 
tables give values of K to 3D forn = 3(1)25(5)50, y = .75, .9, .95, .99, and a = 
.25, .1,:.05, 01, .001, such that a = X — Ks andb = X — Ks, where X is the sample 
mean and S’ is the usual unbiased estimate of o”. For more extensive tables and a 
more complete discussion see [1]. 

C. C. Craic 

University of Michigan 
Ann Arbor, Michigan 


1. D. B. Owen, Tables of Factors for One-sided Tolerance Limits for a Normal Distribution, 
Office of Technical Services, Dept. of Commerce, Washington, D. C., 1958. [See RMT 82.] 


82(K].—D. B. Owen, Tables of Factors for One-sided Tolerance Limits for a Normal 

Distribution, Sandia Corporation, SCR-13, 1958, 131 p., 28 em. Obtainable from 

“the Office of Technical Services, Dept. of Commerce, Washington 25, D. C. 
Price $2.75. 


Given a sample of n from N(y, o), with Z the sample mean and S’ the usual 
unbiased estimate of o°, these tables give values of k for which 


Pr[Pr(z < @ + ks) = P| = +. 


As stated, Table I is a reproduction of one given by Johnson & Welch [1] in which 
values of k are given to 3D for y = .95, n = 5(1)10, 17, 37, 145, « 
and P = 0.7(.05).85, .875, .9, .935, .95, .96, .975, .99, .995, .996, .9975, .999, .9995. 
It is also explained that Table II was obtained from Resnikoff & Lieberman’s 
table of percentage points of the noncentral ¢-distribution [2] appropriately modi- 
fied to give k values to 3D for n = 3(1)25(5)50, © and P = .75, .85, .9, .935, 
.96, .975, .99, .996, .9975, .999 for y = .75, 9, 95. For y = .99, .995, 
n = .6(1)25(5)50, ©, while P has the same range as before. The more extensive 
Table III gives values to 5D obtained by an approximative method due to Wallis 
[3], for nm = 2(1)200(5)400(25)1000, ~, P = .7, 8, .9, .95, .99, .999, and y = 
-75 «8, .9, .95, .99, .999. For small m and the larger values of P and y, the approxi- 
mation breaks down and the entry is left blank or given with a warning sign that 
comparison should be made with neighboring values. (However it looks to the 
reviewer as if this sign has been omitted from the entries for n = 2, P = .99, .999, 
and y = .999.) Finally Table IV is obtained from Bowker’s table of two-sided 
tolerance limits [3] by an approximate procedure suggested by McClung [4] to 
give conservative values of k for one-sided limits. Here values are given to 3D for 
nm. = 2(1)102(2)180(5)300(10)400(25)750(50)1000, ©, P = .875, .95, .975, .995, 
.9995, and y = .75, .9, .99. 

In an appendix auxiliary tables compare values in the four tables for selected 
values of the four parameters. The maximum difference shown between Tables I 
and IT is .01. It is concluded that values in Table III will probably be underesti- 





te 


a 
Cc 
8 
a 
r 
c 





REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 389 


mates for y S .95 and overestimates for y 2 .99, while in Table IV, k is probably 
underestimated for P = .875 and overestimated for the other P values. Differences 
shown between Table II and Table III values in a few cases exceed 20% of the 
presumably more accurate Table II values and differences shown between Table 
II and Table IV sometimes exceed 10% of the Table IT values. 
C. C. Craice 

University of Michigan 
Ann Arbor, Michigan 

1. N. L. Jounson & B. L. Wetcu, “Applications of the non-central ¢-distribution,”’ 
Biometrika, v. 31, 1939, p. 362-389. 


2. G. J. Resnrxorr & G. J. Lizserman, Tables of the Noncentral t-Distribution, Stanford 
University Press, Stanford, Calif., 1957. 


3. C. Ersennart, M. W. Hastay & W. A. Wauuis, Techniques of Statistical Analysis, 
McGraw-Hill Book Co., New York, 1947. 


4. R. M. McCuuna, ‘First aid for pet projects injured in the lab or on the range or what 
to do until the statistician comes,’’ U. 8. Naval Ordnance Test Station Technical Memoran- 
dum No. 1113, October 1955. 


83[K].—K. V. RamMacHAnpRAN, “On the Studentized smallest chi-square,’ Amer. 
Stat. Assn., Jn., v. 53, 1958, p. 868-872. 


Consider the F statistics, . . -- t = 1,2,--- ,k,in which S,, S&, --- , S, and 
S are mutually independent, with each S;;.2 having a x’ distribution under the 
null hypothesis with ¢ degrees of freedom and S/o’ a x’ distribution with m d.f. 


There are numerous applications of statistical methods, a few of which are dis- 
cussed, in which one needs the value of V for which Pr | Smin m =Vi\|=1-a. 


S t 
Sue . . for values of t, m and k as follows: 





The author tabulates lower 5 % points of 


Fort = 1,m 2 5,k = 1(1)8 to 18; fort = 2,5 < m < 10 and m 2 12, k = 1(1)8 
to 3D; for ¢ = 3, 4, 6, m = 5, 6(2)12, 20, 24, ~, k = 1(1)8 to 3D; 
for ¢ = 1(1)4(2)12, 16, 20,m = «~,k = 1(1)8 to 3D; fort = 1(1)4(2)12, 16, 20, 
m = 5, 6(2)12, 20, 24, ©, k = 1, 2,3 to 3D. 

C. C. Craie 
University of Michigan 
Ann Arbor, Michigan 





84(K].—A. E. Sarwan & B. G. Greensere, “Estimation of location and scale 


parameters by order statistics from singly and doubly censored samples. Part 
II.,” Ann. Math. Stat., v. 29, 1958, p. 79-105. 


This paper, a continuation of a previous one [1], is mainly devoted to an ex- 
tension of tables given in the earlier paper to cover samples 11 < n < 15 and to 
a discussion of efficiencies of the estimators used. Samples of n are from N(x, o°); 
r, and rz observations are censored in the left and right tails respectively (rir2 2 0); 
and ~ and o are estimated by the most efficient linear forms in the ordered un- 
censored observations. Table I gives the coefficients for these best linear systematic 
statistics to 4D for all combinations of r; , r2 for nm = 11(1)15. Table II gives vari- 
ances and the covariance of these estimates to 4D for n = 11(1)15 and all pairs of 
T; , T2 values. In Table III efficiencies of the two estimates relative to that for un- 
censored samples are given to 4D for the same range of values of n and 7, , r2 . For 
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n = 12 and 15, variances and efficiencies relative to best linear systematic estimates 
are given for alternate estimates proposed by Gupta [2] for n > 10, and generalized 
in [1] to doubly censored samples, are given to 8D and 4D respectively for all r; , r2 . 
The authors state that extensions of Tables I, II, III to 8D for 16 < n < 20 are 
available upon application. 
C. C. Crate 

University of Michigan 
Ann Arbor, Michigan 

1. A. E. Sarnan & B. G. Greensere, “Estimation of location and scale parameters by 
order statistics from singly and doubly censored samples. Part I. The normal distribution 


up to 5, 2 of size 10,’’ Ann. Math. Stat., v. 27, 1957, p. 427-451. [MT AC, Review 141, v. 12, 
1958 


“rr K. Gupta, ‘‘Estimation of the mean and standard deviation of a normal population 
hoe a censored sample. ”” Biometrika, v. 39, 1952, p. 88-95. 


85(K].—J. M. Sencupra & Nikuitess Buarracuarys, “Tables of random normal 
deviates,”’ Sankhya, v. 20, 1958, p. 250-286. 


As explained by the editor in a foreword, this is a reissue of an original table of 
random normal deviates which appeared in 1934 in Sankhya [1]. Since errors had 
been discovered in the earlier tables, the new set was reconstructed by conversion 
of Tippett’s random numbers [2] to random normal deviates, as was the case before. 
After the present table was prepared, in 1952, as stated by the editor, it was learned 
that an identical table had been constructed in 1954 at the University of California. 
On comparison it was found that the two tables checked perfectly. As discussed 
in the text, rather extensive tests of the hypothesis that the entries were random 
drawings from N(0, 1) were applied with satisfactory results. These tables contain 
10,400 3D numbers. 

C. C. Crate 
University of Michigan 
Ann Arbor, Michigan 

1. P. C. Manaranosis, S. S. Boss, P. R. Roy & S. K. i “Tables of random 

ie ig — a normal cod Bhg ws Sankhya, v. 1, 1934, p. 289-328. 


C. Tipretr, Random Sampling Numbers, Tracts for Computers, No. XV, Cam- 
bridge University Press, London, 1927. 


86[(K].—Minorv Srorani, “Note on the utilization of the generalized Student 
ratio in analysis of variance or dispersion,” Ann. Inst. Stat. Math., v. 9, 1958, 
p. 157-171. 


In samples from a p-dimensional normal universe an important statistic, ap- 
plications of which are discussed in this paper, is Ty) = m tr LV in which L and 
V are two independent unbiased estimates of the population variance matrix with 
n and m degrees of freedom respectively. Tables are given for the 5% and 1% 
points of the distribution of 7,’ to 2D for m = 1(1)10(2)20 and 


= 10(2)30(5)50, 60, 80, 100. 
C. C. Crate 
University of Michigan 
Ann Arbor, Michigan 
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87(K].—Minoru Siorant & Masaru OzAwa, “Tables for testing the homogeneity 
of k independent binomial experiments on a certain event based on the range,” 
Ann. Inst. Stat. Math., v. 10, 1958, p. 47-63. 


Let k series of N trials each of a certain event be performed with the outcome 


of v; occurrences in the i-th series in which the fixed probability of occurrence was 
Pi,t = 1, 2,---, k. To test the null hypothesis of homogeneity: 


Pi = Pp = --* = P= P, 
Siotani had previously proposed the statistic, Ri(N, p), the range of the »; [1]. 
The tables in this paper give for N = 10(1)20, 22, 25, 27, 30; k = 2(1)15; 
p = .1(.1).5; 
a = .001, .005, .01(.01).06, .08, .1, the greatest r, for which 
Pr{Ri(N, p) = ri} < a@ + .0005. 
The cases in which for the r;, given, a < Pr{R.(N, p) = mn} < a + 0005 or 
a — .005 < Pr{R.(N, p) = nj <a 
are indicated by attaching a + or a — respectively to the value of 7; . 


C. C. Crate 
University of Michigan 
Ann Arbor, Michigan 


1. Mrnorvu Srorant, “‘Order statistics for discrete case with a numerical application to 
the binomial distribution,’’ Ann. Inst. Stat. Math., v. 8, 1956, p. 95-104. 


88[K].—P. N. Somervi..e, ‘““Tables for obtaining non-parametric tolerance limits,” 
Ann. Math. Stat., v. 29, 1958, p. 599-601. 


Let P be the fraction of a population having a continuous but unknown dis- 
tribution function that lies between the r-th smallest and the s-th largest values 
in a random sample of n drawn from that population. Then for any r, s = 0 such 
that r + s = m, Table I gives the largest value of m such that with confidence 
coefficient = y we may assert that 100P% of the population lies in the interval 
(r, s) for y = .5, .75, .9, 95, .99 and n = 50(5)100(10)150, 170, 200( 100) 1000. 
Table II gives 7 to 2D for the assertion that 100P % of the population lies within 
the range, (r, s = 1), in a sample of n for P = .5, .75, .9, .95, .99 and 


n = 3(1)20, 25, 30(10)100. 
C. C. Craic 


University of Michigan 
Ann Arbor, Michigan 


89/K|.—G. P. Srecx, “A table for computing trivariate normal probabilities,” 
Ann. Math. Stat., v. 29, 1958, p. 780-800. 


Let X, Y, Z be standardized random variables obeying a trivariate normal dis- 
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tribution law. The author finds Pr (X < h, Y s k, Z S m) in terms of three func- 
tions: 


G(x) = (2")*” exp (- ) dz, T (h, a) 


= (2x) [ fexp {— A°(1 + 2°)/2})(1 + 2”) az, 


and 
h 
S(h,a,b) = / T (as, b)G’(s) ds. 


The T-function has been tabulated by D. B. Owen [1, 2] and a table of S(m, a, b) 
is given in the present paper to 7D for a = 0(.1)2(.2)5(.5)8, b = .1(.1)1 and a 
range of values of m decreasing from 0(.1)1.5, © for a = 0(.1)1.2 to 0(.1).3, « 
for a = 6(.5)8. The tabulated values are believed accurate to 0.6 in the seventh 
decimal place. There is considerable discussion of the main problem, of properties 
of and relations among the functions used, and a numerical example is worked out. 
The method of construction of the table is given and the efficacy of linear inter- 
polation in it is discussed. 


C. C. Craic 
University of Michigan 
An Arbor, Michigan 


1. D. B. Owen, The Bivariate Normal Probability Distribution, Office of Technical Serv- 
ices, apiece of Commerce, Washington, D. C., 1957, [MTAC Review 134, v. 12, 1958, p. 
285-286. 


2. D. B. Owsn, ‘‘Tables for computing bivariate normal probabilities,’’ Ann. Math. Stat., 
v. 27, 1956, p. 1075-1090. [MT'AC Review 135, v. 12, 1958, p. 286.] 





90[K].—G. Tacuti, “Tables of tolerance coefficients for normal populations,’ 
Union of Japanese Scientists and Engineers, Reports of Statistical A pplication 
Research, v. 5, 1958, p. 73-118. 


The tolerance limits T,, T2 are to be determined so that with probability 
1 — a@ the interval (7; , T:) includes a given fraction, P, of the population. Follow- 
ing the method of Wald & Wolfowitz [1] for a sample from N(u, 0°), T; and T: are 
found by T; = ~ — kv/S,/v and T; = ~ + kv/S,/v, in which @ is an unbiased 
estimate of » with variance o'/n and S, is an independent error sum of squares 
with v degrees of freedom. As illustrated by the author this permits useful applica- 
tions in which n is not simply the sample size and v = n — 1 as is the case for the 
tables of Bowker [2]. The present tables give k to 3S for P = .9, .95, 99,1 -— a= 
9, .95, 99, n = .5(.5)2(1)10(2)20(5)30(10)60(20)100, 200, 500, 1000, © and 
vy = 1(1)20(2)30(5)100(100)1000, ©. The calculations were done with a slide 
rule and the author fears there may be errors up to one per cent. Some cursory 
comparisons with Bowker’s tables for » = n — 1 showed frequent differences in 
the third significant figure. 

C. C. Crate 

University of Michigan 
Ann Arbor, Michigan 
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1. A. Wap & J. Wo.Fowr17, ‘‘Tolerance limits for a norma! distribution,” 4nn. Math. 
a v. 17, 1946, 


2. AR nD, mae wall M. W. Hastay & W. A. Wa1is, Techniques of Statistical Analy- 
sis, McGraw-Hill Book Co., New York. 1947. (See p. 102-107.) 


91[K].—R. F. Tate & R. L. Goan, “Minimum variance unbiased estimation for 
the truncated Poisson distribution,” Ann. Math. Stat., v. 29, 1958, p. 755-765. 


For a sample of n from a population with the density function, e~r*/(1 — ¢~), 
z = 1,2, --- ,i.e., a Poisson distribution truncated on the left at z = 1, the authors 
derive the minimum variance unbiased estimation of 


Pee -S)-! 
Azo(Z) -_ n ( S — n C(n, t), 
in which ¢ is the sample sum and ©,” is a Stirling number of the second kind. 
Using an unpublished table of F. L. Micksa [1] of ©," for n = 1(1)t,t = 1(1)50, 
this paper contains a table of C(n, t) to 5D for n = 2(1)t — 1, t = 3(1)50. 





C. C. Craic 
University of Michigan 
Ann Arbor, Michigan 


1. Francis L, Mrxsa, Stirling numbers of the second kind, RMT 85, MTAC v. 9, 1955, p. 198. 


92(K, P].—P. A. P. Moran, The Theory of Storage, John Wiley & Sons, Inc., New 
York, 1960, 111 p., 19 cm. Price $2.50. 


This is a book about dams. Prof. Moran is at the Australian National University 
at Canberra, and I imagine that dams have great practical interest there. For 
many years he has been interested in estimating the probability that a dam will 
go dry or that it will overflow. He is also interested in how one finds a program of 
releasing water from a dam in such a way as to optimize the operations of a hydro- 
electric plant. 

The first chapter contains some basic information about statistics and probabil- 
ity. To spare 14 pages for this from a total of a mere 96 shows how necessary Prof. 
Moran considered it to be. 

The second chapter considers various general inventory and queueing problems 
analogous to dam problems. 

In the third chapter the author plunges into his favorite topic, dams. First 
he considers discrete time—he looks at his water level only oncea day. Under certain 
conditions distributions for the amount of water can be found, but two troublesome 
conditions occur which limit the regions of analyticity of the distributions. One is 
overflow. The other is running dry. If one ignores either or both of these, then he 
is dealing with an imaginary “infinite dam’’. Some queueing is analogous to an 
infinite dam, since there is no law limiting the lengths of queues. 

Another chapter is devoted to dams which have as input a continuous flow, 
and from which the release is continuous. 

In practice the inputs do not satisfy the assumption of independence, dry weeks 
tend to come in succession, so the results of the first four chapters are of limited 
applicability. Monte Carlo methods get estimates of the probabilities without 
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these restrictions, and in addition can be applied to configurations of dams com- 
pletely beyond other methods of analysis. Of course Monte Carlo has disadvantages 
of its own. An example is given of a complex configuration for which probabilities 
were urgently wanted. A large retaining wall of earth was to be built. Overflow 
would ruin it, so a diversion tunnel was to be built large enough to insure against 
this contingency. During the building the tunnel is closed to permit pouring the 
concrete at its mouth. If water accumulates too high behind the wall there will be 
danger of overflow, ruining the wall. This can be prevented by opening the tunnel, 
ruining its outworks but preferable to damaging the main wall. The critical height 
changes each day as the wall is built up. What are the chances of this decision 
being forced? 

The last few pages are devoted to ways of finding an optimum strategy for 
operating a hydroelectric system, or other program of releasing, replenishing, or 
otherwise tending the locks. The recommended solution is a method of successive 
approximations, which would probably be feasible only on a digital computer. 
The author suggests that a special analog device would be in order for the more 
complicated configurations. 

The analogies between dams and queues or inventories are not pursued beyond 
the third chapter, in which it is merely mentioned. If these analogies are indeed 
valid they deserve more treatment. Without this treatment the title is misleading, 
for we find we are storing only water. 


H. H. CAMPAIGNE 
National Security Agency 
Washington 25, District of Columbia 


93([M, X].—Jaxos Horn & Hans Wirticu, Gewdéhnliche Differentialgleichungen, 
Walter de Gruyter & Co., Berlin, 1960, 275 p., 24 em. Price DM 32. 


This book is the sixth completely revised edition of Jakob Horn’s Gewéhnliche 
Differentialgleichungen, which was published first in 1905. Like the previous editions, 
this book is intended for mathematicians, physicists, and engineers. In the selec- 
tion of the material somewhat greater emphasis has been given to subjects that 
lend themselves to applications. Nevertheless, this book is primarily an introduc- 
tion to the theory of ordinary differential equations. The text contains existence 
proofs and a comparatively detailed presentation of differential equations in the 
complex domain. 

Considerable space is devoted to special functions which arise from differential 
equations. Numerical and graphical methods of solution are treated in a brief 
chapter. Besides a thorough knowledge of differential and integral calculus on the 
part of the reader, a familiarity with the basic concepts of the theory of functions 
of a complex variable is assumed. 

No problem sections appear in the book, but numerous illustrative examples 
are provided. 

F. THEILHEIMER 
Applied Mathematics Laboratory 
David Taylor Model Basin 
Washington, District of Columbia 
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94 (Q, S).—Garrerr Birxuorr anv R. E. Lanoer, Editors, Proceedings of Sym- 
posia in Applied Mathematics, Vol. IX, “Orbit Theory,” (Proceedings of the 
Ninth Symposium In Applied Mathematics of the American Mathematical So- 
ciety, held at New York University April 4—-6, 1957, cosponsored by The Office 
of Ordnance Research, Ordnance Corps, U. 8. Army) American Mathematical 
Society, Providence, R. I., 1959, v + 195 p., 26 em. Price $7.20. 


The purpose of the book is, paraphrasing the words of the editors, to direct the 
attention of mathematicians to recent advances in celestial mechanics and, more 
importantly, to inform them of the problems that remain to be solved. Celestial 
mechanics owes its present form very largely to analysis, as it was developed in 
the eighteenth and nineteenth centuries. Whether modern mathematics can con- 
tribute anything important to the subject is a question that has hardly been ex- 
plored, and it is high time that it should be. 

Of the ten contributions by as many authors the first three deal with the mo- 
tions of particles in magnetic fields, the remaining seven with motions of particles 
in gravitational fields. The magnetic fields considered are those in particle accelera- 
tors, in the galaxy, and about a laboratory model of the earth. The gravitational 
fields are principally those of the earth and of the solar system, although one paper 
deals generally with the field about any massive particle, and one with a general 
planetary system. 

The various contributions are very uneven, ranging from rather trivial special 
applications of general formulae, through adaptations and modifications that are 
not trivial, to some important original contributions, both general and particular. 
Some authors describe what they have done themselves, some what others have 
done, and some what has not been done. Brouwer, Courant, and Olbert give special 
attention to unsolved problems; the references will be valuable to a mathematician 
not previously acquainted with their subjects. Herget and Eckert deal with prac- 
tical computation. 


G. M. CLEMENCE 
U. 8. Naval Observatory 


Washington, D. C. 


95[X].—Rvupours E. Lancer, Editor, Boundary Problems in Differential Equations, 
Proceedings of a symposium conducted by the Mathematics Research Center, 
University of Wisconsin, Madison, Wisconsin, The University of Wisconsin 
Press, Madison, 1960, x + 324 p., 24 cm. Price $4.00. 


This volume contains the nineteen papers presented at the symposium on 
‘Boundary Problems in Differential Equations” conducted by the Mathematics 
Research Center at Madison, Wisconsin during the period April 20-22, 1959. The 
papers are quite varied in nature and subject matter, as is clear from the table of 
contents given below: 

Boundary Problems of Linear Differential Equations Independent of Type 

K. O. Friedrichs, Institute of Mathematical Sciences, New York University 

Numerical Estimates of Contraction and Drag Coefficients 

Paul R. Garabedian, Stanford University 
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Complete Systems of Solutions for a Class of Singular Elliptic Partial Differential 
Equations 
Peter Henrici, University of California 
Application of the Theory of Monotonic Operators to Boundary Value Problems 
Lothar Collatz, University of Hamburg, Germany 
Upper and Lower Bounds for Quadratic Integrals and, at a Point, for Solutions 
of Linear Boundary Value Problems 
J. B. Diaz, Institute for Fluid Dynamics and Applied Mathematics, Univer- 
sity of Maryland 
Error Estimates for Boundary Value Problems Using Fixed-Point Theorems 
Johann Schroder, University of Hamburg, Germany 
On a Unified Theory of Boundary Value Problems for Elliptic-Parabolic Equa- 
tions of Second Order 
Gaetano Fichera, The Mathematical Institute, University of Rome, Italy 
Factorization and Normalized Iterative Methods 
Richard 8. Varga, Westinghouse Electric Corporation, Bettis Atomic Power 
Division, Pittsburgh 
Some Numerical Studies of Iterative Methods for Solving Elliptic Difference 
Equations 
David Young and Louis Ehrlich, The University of Texas. 
Presented by David Young 
Albedo Functions for Elliptic Equations 
Garrett Birkhoff, Harvard University 
A Numerical Method for Analytic Continuation 
Jim Douglas, Jr., The Rice Institute, Texas 
Stress Distribution in an Infinite Elastic Sheet with a Doubly-Periodic Set of 
Equal Holes 
W. T. Koiter, Technical University, Delft, Holland 
Some Stress Singularities and Their Computation by Means of Integral Equa- 
tions 
Hans F. Bueckner, Mathematics Research Center, U. S. Army 
Boundary Value Problems in Thermoelasticity 
Ian N. Sneddon, The Univerity, Glasgow, Scotland 
Some Numerical Experiments with Eigenvalue Problems in Ordinary Differen- 
tia Equations 
Leslie Fox, University Computing Laboratory, Oxford, England 
Dynamic Programming, Invariant Imbedding, and Two-Point Boundary 
Value Problems 
Richard Bellman, The Rand Corporation, California 
Remarks about the Rayleigh-Ritz Method 
Richard Courant, Institute of Mathematical Sciences, New York University 
Free Oscillations of a Fluid in a Container 
B. Andreas Troesch, Space Technology Laboratories, Inc., California 
A Variational Method for Computing the Echo Area of a Lamina 
Calvin H. Wilcox, Mathematics Research Center, U.S. Army, and California 
Institute of Technology 
Workers in Numerical Analysis will be particularly interested in the papers 
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of Friedrichs, Garabedian, Collatz, Varga, Schroder, Young and Ehrlich, Douglas, 
and Fox. The first author has a very short paper in which he gives an interesting 
outline of a unified approach to the numerical treatment of linear partial differ- 
ential equations irrespective of their type. The unified approach is said to also 
cover certain equations of mixed type. Unfortunately, the author did not have 
space to completely describe the conditions he must impose on the equations he 
treats. 

Garabedian describes a method that has been used to calculate axially sym- 
metric flows with free streamlines. In particular he discusses methods for calculating 
the contraction coefficient in the vena contracta. The method involves generalizing 
the differential equation governing the flow by introducing a parameter \ and 
studying the dependence of the solution as a function of X. 

Collatz’s paper is an expository one in which he discusses various definitions of 
monotonic operators and applies such definitions to the determination of bounds 
on the solutions of various problems. 

Schroder uses monotonic operators which satisfy a fixed point theorem to prove 
the existence of solutions to problems involving differential equations and boundary 
conditions. He also determines approximate solutions and error bounds by solving 
a so-called comparison problem. 

Varga discusses a class of iterative methods for solving a system of linear equa- 
tions which depend on the direct solution of matrix equations of matrices more 
general than tridiagonal matrices. He shows how such matrix equations can be 
directly and efficiently solved and, in addition, applies standard methods for ac- 
celerating convergence. 

Young and Ehrlich report on numerical experiments which attempted to deter- 
mine the extent to which theoretical results on the rate of convergence of the suc- 
cessive over-relaxation method for solving linear equations and for the Peaceman- 
Rachford method would apply for non-rectangular regions. The theoretical results 
are known for the latter method only in the rectangular case. In nearly every case 
it was found that the number of iterations using the Peaceman-Rachford method 
was less than was required using the successive over-relaxation method. However, 
approximately three times as much computer time is required for a double sweep 
of the former method as is required for a single step of the latter method. 

Douglas discusses the determination of an approximation of an analytic function 
of a complex variable inside the disk 0 < | z| < 1 when bounds on the function 
and its first two derivatives are known and when approximate values of the func- 
tions are known at p points equidistributed on the circle | z| = 1. An estimate of 
the error of the approximation is also obtained. 

Fox discusses a method for the determination of approximate proper values 
and proper solutions to single or systems of ordinary differential equations for which 
a reasonable approximation is already known for the proper value. The method 
involves the introduction of parameters such as initial values of the solution at 
one of the boundary points and the determination of improved values for these 
parameters by the Newton process. The method described is not new, but the 
applications made by Fox to fairly difficult problems give an impressive demon- 
stration of its power. 

Space limitations prevent the reviewing of the remaining papers in this volume. 
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They are of high quality. The organizers of the conference are to be congratulated 
on the papers solicited. The University of Wisconsin Press has produced a handsome 
volume by a photographic process which makes a very readable page. The rela- 
tively low cost of the volume is especially noteworthy. 

A. H. T. 





, “Operators for solution of discrete Dirichlet and 
Plateau problems over a regular triangular grid,’ May 1959, 29 em., 191 p. 
Deposited in UMT File. 


These tables list to 10D coefficients of a matrix operator for conversion of 
boundary values over an equilateral triangle to a discrete harmonic function over 
a regular triangular grid of 190 points in this triangle [1]. Sixty-three boundary 
values are involved, of which the three at the vertices do not influence the interior 
values of the function. The tables are useful in the approximate numerical solution 
of the Laplace equation over this triangular region. 

Solutions for smaller triangles have been placed in the UMT File by the same 
author [2]. 

Also included are tables giving 10D coefficients of the analog of the Douglas 
functional over this same grid. Specifically, these are coefficients of a quadratic 
form (using scalar multiplication) of vector functions from the grid points of the 
bounding equilateral triangle to some euclidean space such that the value of the 
form is the Dirichlet integral 


=5[€+@d 


where E and G are coefficients of the first fundamental form of the surface got by 
linear interpolation of the discrete harmonic vectors resulting from application of 
the operator described above to the boundary values. This is a discrete analog 
of the functional used by J. Douglas [3] in his solution of the Problem of Plateau; 
it has application in the approximate numerical solution of that problem. 


C. B. Tompkins 
University of California 
Los — California 


L. V. Kanrorovicu & K. I. Kryxov (translated by Curtis D. BensteEr) Aagegrinnsts 
Methods aa Hi [Biger Analysis, Nordhoff, Gronigen, Interscience, New York, 1958, p. 187-188. 
ILSON, JR., “Tables of inverses to Laplacian operators over triangular grids,”’ 
umt File MMTAC. No. 58, v. XI, 1957, p. 108. 
¥ Dove.as, “Solution of the pro lem of Plateau,’’ Amer. Math. Soc. Trans., v. 33, 
1931" p. 263-321. 





97(Z]. , Mathematical Programming and Electrical Net- 
works, John Wiley & Sons, Inc., New York, 1959, vi + 186 p., 24 em. Price 
$4.50. 


As the title indicates, the purpose of this little monograph is to explore the 
relationships of general programming problems and corresponding electrical net- 
works, with a view towards gaining physical insight and developing computational 
algorithms. The contents of the book essentially comprise the author’s doctoral 
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aissertation in the department of Electrical Engineering at M.I.T. The pages are 
offset reproductions of typescript. In a foreword by J. A. Stratton, it is stated 
that “there has long been a need for publication of research studies larger than a 
journal article but less ambitious than a finished book,” and with this apology the 
present volume is put forth. 

After an introductory Chapter 1, the general programming problem is presented 
in Chapter 2, along with discussions of convexity and concavity, the generalized 
method of Lagrangian multipliers due to Kuhn and Tucker, and duality. Chapter 
3 consists of basic material on electrical networks containing resistors, diodes, 
ideal transformers, voltage sources, and current sources. The electrical network 
problem, which is a set of linear equations with side conditions in the form of in- 
equalities due to the presence of diodes, is stated and shown to be equivalent to a 
quadratic programming problem (and its dual), viz., to find a feasible current dis- 
tribution which minimizes the power absorbed by the voltage sources plus one- 
half the power absorbed by the resistors, with a corresponding statement for the 
dual. The concept of the two-terminal, or terminal-pair, network is introduced in 
this chapter, with a discussion of the set of solutions (¢, 5), where « is the voltage 
between the terminals when current 6 enters one terminal and leaves the other. 
The set of all (e, 5) forms a “break-point curve”’ in the «5 plane, i.e., a non-decreasing 
polygonal graph. 

Chapter 4 is devoted to the problem of flow in a network, which includes alloca- 
tion, distribution, and assignment problems. It is shown that every flow problem 
can be realized by an electrical network containing only diodes, voltage sources, 
and current sources. Existence conditions from the theory of programming and 
non-redundancy assumptions are stated here in electrical network terms. Two 
algorithms are presented for the solution of diode-source networks, one correspond- 
ing to the primal and the other to the dual problem. They are similar to but more 
general than the procedure of Ford and Fulkerson for the transportation problem. 

Chapters 5 and 6 treat the general linear and quadratic programming problems. 
A procedure is described for “tracing” the break-point curve corresponding to a 
pair of terminals, which is very similar to what goes on in the simplex method of 
Dantzig. This procedure is used with electrical models of the general quadratic 
(including linear) problems, and two types of algorithms for their solution are 
described. Chapter 7 contains a brief and incomplete discussion of the general 
programming problem. An algorithm is proposed which is based on the method of 
steepest descents. The main body of text is followed by eight appendices in which 
proofs are given of the main theorems of programming. 

There are several typographical errors of a minor nature. The number of state- 
ments that are only partially true seems to be about par for a book on a mathe- 
matical subject written by an engineer. 


J. Bram 
Bureau of Ships 
Navy Department 
Washington 25, District of Columbia 











TABLE ERRATA 


286.—Car. Burrav, Tafeln der Funktionen Cosinus und Sinus, etc., Reimer, Berlin, 
1907. 


This table is one of the relatively few that give e* at .001 interval of z for values 
of z in the neighborhood of 9. (The range covered is 8.000 (.001) 9.809.) The values 
of e* are not all in one place; some are on pages 8-10, and the remainder are on 
page 44. 

A complete check of the 7-figure values of e* for z from 8.000 to 9.210 revealed 
65 last-figure unit errors. 


CHARLES R. SExTON 


2947 Elmwood Court 
Berkeley 5, California 


287.—F. E. Fowx, Smithsonian Physical Tables, First Reprint of Eighth Revised 
Edition, The Smithsonian Institution, Washington, D. C., 1934. 


On p. 68, in Table 35, the following corrections should be made: Jo(.62) should 
read .906184 instead of .905184; and Jo(1.89) should read .287631 instead of 
.286631. 


CHARLES R. SEXTON 


288.—J. W. Heap & W. P. Witson, Laguerre Functions: Tables and Properties, 
Monograph No. 183 R, The Inst. of Elec. Engineers, London, 1956. See R 47, 
MT AC, v. 12, 1958, p. 89. 


A partial examination of Table 1, which tabulates 4D values of Laguerre func- 
tions \,(z) for x = 0(.1)1(.2)3(.5)6(1)10 and n = 0(1)10, has disclosed 14 last- 
figure errors and two errors in sign in tabular entries corresponding to n = 5, 
z = 14and n = 6, x = 1.2. These last should read —.0088 and .0086, respectively 


CHARLES R. Sexton 


289.—Her Magesty’s NavuticaL Atmanac Orrice, Interpolation and Allied 
Tables, Her Majesty’s Stationery Office, London, 1956. 


On page 66, in Gregory’s interpolation formula, at the end, for — --- read 
eon 


H. M. Nautical Almanac Office 
Royal Greenwich Observatory 
Herstmonceux Castle 
Hailsham, Sussex 

England 
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290.—Her Masesty’s Navuticat Atmanac Orrice, Subtabulation, A Companion 
Booklet to Interpolation and Allied Tables, Her Majesty’s Stationery Office, 
London, 1958. 


On page 31, the tabular entry corresponding to a = 0, b = 175, r = 5 should 
read 45 instead of 25. 


H. M. Nautical Almanac Office 
Royal Greenwich Observatory 
Herstmonceux Castle 
Hailsham, Sussex 

England 


291.—Puitie M. Morse, Vibration and Sound, Second Edition, McGraw-Hill 
New York, 1948. 


The following corrections should be made in Table V (p. 444): the terminal 
digits of the 4D values for Jo(6.0), No(6.4), and Ni(6.6) should be decreased by a 
unit; the final digits of the values for Jo(8.0), No(3.0), No(5.2), and J2(5.8) should 
be increased by a unit. For J2(7.4) the approximation should read —0.2490 instead 
of —0.2487. 

All of these errors are present in the corresponding table appearing on p. 333 
of the first edition (1936). 


CHARLES R. Sexton 


292.— Cours Rosin, Tables Numériques des Functions Associées de Legendre, Edi- 
tions de la Revue d’Optique, Paris, 1959. 


This table contains a number of entries which have defects due to printing 
difficulties. The publisher has made available a listing containing these errors. 
Copies may be obtained by writing to the author. 


Louis Rosin 


Centre National d’Etudes des Telecommunications 
3, Avenue de la République 

Issy-les-Moulineaux (Seine) 

France 


293.—I. M. Rysurx & I. 8. Grapstern, Summen-Produkt-und Integral-Tafeln: 
Tables of Series, Products, and Integrals, Deutscher Verlag der Wissenschaften, 
Berlin, 1957. 


An examination of these tables has revealed the following errors, some of which 
can be traced to sources cited in the authors’ list of references. 

On p. 416 the final decimal shown in the 16D approximation to Euler’s constant 
should read 9 (when rounded) instead of 5. 








TABLE ERRATA 






































Page Formula For Read 
2 0.126 k 1 k=0 
19 19 
0.131 A,= 30 Ay= in 
7 0.234 2 k=0 k=1 
0.234 3 k=1 k=0 
(2k + 1) (2k + 1)? 
0.234 4 (2k — 1)? (2k — 1)3 
7x4 7x* 
24 1.216 2 +? a ale 
34 1.393 1 2k 2k 
1.393 2 m” i 
44 1.621 4 xz — nz z— nr 
149 3.235 1 1.171 953 619 4 1.171 953 619 3 
186 3.622 7 —1.171 953 619 35 —1.171 953 619 34 
3.622 8 —0.157 660 149 15 —0.157 660 149 17 
274 6.129 1 Wt V2 
T T 
6.129 3 7 Ta 
301 6.338 5 Te, nm, 
303 6.342 2 1 — ¢(2n + 1) ¢(2n + 1) 
330 6.514 8 Z pact (z) Z544(ciz) 
6.514 9 
Page Section For Read 
854 513 854 513 
413 8.21 Bo = 123 Ba = 138 
2 577 867 858 367 2 577 687 858 367 
By = 6 Bus aed 
6 
8.3 (11) = 1.000 494 183 6 ¢(11) = 1.000 494 188 6 
414 8.41 0.154 398 101 8 0.154 981 017 1 
0.149 446 010 5 0.149 445 980 8 
0.144 479 346 3 0.144 464 448 1 
67 108 684 67 108 864 
8.42 0.007 312 527 5 0.007 312 525 9 
0.006 447 210 5 0.006 447 210 3 
0.005 740 026 5 0.005 740 037 7 
0.004 660 148 3 0.004 660 143 5 
415 8.44 1125 1152 
0.000 229 601 1 0.000 229 601 5 
0.000 160 694 8 6 


0.000 160 694 





With the exception of the errors found on pages 7 and 44 and the error of trans- 
position of digits detected on page 415, all the errors noted above appear also in 


the Russian third edition. 


The errors noted in formulas 0.131 and 1.216 2 appear in both Adams [1] and 
Jolley [2], whereas the error noted in formulas 1.393 1 and 1.393 2 occur only in 
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the latter. The errors in formulas 6.129 1 and 6.129 3 apparently were reproduced 
from Magnus and Oberhettinger [3]. 

The terminal-digit errors occurring in definite integrals 3.235 1, 3.622 7, and 
3.622 8 are due to unjustified retention of guard figures from data in Lindman (4], 
Table 113 (p. 61). 

The tabular errata noted on p. 414-415 appear in all the previous editions of 
this book. The exact values of (2m — 1)!!/(2n)!! and of (2n — 1)!!/(2n + 2)! 
for n = 1(1)15 and n = 1(1)14 were given by Lambert [5]. I have verified by 
independent calculation that all of Lambert’s values are free from error. 

Confused notation renders formula 6.362 (p. 307) incorrect. It should read 


n—l 


1 1 1 
C= Bo g-be+5 + - met oe 
1 B., 1 Bo, 0 
-soet cote ra) wr (0 <0 < I. 
J.W.W. 


E. P. Apams, Smithsonian Mathematical Formulae and Tables of Elliptic Functions, 

The Smithsonian Institution, Washington, 1947. 

2. L. B. W. Jouuey, Summation of Series, Chapman & Hall, Ltd., London, 1925. 

3. W. Maonus & F. OBERHETTINGER, Formulas and Theorems for the Functions of Mathe- 
matical Physics, Chelsea Publishing Company, New York, 1954. 

4. C. F. Linpman, Examen des Nouvelles Tables d’ Intégrales Définies de M. Bierens de 
Haan, Stockholm, 1891. Reprinted by G. E. Stechert & Co., New York, 1944. 

&. J. H. LAMBERT, Zusdtze zu den Logarithmischen und Trigonometrischen Tabellen zur 
Erleichterung und Abkitirzung der bei Anwendung der Mathematik vorfallenden Berechnungen, 
Berlin, 1770. 


294.—L. SitBeRSTEIN, Synopsis of Applicable Mathematics with Tables, Bell, 
London, also Van Nostrand, New York, 1923. 


The following corrections should be made in the tables of Bessel functions, 
p. 143. 























Jo (x) 
’ eae ers pee yy pep | 
| For | Read 
.62 .90518 90618 
1.89 . 28663 28763 
2.70 | — .11424 | — .14245 
5.90 .11203 | . 12203 
6.92 | 29873 | . 29874 
| 
Ji (x) 
5.87 — .30109 — .30019 
6.32 — .20291 — .20292 
7.87 21401 21407 





CHARLES R. SEXTON 


























404 TABLE ERRATA 
295.—G. W. Spencetey & R. M. Spencetzy, Smithsonian Elliptic Functions 
Tables, The Smithsonian Institution, Washington, D. C., 1947. 
The following two corrections should be made in the tabular values of ¢ ap- 
pearing on page 187: 
6 r For | Read 
47° 36 0.72012 80371 23 0.72012 80370 23 A 
47° 44 0.82824 23463 76 0.86314 89314 15 
E 
CuHarLes R. SExTon E 
E 
C 
[ 
F 
F 
I 
C 
C 
C 
C 
] 
E 
I 
k 
k 
k 
L 
L 
L 
L 
L 
L 
L 


$< _- —————— _- 














Author 
Ascuer, Marcia 


Bauer, F. L. 
Boersma, J. 


VOLUME XIV 
AUTHOR INDEX 


Papers and Technical Notes 
Short papers and notes are marked (N) in this index 

Title Page 

Explicit solutions of the one-dimensional heat equation for a 
ee ee Pee ns ee: ie ae pre 346 
The formula-controlled logical computer “Stanislaus”... ... 64 
Computation of Fresnel integrals. .......................... 380 
On the factors of certain Mersenne numbers............... . 365 


BRILLHART, JOHN & 
Jounson, G. D. 
Comér, Stie¢ 


Dempsey, E. & Benson, 
G. C. 


Fertis, Henry E. 


Frievper, Dante. C. 
FLoyp, Rosert W. 
Gruuis, Josern & 
Weiss, GEORGE 
GLODEN, A. 
Go.LpBErRG, MICHAEL 
GREENBERG, IRWIN 
HAMMER, Preston C. & 
Wickxe, Howarp H. 


Harumi, Kasasuro, 


KatsurRA, SHIGETOSHI, 


& Wrencu, Joun W., 
JR. 


IsEMANGER, K. R. 
KILLGROVE, 

Raymonp B. 
KIRKPATRICK, 

Epwarp T. 
Kravitz, SIDNEY 
Lane, H.. A. & 

Stevens, D. F. 
LexuMan, R. SHERMAN 
Lino, CutH-Bine 


Loneman, I. M. 
LonemaN, I. M 


Lorn, R. D. 
Lowan, ARNOLD, N. 


Lowan, ARNOLD N. 


Improved methods to calculate the characters of the sym- 


a EF eer rerreeyT terres: seers oe 104 
Asymptotic expansion of the modified Bessel function of the 
ee Te a ee 362 


Note on [ ed = 1 > \rdz EA eters ontario eh 372 
0 g gE 


Tabulation of coefficients for operations on Taylor series..:. 339 
A note on rational approximation....... ; Sr 
Products of Laguerre polynomials... ERA: 60 
oe te eee ee 278 
Rotors in polygons and polyhedra. ... 229 
The evaluation of certain probability integrals. 376 


Quadrature formulas involving derivatives of the integrand.. 3 


2 6° /si 
Values of — [ (“Ya 
™ Jo t 


Res BOARD 6 600k sk ccicexsbnapekeod 379 
The complete factorization of 2'** + 1 (N). ere” 73 
A note on the nonexistence of certain projective planes of 
EES FOE EOE LE CELE LE 70 
Tables of values of the modified Mathieu functions.......... 118 
The congruence 27-1 = 1 (mod p*) for p < 100,000 (N).. 378 
On the evaluation of certain complex elliptic integrals (N)... 195 
On Liouville’s function. . .... 311 
Evaluation at half periods of Weierstrass’ elliptic funetion with 
rectangular primitive period-parallelogram (N)...... 67 
A method for the numerical evaluation of finite integrals of 
oscillatory functions............... 53 
On the utility of Newton’s method for computing complex 
roots of equations (N)........... 187 
Integrals of products of Laguerre poly nomials (N) 375 


On the propagation of round-off errors in the numerical in- 
tegration of the heat equation. ... 139 

On the propagation of round-off errors in the numerical treat- 
ment of the wave equation 


405 


223 











406 INDEX OF REVIEWS BY AUTHOR OF WORK REVIEWED 


Author 


Lowan, ARNOLD N. 
LowaN, ARNOLD N. 


Macuta, D. A. 
Miter, J. C. P. 


Miter, J. C. P. 
Miter, J. C. P. 


Nixoual, Pau J. 
Paty, Gorvon & 
Se1pEen, EstHer 
Peru, I. E. & 
GARRETT, JAMEs R. 
Pops, Davin A. 


RaBINOWITz, PHILIP 
Rapok, J. R. M. 
Rose, Miuton E. 


RorTeNnBERG, A. 
Satzper, HERBERT E. 


Title Page 


On the numerical treatment of heat conduction problems with 
mixed boundary conditions (N)........................... 266 
On the propagation of errors in the inversion of certain tri- 
EE RT. Se ee 333 
A conjugate factor method for the solution of a : cuble (N)... 281 
Numerical quadrature over a rectangular domain in two or 


more dimensions, Part I............ 13 
Numerical quadrature over a rectangular domain in two or 
more dimensions, Part II......................... . 130 
Numerical quadrature over a rectangular domain in two or 
more dimensions, Part III........... Decteiee 0 dri, wales .. 240 
Permanents of incidence matrices (N)...... Bs +biss eA 262 
A problem in abelian groups, with application to the transposi- 
tion of a matrix on an electronic computer (N)............ 189 
High precision calculation of aresin xz, arccos x, and arctan 2. 
Beer eee oe: ees . 270 
A method of “Alternating Corrections’’ for the numerical 
solution of two-point boundary value problems... . . 354 


Abscissas and weights for lobatto quadrature of high order. . 47 
Transcendental equation for the Schrédinger equation (N)... 276 
A method for calculating solutions of parabolic equations 
with a free boundary.................... «sty ae 
The calculation of Toroidal harmonics (N).......... : . 274 
Alternative formulas for osculatory and hy perosculatory j in- 


SN nice sre na nba vine iene rdscees ckiveies 257 
SauzeR, HERBERT E. A note on the solution of quartic equations (N)...... . 279 

Scuerrier, D. & On the numerical evaluation of the 18th perfect number 
ONDREJEA, R. Ma. 8 oe ees 199 
SHanks, DaNTEL A note | on Gaussian twin primes (N)............. 201 

SHanxks, DANIEL On the conjecture of Hardy & Littlewood concerning ‘the 
number of primes of the form n? + a...................... 321 
Srrovup, A. H. Numerical integration formulas of degree two 21 

StTRUBLE, GEORGE Tables for use in quadrature formulas involving derivatives 
Of the mnteewand................. 8 

Swirt, J. D. Construction of Galois fields of characteristic two and ir- 
reducible polynomials.................... . 9 

Watmstey, Cuartes & The Fourier transform of [ce*# + (x — x 1)%}1 arising oo 
Grant, Artur S. G. study of tuned circuit spectra (N).. tual oo 
Watton, Tuomas 8S. & Calculation of the transient motion of subenerged cables Sees 

Po.LAcHEK, HARRY 

Wrenca, Joun W., Jr. Further evaluation of Khintchine’s constant (N)......... .. 370 

Wynn, P. The rational approximation of functions which are formally 
defined by a power series expansion....................... 147 

INDEX OF REVIEWS BY AUTHOR OF WORK REVIEWED 
Review Table 

Author Number Classification Page 
AGaRWALA, S. P. 60 ({K] 298 
Aut, Franz L. 25 [Z] 91 
ARBENZ, Kurt 68 [X] 307 
Armour RESEARCH FouNDATION 34 [P, W] 210 
Arrow, KENNETH J. 21 [Ww] 89 








BHARAQADMOADMHH O0OdOWtWwWwwww wm PS 


— 
LJ 


ORO hwheohohohohohes| 


1 


1 
1 


3 


10 


47 


ge 
98 
91 
07 
10 
89 








INDEX OF REVIEWS BY AUTHOR OF WORK REVIEWED 


Author 


AsubBy, Vau J. 
BaRTHOLOMEW, G. A. 
Barton, D. E. 
BERKSON, JOSEPH 
Bernprt, G. D. 
BHATTACHARYA, NIKHILESH 
BrrkHorr, GARRETT 
Brvins, R. 
Bross, I. D. J. 
CasHWELL, E. D. 
Catron, Henry C. 
Centre Nationa D’Erupes pes T&éL£- 
COMMUNICATIONS 
CuHakKRABorty, P. N. 
CuEeRNOFF, HERMAN 
1IsTova, E. A. 
CGHisrova, E. A. 
CiarkK, F. E. 
Ciatwortny, W. H. 
CULBERTSON, James T. 
Davin, F. N. 
Dennis, JACK BONNELL 
DEPARTMENT OF THE ARMY 
Drxon, W. J. 
Donee, H. F. 
Doornsos, R. 
Drang, C. J., JR. 
DuNNETT, CHARLES W. 
Dunnett, C. W. 
EMERSLEBEN, OTTo 
Epstein, L. Ivan 
Everett, C. J. 
Frereuson, Rosert O. 
Freier, E. C. 
Foster, F. G. 
FREUDENTHAL, H. 
FREVEL, L. K. 
GANTMACHER, F. R. 
GARDNER, MARTIN 
Gawiik, H. J. 
GeRBEs, W. W. 
GERHAUSER, JOHN M. 
GLASSNER, ALVIN 
GLODEN, A. 
Goan, R. L. 
Gooner, J. N. 
GrapsTEIN, I. 8. 
GREENBERG, B. G. 
JREENBERG, B. G. 
Hart ey, H. O. 
HarTree, D. R. 
HeEnsMAN, R. 
Hersuey, A. V. 
Hetz, W. 


Review 
Number 
35 
17 
73 
44 
74 
85 
94 
71 
45 
62 
35 
63 


60 
75 

7 
32 
46 
47 
26 
73 
97 


77 


36 
42 
91 
23 
69 
51 
12 


78 


Table 


Classification 


[S} 
[Ss] 
(K] 
[K] 
[K] 
[K] 
[Q, 8] 
{G] 
(K] 
[K, X, Z] 
[Ss] 
{L] 


(K} 
(K] 

[L] 

(L] 
(K] 
(K] 
(Z] 
(K] 

(Z) 

(2) 

(K] 
(K} 
(K] 
{L, MJ 
(K] 
(K] 

{L) 

[X] 
(K, X, Z] 
[W, X] 
(K] 
[K) 

(Z) 

[D, L] 
(G, X] 
(Z) 

[H, M} 
[L, MJ 
[S, T] 
(T} 

(F] 

(K) 
[X] 
[A, J, L, M] 
(K) 
(K) 
(K] 

(s] 

[L, M] 
[L, M] 
(K) 


P. 


SERSSRRZeE s 


388 85 


SBR eF Fess 


e8ue 


212 


211 


389 








408 INDEX OF REVIEWS BY AUTHOR OF WORK REVIEWED 


Author 


Hiees, L. A. 
Hopeg, P. G., Jr. 
Hogs, M. R. 
Horn, JaKkos 
HorsnE.t, G. 
Howe tt, K. M. 
Hourtson, A. 
JEENEL, JOACHIM 
Jenkins, D. P. 
Jounson, N. L. 
Kamat, A. R. 
Kanrorovicu, L. V. 
KARLIN, SAMUEL 
Karmazina, L. N. 
Karpov, K. A. 
Kasten, E. L. 
Katsura, 8. 
Kimsro, GENEVIEVE M. 
Kurncer, H. 
Krytov, V. I. 
Kou.ipack, SOLOMON 
Lama, R. A. 
Lanepon, Lyte R. 
Lancer, R. E., Editor 
Lancer, Rupotes E., Editor 
Les, Tsai-Hwa 
Lexumann, E. L. 
Levens, ALEXANDER 8S. 
LIEBERMAN, G. J. 
Lowan, A. N. 
LowE.., Herman H. 
McCracken, DaniEt D. 
MatseEn, F. A. 
MaxFIELp, J. E. 
Merropo.uis, N. 
MILLER, JAMES 
Mong, D. T. 
Moorg, P. G. 
Moran, P. A. P. 
Mosgs, Lincoun E. 
NBS, Appiiep MatTHematics SERIES 
No. 50 
NIEVERGELT, E. 
Nisan, G. J. 
Nomura, Y. 
NumericaL Computation BurEau 
Owen, D. B. 
Owen, D. B. 
Ozawa, Masaru 
Paxton, F. A. 
Pearson, E. S. 
Petersen, D. R. 
PRraNGE, EUGENE 
Prins, H. J. 


Review 
Number 


17 


93 


Table 


Classification 


[S] 
[X] 
{L, M] 
[M, X] 
(K] 
(G] 
[K] 
(Z] 

{L, M] 
(K] 
[K] 
[X] 
(W) 
{L] 
{L, M] 
(K] 
({L, M, 8] 
tT) 

(K] 
(X] 
(K] 
(K] 
[X, Z] 
(Q, 8] 
[X] 
[W, Z] 
(K] 
[X] 
(K] 

(T] 

[L] 
[W, Z] 
[s, T] 
{L, M] 
[G] 

[S, T] 
(K] 
(K] 
[K, P] 
[K] 
(K] 


{K] 

[B, C, K, 8] 
[L, M, 8] 
[L] 


=S8skee8e 


w 
R=) 
o 


295 


4a 


81 
295 


391 


291 


16 

















RAMACHANDRAN, K. V. 
Reyno3ps, G. E. 
REYNOLDs, GeorcE E. 
Rosin, Louis 
Roun, J. E. 
Romie, H. G. 

Ross, M. E. 
RorTensBere, M. 
Rysurk, I. M. 
Sauzer, Herpert E. 
SarGENT, LAUREN F. 
Sarwan, A. E. 
Sarnan, A. E. 
Satue, Y. 8. 

Scarr, HERBERT 
Sexar, C. C. 
SELFRIDGE, JoHN L. 
SEe.rrinpcg, R. G. 
Senoupta, J. M. 
SHEerman, B. 
Srorani, Minoru 
SroTani, Minoru 
SoMERVILLE, P. N. 
Steck, G. P. 
Tacout1, G. 

Tate, R. F. 
THORNE, CHARLES J. 
Turtey, J. W. 

vAN LiesHovtT, R. 
Wapstra, A. H. 
Weiss, Haroup 
Westcott, C. H. 


INDEX OF TABLE ERRATA 


Wiuson, WatTER L., Jr. 


Witticu, Hans 
Wooten, J. K., Jr. 


No. Author 
274 Arnot, L. 


286 Burrav, Car. 
275 Canen, E. 

284 Cams, E. 

276 CHEBYSHEV, P. L. 
277 Date, J. B. 

278 Dwient, H. B. 


287 Fow1sz, F. E. 


TABLE ERRATA 


Recherches sur le calcul des forces perturbatrices dans la théorie 
des perturbations séculaires 


FPSSSRe- ws lSSRQARSRESSRSESASNSLPSSERESVER 


[K) 
[L, M] 
[B] 
{L} 
{L} 
[K] 
[Ss] 
[G] 
IA, J, L, M) 
(1) 
[W, X} 
(K] 
(K] 
(K) 
[Ww] 
[K] 
[G] 
[L, M} 
[K] 
[K]} 
(K] 
[K} 
[K] 
[K] 
(K] 
[K] 
[P} 
[D, L) 
[B, C, K, 8) 
[B, C, K, 8) 
[W, Z) 
{s} 
[X] 
[M, X] 
IG] 


Tafeln der Funktionen Cosinus und Sinus, etc. 


Théorie des Nombres 


Eleven and Fifteen-Place Tables of Bessel Fuctions of the First 
Kind, to all Significant Orders 


Teoria delle Congruenze 


Five-Figure Tables of Mathematical Functions comprising Tables of 
Logarithms, Powers of Numbers, Trigonometric, Elliptic, and 
other Transcendental Functions 

“Table of Bessel functions and derivatives J: , J:', J2', Nz, Ni’, 


ye 
N,’” 


Smithsonian Physical Tables, 1st reprint, 8th rev. ed., 1934 


s 


BSSuennseSBSSSSSSSESBPePSSSSBRunBSSSank® 


Page 
218 
400 
218 
308 
218 
219 


219 








410 


Author 


Heap, J.W. & 
Witson, W. P. 
H. M. Navticau 
ALMANAC Or- 
FICE 
290 H. M. Nautica 
ALMANAC OFr- 
FICE 
273 JAHNKE, E. & 
Empeg, F. 
276 Kuurkx, J. P. 
285 MEYER ZUR 
CaPELLEN, W. 
291 Mors, Pxri.ip 
M. 
279 Ritey, J. A. & 
BILinGs, C. 
292 Rosin, Louis 
293 Rysuix, 1. M. & 
GRADSTEIN, 
I. 8. 
Sauzer, Her- 
BERT E. 
SILBERSTEIN, L. 
SPENCELEY, G. 
W. & SPENCE- 
LEY, R. M. 
295 SpENcELEY & 
SPENCELEY 
283 SPENCELEY, 
SPENCELEY & 
EPPERSON 
282 Watson, G. N. 


BS 


8 


2 


88 





CORRIGENDA 


Title 
Laguerre Functions: Tables and Properties 


Interpolation and Allied Tables 

Subtabulation, A Companion Booklet to Interpolation and Allied 
Tables 

Tables of Functions with Formulae and Curves 

“Uber die Tafel primitiver Wurzeln” 

Integraltafeln. Sammlung unbestimmter Integrale elementarer 
Funktionen 

Vibration and Sound 

‘‘Gaussian quadrature of some integrals involving Airy functions”’ 

Tables Numériques des Functions Associées de Legendre 


Summen-Produkt-und Integral-Tafeln: Tables of Series, Products, 
and Integrals 


“Orthogonal polynomials arising in the numerical evaluation of 
inverse Laplace transforms’’ 

Synopsis of Applicable Mathematics with Tables 

Smithsonian Elliptic Functions Tables, 1947 

Smithsonian Elliptic Functions Tables, 1947 


Smithsonian Logarithmic Tables to Base e and Base 10, 1952 


A Treatise on the Theory of Bessel Functions 


NOTES 
Title 


Mathematics of Computation 


New Journals 


Automatic Programming of Digital Computers—National Information Centre 


Author 
FrevvER, D. C. 
HensMaN, R. & 

JENKINS, D. P. 
Her Magessty’s Nav- 
TICAL ALMANAC OrF- 
FICE 
Line, Cu1H-BING 
Parker, E. T. & 
NrKo.al, Pavut J. 
Rosin, L. 


CORRIGENDA 


Title 


*‘A note on summation formulas of powers of roots’’ 


) 
Tables of e af eat for Complex z 
T x 


Subtabulation, A Companion Booklet to Interpolation and Allied 
Tables 


“Tables of values of 16 integrals of algebraic-hyperbolic type”’ 
‘A search for analogues of the Mathieu groups”’ 


Table Erratum, 261, MTAC, 1958 


Page 


401 


219 


401 


308 


221 


Page 


309 
309 


Page 
97 


310 


97 


222 
222 


222 








— 


01 


97 


310 


222 
222 


222 








REVIEWS BY SUBJECT 411 


REVIEWS BY SUB)ECT 


A. Arithmetical Tables, Mathematical Constants 


Author /s 
I. M. Rysurx & 
I. S. GRADSTEIN 


Georce E. ReYNoLps 


A. H. Wapstra, 
G. J. Nueu & 
R. van LiEsHOUT 


A. H. Wapstra, 
G. J. Nueu & 
R. van LiEsHOUT 


L. K. FREevEt, 
J. W. Turtey & 
D. R. PETERSEN 


A. GLODEN 


F. R. GANTMACHER 
K. M. Howe. 
EvuGENE PRANGE 


M. RoTeNBERG, 
R. Brvins, 
N. Merropouts, & 
J. K. Wooten, Jr. 
Joun L. SELFRIDGE 


H. J. GAWLIK 


A. N. Lowan 


Hersert E. Sauzer & 
GENEVIEVE M. 
KIMBRO 


Review 
Number Title Page 
69 Summen-, Produkt- und Integral-Tafeln: Tables 381 
of Series, Products, and Integrals 
B. Powers 


29 A New Method of Cube Root Extraction on Desk 204 
Calculators 
1 Nuclear Spectroscopy Tables 75 


C. Logarithms 
1 Nuclear Spectroscopy Tables 75 


D. Circular Functions 


2 Seven-Place Table of Iterated Sine 76 


F. Theory of Numbers 


42 Table des Solutions de la Congruence z*+1=0 284 
(mod p) pour 800 000 < p < 1 000 000 


G. Higher Algebra 


43 Applications of the Theory of Matrices 284 

3 Revised Tables of 6j-Symbols 76 

70 An Algorism for Factoring X" — 1 over a Finite 382 
Field 

71 The 3 — j and 6 — j Symbols 382 

30 On Finite Semigroups 204 


H. Numerical Solutions of Equations 


4 Zeros of Legendre Polynomials of orders 2-64 and 77 
weight coefficients of Gauss quadrature formulae 


I. Finite Differences, Interpolation 


5 The Operator Approach to Problems of Stability 78 
and Convergence of Solutions of Difference Equa- 
tions and the Convergence of Various Iteration 
Procedures 

Table for Bivariate Osculatory Interpolation Over a 384 
Cartesian Grid 


“1 
i) 








412 


Author/s 


I. M. Rysuix, & 
I. S. GRADSTEIN 


D. E. Barton & 
F. N. Davip 
JosEPH BERKSON 


G. D. Bernpt 

I. D. J. Bross & 
E. L. Kasten 

E. D. CasHwe.it & 

C. J. Everett 
Herman CHERNOFF & 
Lincotn E. Moses 

F. E. Cuark 
W. H. Ciatwortuy 


W. J. Drxon 


H. F. Dopgse & 

H. G. Romie 
R. Doornzos & 

H. J. Prins 
Cuar.es W. DuNNETT 


C. W. Dunnetr & 
R. A. Lamu 

E. C. Freier, 
H. O. Hartiey & 
E. S. Pearson 

F. G. Foster 


W. Herz & H. Kuincer 


G. HorsNneLu 
A. Hurtson 


N. L. Jounson 

Sotomon KvULLBACK 

E. L. LEaMANN 

G. J. LizBpeRMAN 

P. G. Moore 

P. A. P. Moran 

NBS, Appirep Mats. 
Serres, No. 50 

E. NIEVERGELT 

D. B. OwEn 


D. B. Owren & 

D. T. Monx 
K. V. RAMACHANDRAN 
A. E. Saruan & 

B. G. GREENBERG 


Review 
Number 
J. Summation of Series 


69 


73 


74 


62 


75 


47 


48 


49 


76 


51 


77 


78 


52 
79 


31 
81 


92 


82 


57 


& & 





REVIEWS BY SUBJECT 


Title 


Summen-, Produkt- und Integral-Tafeln: Tables 
of Series, Products, and Integrals 
K. Statistics 


“‘A test for birth order effect’’ 


“Tables for use in estimating the normal dis- 
tribution by normit analysis” 

‘Power functions of the gamma distribution”’ 

‘Rapid analysis of 2 x 2 tables’’ 


A Practical Manual on the Monte Carlo Method for 
Random Walk Problems 
Elementary Decision Theory 


“Truncation to meet requirements on means’’ 

Contributions on Partially Balanced Incomplete 
Block Designs with Two Associate Classes 

‘Estimates of the mean and standard deviation of 
a normal population’”’ 

Sampling Inspection Tables 


“On slippage tests. I’’ 


‘*Tables of the bivariate normal distribution with 
correlation 1/+/2” 

“Some tables of the multivariate normal prob- 
ability integral with correlation coefficients 1/3” 

‘Tests for rank correlation coefficients. I” 


‘‘Upper percentage points of the generalized beta 
distribution. ITI” 

“Untersuchungen zur Frage der Verteilung von 
Objekten auf Platze”’ 

‘Economical acceptance sampling schemes”’ 

‘Further critical values for the sum of two var- 
iances”’ 

“Optimal sampling for quota fulfillment’’ 

Information Theory and Statistics 

Testing Statistical Hypotheses 

‘Tables for one-sided statistical tolerance limits”’ 

“The two-sample t-test based on range’”’ 

The Theory of Storage 

Tables of the Bivariate Normal Distribution Func- 
tion and Related Functions 

“Die Rangkorrelation U”’ 

Tables of Factors for One-sided Tolerance Limits 
for a Normal Distribution 

Tables of the Normal Probability Integral 


“On the Studentized smallest chi-square”’ 

‘Tables for best linear estimates by order statis- 
tics of parameters of single exponential distribu- 
tions from singly and doubly censored samples’’. 


R28 8 88 es 


B & 


g 


295 


389 
297 








ge 


81 


S28 $$ 8& & FE 


B & 


z 


295 


389 
297 








Author/s 
A. E. Saruan & 
B. G. GREENBERG 


Y. 8S. Satne & 
A. R. Kamat 


C. C. SEKAR, 
S. P. Acarwata & 
P. N. CHAKRABORTY 

J. M. Sencupta & 
NIKHILESH 
BHATTACHARYA 

B. SHERMAN 

Minoru SIoTANI 


Minoru SroTranit & 
Masaru Ozawa 


P. N. SoMERVILLE 
G. P. Steck 
G. Taevut! 


R. F. Tate & 
R. L. Goan 

A. H. Wapstra, 
G. J. Niuscu & 
R. van LiEsHOUT 


CENTRE NATIONAL 
D’ErupEs DEs 
T£ELECOMMUNICATIONS 

E. A. GBrstova 


Orro EMERSLEBEN 


L. K. FRevEt, 
J. W. Turtey & 
D. R. PETERSEN 
W. W. GERBEs, 
G. E. ReyNo.ps, 
M. R. Hogs & 
C. J. Drang, Jr. 
R. Hensman & 
D. P. JENKINS 


A. V. HersHEey 


L. N. Karmazina & 
f. A. CHistova 


K. A. Karpov 


REVIEWS BY SUBJECT 413 


Review 
Number Title Page 
84 “Estimation of location and scale parameters by 389 


order statistics from singly and doubly censored 
samples. Part II’’ 


59 “Approximations to the distributions of some 297 
measures of dispersion based on successive dif- 
ferences” 

60 **On the power function of a test of significance for 298 
the difference between two proportions”’ 

85 “Tables of random normal deviates” 390 

61 “Percentiles of the w, statistic’”’ 299 

86 ‘‘Note onthe utilization of the generalizedStudent 390 
ratio in analysis of variance or dispersion”’ 

87 ‘Tables for testing the homogeneity of kindepend- 391 


ent binomial experiments on a certain event 
based on the range’”’ 


88 “Tables for obtaining non-parametric tolerance 391 
limits’’ 

89 “A table for computing trivariate normal prob- 391 
abilities” 

90 “Tables of tolerance coefficients for normal popu- 392 
latons”’ 

91 ‘Minimum variance unbiased estimation for the 393 
truncated Poisson distribution”’ 

1 Nuclear Spectroscopy Tables 75 


L. Higher Mathematical Functions 


63 Tables numériques des fonctions associées de te- 300 
gendre. 


~~ 


Tables of Bessel functions of real argument and of 79 
integrals involving them 


8 ‘“‘Werte einer Zetafunktion zweiter Ordnung mit 80 
Argument s = 2” 

2 Seven-Place Table of Iterated Sine 76 

11 Table of S(x) and its First Eleven Derivatives 82 


os 
12 Tables of a af edt for Complex z 83 
Vr J, 

13 Computing Programs for the Complex Exponential 83 
Integral 

32 Tables of Bessel functions of imaginary argument 208 
and of integrals involving them 

14 Tables of the function F(z) = [ e** dz in the com- 84 


plex domain 








414 


Author/s 
Herman H. Lowe. 


Y. Nomura & 

S. Katsura 
NumeEricaL Computa- 
TION Bureau, RE- 

Port No. 11 
F. A. Paxton & 
J. E. Roun 
Louis Rosin 


I. M. Rysaix & 
I. S. GRADSTEIN 
R. G. SELFRIDGE & 
J. E. MaxFretp 


H. J. Gaw1ikx 


W. W. GERBEs, 
G. E. REyYNo.ps, 
M. R. Hogs & 
C. J. Drang, Jr. 
R. Hensman & 
D. P. JENKINS 
A. V. HersHey 


Jaxos Horn & 
Hans WIttTIcH 


K. A. Karpov 


Y. Nomura & 
S. Katsura 
I. M. Rysarx & 
I. S. GRADSTEIN 
R. G. SELFRIDGE & 
J. E. MAxFIeLp 


Armour RESEARCH 
FouNDATION 
CHARLES J. THORNE 


P. A. P. Moran 
Garrett Birkuorr & 


R. E. Lancer, 
Editors 


Vat J. Asupy & 
Henry C. Catron 





REVIEWS BY SUBJECT 


Review 
Number 


9 


15 


10 


69 


ll 


12 
13 


93 


15 


69 


34 


16 


92 


94 


Title 

Tables of the Bessel-Kelvin Functions Ber, Bei, Ker, 
Kei, and their Derivatives for the Argument Range 
0(0.01)107 .50 

“Diffraction of electric waves by circular plate and 
circular hole’’ 

Tables of Whittaker Functions (Wave Functions in 
Coulomb Field) Part 2 


Tables of the Incomplete Elliptic Integrals of the 
First and Third Kind 

Fonctions sphérique de Legendre et fonctions sphé- 
roidales 

Summen-, Produkt- und Iniegral-Tafeln: Tables 
of Series, Products, and Integrals 

A Table of the Incomplete Elliptic Integral of the 
Third Kind 
M. Integrals 


Zeros of Legendre Polynomials of orders 2-64 and 
weight coefficients of Gauss quadrature formulae 
Table of S(x) and its First Eleven Derivatives 


2 x 
Tables ST e" | e~*? dt for Complex z 
Lal z 


Computing Programs for the Complex Exponential 
Integral 
Gewodhnliche Differentialgleichungen 


z 
Tables of the function F(z) = e*” dz in the com- 
0 
plex domain 
‘Diffraction of electric waves by circular plate and 
circular hole”’ 
Summen-, Produkt- und Integral-Tafeln: Tables 
of Series, Products, and Integrals 
A Table of the Incomplete Elliptic Integral of the 
Third Kind 
P. Engineering 
Proceedings of the Fourth Annual Computer Ap- 
plications Symposium, 1957 
Temperature Tables: Part 1. Cne-Layer Plate, One- 


Space Variable, Linear 
The Theory of Storage 


Q. Astronomy 


Proceedings of Symposia in Applied Mathematics, 
Vol. IX, “Orbit Theory”’ 


S. Physics, Geophysics, Crystallography 


35 


Tables of Nuclear Reaction Q Values 


Page 
81 


81 


301 


381 


302 


~J 
~I 


82 


83 


83 


394 


381 


302 


211 








81 


01 


381 


302 


~I 


~I 


83 


83 


302 


210 


85 


393 


395 


211 








G. A. BarTHotomew & 
L. A. Hiees 

GARRETT BirkHorr & 
R. E. Lancer, 
Editors 

D. R. HarTREE 

JAMES MILLER, 
JoHn M. GERHAUSER 
& F. A. MaTsen 

Y. Nomura & 
S. Katsura 

M. E. Rose 

A. H. Wapstra, 
G. J. Nisan & 
R. van LiEsHoUT 

C. H. Westcotrr 


ALVIN GLASSNER 


James MILLER, 
JoHN M. GERHAUSER 
& F. A. MaTsen 


Armour RESEARCH 
FOUNDATION 

KENNETH J. ARROW, 
SAMUEL KARLIN & 
HERBERT SCARF 

Rosert O. Ferauson 
& LaureN F. SARGENT 

Danie, D. McCracken, 
Haroip Weiss «& 
Tsat-Hwa LEE 


REVIEWS BY SUBJECT 


_ 


7 Compilation of Thermal Neutron Capture Gamma 
Rays 


94 Proceedings of Symposia in Applied Mathematics, 
Vol. IX, “Orbit Theory”’ 


66 The Calculation of Atomic Structures 
36 Quantum Chemistry Integrals and Tables 


15 ‘Diffraction of electric waves by circular plate and 
circular hole’’ 
18 Internal Conversion Coefficients 
1 Nuclear Spectroscopy Tables 
19 Effective Cross Section Values for Well-Moderated 


Thermal Reactor Spectra 


T. Chemistry 


20 The Thermochemical Properties of the Orides, 
Fluorides, and Chlorides to 2,500° K 
36 Quantum Chemistry Integrals and Tables 


W. Economics and Social Sciences 


34 Proceedings of the Fourth Annual Computer Ap- 
plications Symposium, 1957 


21 Studies in the Mathematical Theory of Inventory 
and Production 

67 Linear Programming 

22 Programming Business Computers 


X. Numerical Analysis and Applied Mathematics 


Kurt ARBENZ 


E. D. CasHwELL & 
C. J. Everett 
L. Ivan EpsTEeIN 
Rospert O. FerRGuson 
& LaurReEN F. SARGENT 
F. R. GANTMACHER 
J. N. Goopvrer & 
P. G. Hopes, Jr. 
Jaxos Horn & 
Hans WITTICH 
L. V. KaNTOROVICH 
& V. I. Kryiov 
LyLe R. LANGDON 
Rupoupu E. LANGER, 
Editor 
ALEXANDER S. LEVENS 


68 ‘‘Integralgleichungen fiir einige Randwert probleme 
fiir Gebiete mit Ecken”’ 

62 A Practical Manual on the Monte Carlo Method for 
Random Walk Problems 

37 Nomography 

67 Linear Programming 

43 Applications of the Theory of Matrices 

23 Elasticity and Plasticity 

93 Gewohnliche Differentialgleichungen 

24 Approximate Methods of Higher Analysis 

39 Approximating Functions for Digital Computers 

95 Boundary Problems in Differential Equations 

38 Nomography 


415 


g 


304 
211 


87 
75 


211 


210 


89 


305 


89 


394 


90 


214 


395 





416 REVIEWS BY SUBJECT 


Water L. Wison, Jr. 96 ‘Operators for solution of discrete Dirichlet and 
Plateau problems over a regular triangular grid”’ 


Z. Calculating Machines and Mechanical Computations 


Franz L. ALT 25 Electronic Digital Computers 
E. D. CasHwe.t.t & 62 A Practical Manual on the Monte Carlo Method for 
C. J. Everetr Random Walk Problems 
James T. CULBERTSON 26 Mathematics and Logic for Digital Devices 
Jack BoNNELL DENNIS 97 Mathematical Programming and Electrical Net- 
works 
Dept. OF THE ARMY 27 Catalog of Commercially Available Automatic Data 
Processing Systems 
H. FREUDENTHAL 40 ‘‘Logique Mathématique Appliquée”’ 
Martin GARDNER 28 Logic Machines and Diagrams 
JOACHIM JEENEL 41 Programming for Digital Computers 
Lyte R. LANGDON Approximating Functions for Digital Computers 
Danret D. McCracken, 22 Programming Business Computers 
Harotp Weiss & 
Tsat-Hwa LEE 


Compiled by 

Elizabeth Harper 

Applied Mathematics Laboratory 
David Taylor Model Basin 
Washington 7, District of Columbia 








q¢€mga50Amk oe 





—] 


J. 
K. 
L. 
M. 
N. 
0. 
P. 
Q. 
R.z 
8. 
mT. 
U. 
. 
W. 
X. 


rmons so oO PP 


CLASSIFICATION OF TABLES 


Arithmetical Tables, Mathematica! Constants 
Powers 

Logarithms 

Circular Functions 

Hyperbolic and Exponential Functions 
Theory of Numbers 

Higher Algebra 

Numerical Solution of Equations 
Finite Differences, Interpolation 
Summation of Series 

Statistics 

Higher Mathematical Functions 
Integrals 

Interest and Investment 

Actuarial Science 

Engineering 

Astronomy 

Geodesy 

Physics, Geophysics, Crystallography 
Chemistry 

Navigation 

Aerodynamics, Hydrodynamics, Ballistics 


. Economics and Social Sciences 


Numerical Analysis and Applied Mathematics 
Calculating Machines and Mechanical Computation 





Mathematics of Computation 


TABLE OF CONTENTS 


OcroBER 1960 


On Liouville’s Function R. SHerMan LEHMAN 
On the Conjecture of Hardy & Littlewood concerning the Number of Primes 
DaNtEL SHANKS 
On the Propagation of Errors in the Inversion of Certain Tridiagonal 
ArnoLtp N. Lowan 

Tabulation of Coefficients for Operations on Taylor Series 
DanreL C. FIELDER 
Explicit Solutions of the One-dimensional Heat Equation for a Composite 
Wall Marcia ASCHER 
A Method of ‘Alternating Corrections” for the Numerical Solution of Two- 
point Boundary Value Problems Davin A. Pore 


TrecHNICAL Notes AND SHORT PAPERS 


Note on the Asymptotic Expansion of the Modified Bessel Function 
of the Second Kind E. Dempssty & G. C. BENSON 

On the Factors of Certain Mersenne Numbers 
JoHN Brit~tHart & G. D. JoHNSON 
Further Evaluation of Khintchine’s Constant. .Joun W. Wrencu, Jr. 


Note on [ e Jo (“) J; (2) 2" Henry E. Ferris 


Integrals of Products of Laguerre Polynomials 
The Evaluation of Certain Probability Integrals... . Irwin GREENBERG 
The Congruence 2?-! = 1 (mod p’) for p < 100,000. . .Sipney Kravirz 


Values of — -[ (=! 2) dt KasaBurO HarvuMt, 


SHIGETosHI Katsura, & JoHN W. WRENCH, JR. 
Computation ee” J. BormrsMA 
REVIEWS AND DeEscriPTIONS OF TABLES AND BOOKS.................... 
Rysuik & GRADSTEIN 69, PRANGE 70, RoTENBERG, Bivins, METROPOLIS 
& Wooten 71, Sauzer & Kimsro 72, Barton & Davin 73, BERNDT 
74, CHERNOFF & Moses 75, Doornsos & Prins 76, Foster 77, Herz 
& Kuincer 78, Hurrson 79, KuLLBACK 80, ee 81, OWEN 82, 
RAMACHANDRAN 83, SARHAN & GREENBERG 84, SENGUPTA & Buarra- 
CHARYA 85, SIoTANI 86, S1orani1 & OzAwa 87, Sosmavitae 88, STECK 
89, Tacuti 90, Tare & Goan 91, Moran 92, Horn & WirticH 93, 
BrrkHorr & LANGER 94, LANGER 95, WILSON 96, DENNIS 97 
TABLE ERRATA.. 


Burrav 286, Suma 287, Heap & WiLson 288, Her Masesty’s Nav- 
TICAL ALMANAC OFFICE 289, Her Masesty’s NavutTicaL ALM .~ AC 
OFFicE 290, Morse 291, Rosin 292, Rysurk & GRADSTEIN 293, SIL- 
BERSTEIN 294, SPENCELEY & SPENCELEY 295 

InpicEs TO VoLUME XIV Fi ta 
Author Index, Papers and Technical Notes. anasceks 
Reviews by Author of Work Reviewed 
Table Errata... . 


Corrigenda 
Reviews by Subject 


Published Quarterly by the 


National Academy of Sciences-National Research Council 








